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The masses of the particles have been determined by 
measurement of magnetic rigidity and range of tracks 
originating from bombardment with 100-Mev x-rays of the 
glass wall of a cloud chamber. Range was determined by 
observing the stopping in a series of aluminum foils 
mounted in the cloud chamber and the magnetic rigidity 
was determined from measufement of curvature of the 
track before it entered the foils. The magnetic field of 
2650 gauss permits particles of intermediate mass to be 
readily distinguished from protons. About 200 heavily 
ionizing particles were observed in a total of 3100 pictures. 
102 of these tracks stopped either in the gas of the chamber 


or in one of the foils. Photographs of typical tracks are 
presented. None of these 102 tracks could be identified 
as caused by a particle having a mass in the neighborhood 
of, or less than, the meson mass (i.e., 200 electron masses). 
All of these tracks showed positive curvature. The ma- 
jority could definitely be identified as protons. The few 
deviations can be accounted for by considering the effects 
of multiple scattering in the gas of the cloud chamber. 
Several particles were observed to penetrate all of the 
foils; protons with energies in excess of 28 Mev are 
indicated. The energy distribution of the protons is given. 





N previous work,'? it has been shown that 

disintegrations of nuclei with the emission of 
heavy charged particles are induced by 100-Mev 
x-radiation. The present work was undertaken 
to establish more certainly the identity of the 
heavy particles emitted in these disintegrations. 

The Wilson cloud chamber affords a conveni- 
ent means of determining the mass of these 
particles. There are four properties of a particle, 
namely, the range, magnetic rigidity (Hp), dens- 
ity of ionization, and rate of change of Hp with 
range, which can be determined with the cloud 
chamber. As pointed out by Corson and Brode,* 


* Presented in part at the September, 1946 meeting of 
the American Physical Society in New York City. 
aaa C. Baldwin and G. S. Klaiber, Phys. Rev. 70, 259 


*M. Schein, A. J. Hartzler, and G. S. Klaiber, Phys. Rev. 
70, 435 (1946). 
(1938) R. Corson and R. B. Brode, Phys. Rev. 53, 773 


any two of these quantities are sufficient to 
determine the mass of a particle, but of the 
various combinations of these four quantities, 
the use of range and magnetic rigidity affords the 
most accurate determination of the mass. These 
are the two quantities which were used for the 
mass determinations in this work. The ranges 
of the particles were determined by observing 
their stopping in a series of foils mounted in the 
cloud chamber. The magnetic rigidity was deter- 
mined from a measurement of the curvature of 
that part of the track generated by the particle 
before it entered the stopping foils. 

A schematic drawing of the experimental set- 
up is shown in Fig. 1. The x-ray beam from the 
100-Mev betatron‘ is collimated by an aperture 
1} inches high and 1? inches wide in a lead 


*W. F. Westendorp and E. E. Charlton, J. App. Phys. 
16, 581 (1945). J. App 
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Fic. 1. Schematic diagram of experimental arrangement. 


shield having a total thickness of 16 inches. The 
collimated beam of x-rays is then allowed to 
strike one edge of the glass wall of the cloud 
chamber. Heavy particles emitted by the photo- 
disintegrations in the glass can then be observed 
traveling across the chamber, first through an 
unobstructed part of the chamber where curva- 
ture measurements can readily be made, and 
finally inte the series of aluminum foils in which 
they are stopped. 

The cloud chamber used in this work had a 
diameter of 30 cm and a depth of 3 cm of the 
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chamber was illuminated. The source of illumi. 
nation was two xenon capillary lamps mounted 
on opposite sides of the chamber. The aluminum 
foils in the chamber increased in thickness jp 
logarithmic steps from 0.006” for the first foil 
to 0.080” for the last. The chamber was filled 
with air or hydrogen and a mixture of aleohg! 
and water vapors to a pressure of one atmo. 
phere. A pair of Helmholtz coils supplied 4 
magnetic field of 2650 gauss perpendicular tp 
the plane of the chamber. The chamber was 
operated once a minute in synchronism with 
the betatron. 

The alignment of the hole in the lead shield 
with the very narrow 100-Mev x-ray beam was 
accomplished using radiographic techniques, Ap 
8-X10-inch x-ray film wrapped in thin black 
paper was placed against the side of the col 
limator shield facing the betatron and in such 
a position that the aperture of the collimator was 
at the center of the film. A }-inch sheet of lead 
was then placed in front of the film to act as an 
intensifier. Figure 2 shows the radiograph ob. 
tained. Here one can seé the circular spot due to 
the x-ray beam, the outline of the aperture, and 
their position relative to one another. The 
accuracy of the line-up used in this work is 
apparent from this picture. 

The degree of collimation obtained with the 
16 inches of lead used for this purpose can be 


Fic. 2. Radiograph of 100-Mev x-ray beam incident 
on lead shield aperture. 
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Fic. 3. Cloud-chamber photograph taken without magnetic 
field showing degree of collimation of x-ray beam. 


seen in Fig. 3. For this picture, the chamber was 
placed so that the x-ray beam passed through 


the chamber instead of just hitting the edge. 
This picture was taken without a magnetic field. 
Since at high quantum energies, the pair and 
compton electrons produced in the chamber wall 
travel predominantly in a forward direction, they 
give a good indication of the extent of the x-ray 
beam. Figure 4 is similar to Fig. 3 except for a 
magnetic field of 2650 gauss applied perpen- 
dicular to the plane of the chamber. Here one 
can see very clearly the separation of the elec- 
trons and positrons in the electron beam from 
the glass wall. These pictures indicate that quite 
good collimation of the x-ray beam was achieved. 

In selecting the thicknesses of the stopping 
foils used in this experiment, one criterion had 
to be satisfied. The foils had to be of such a 
thickness that for a particle stopping in any 
given foil, the uncertainty in its range, caused 
by the thickness of the foil, must be small enough 
compared to the total range of the particle to 
allow sufficient resolution in its mass determina- 
tion to distinguish between protons and particles 
having masses of the order of the meson mass. 
That the foils used in this experiment satisfy 
this criterion can be seen in the following illus- 
trations. Figure 5 is a diagram illustrating how a 
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Fic. 4. Cloud-chamber photograph, similar to Fig. 3, taken 
with magnetic field of 2650 gauss. 


proton, a meson of mass 200 times the electron 
mass, and a particle of 50 times the electron 
mass should appear in the chamber with a 
magnetic field of 2650 gauss, if they stop in the 
first foil. The shaded area for each particle indi- 
cates the possible variations in the radius of 
curvature of that particle because of the uncer- 
tainty of its range owing to foil thickness. It 
can be seen that for particles stopping in the first 
foil, there should be no uncertainty as to the 
identity of the particle. Figure 6 is a diagram 
similar to Fig. 5 showing the resolution for 
particles stopping in the last foil. Here, as for 
particles stopping in the first foil, the resolution 
is still good. In a similar manner, it can be shown 
that for a particle stopping in any of the other 
foils, there is sufficient resolution in the mass 
determination to identify the particle. 

A general summary of the resolution in the 
mass determination obtained in this experiment 
is shown in Fig. 7. Here each heavy diagonal 
line is a plot of range against magnetic rigidity 
for the particular particle indicated on the curve. 
The portions of the ranges corresponding to the 
various foils are indicated by the horizontal 
shaded areas on the graph. It can be seen from 
this diagram that for particles stopping in any 
foil there is no overlapping in the values of Hp 
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Fic. 5. Diagram illustrating how 


rticles of mass indicated would appear in the cloud chamber with a mag- 


netic field of 2650 gauss, if they stop in the first foil. The extent of the shaded area shows the variation of curvature 


with range due to foil thickness. 


for any of the particles shown. It should also be 
noted that a meson of mass 200 with a range 
corresponding to stopping in the last foil has a 
smaller p than a proton stopping in the first foil. 

A total of 3100 pictures were taken with the 
foil arrangement described above and with a 
magnetic field of 2650 gauss perpendicular to the 
plane of the chamber. There were about 200 
heavily ionizing particles observed in this series 
of pictures. 102 of these observed tracks stopped 
either in the gas of the chamber or in one of the 
foils. Examples of some of these tracks are shown 
in the accompanying figures. 

Figure 8 shows a particle having a range of 
14.5 cm in air and a measured radius of curva- 
ture of 75 cm. These values of range and curva- 
ture are consistent for a particle of protonic 
mass having an energy of 3.1 Mev. A meson of 
mass 200 times the electron mass and having 
this range of 14.5 cm would generate a track 
having a radius of curvature of only 16 cm. 

Figure 9 is of a particle stopping in the first 
foil. The range of this particle is between 13.5 
and 44 cm of air (i.e., the equivalent thickness 


of the foil). The radius of curvature of this track 
is 100 cm. From Fig. 7, it can be seen that thes 
data are consistent for a proton. A 200-mas 
meson with the extreme range of 44 cm would 
produce a track of only 26-cm radius of curva 
ture. This particle is a proton having an energy 
somewhere in the range of 3'to 5.8 Mev. 

In Fig. 10 a particle is shown stopping in the 
second foil. This particle has a range between 
47 and 92 cm of air and its track exhibitsa 
curvature with a radius of 120 cm. These data 
identify this particle to be a proton. The track 
of a meson of mass 200 having a range of 92 cm 
would have a radius of curvature of only 32 cm 
This proton has an energy between 6 and 88 
Mev. 

Figure 11 shows a particle with a range some 
where between 105 and 195 cm of air and whow 
track has a radius of curvature of 200 cm. This 
identifies the particle as a proton with an energy 
between 9.5 and 13.5 Mev. The track of a meson 
of mass 200 with a 195-cm range should havea 
radius of curvature of 40 cm. 

A particle stopping in the fourth foil is show 
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Fic. 6. Diagram illustrating how particles of mass indicated would apeees in the cloud chamber with a mag- 


netic field of 2650 gauss, if they stop in fifth plate. The extent of the sha 


with range due to foil thickness. 


in Fig. 12. The range of this particle is between 
190 and 380 cm of air, and the radius of curvature 
of its track is about 200 cm. A 200-mass meson 
with a range of 380 cm would generate a track 
with only a 49-cm radius of curvature. This 
particle is a proton with an energy between 13.2 
and 19.5 Mev. 

Figure 13 shows a particle stopping in the 
fifth foil and thus having a range between 400 
and 800 cm of air. The curvature of this track 
is so slight that it was impossible to get an 
exact measurement of its radius of curvature 
with the measuring device used for this purpose. 
The best that can be said is that the radius of 
curvature is greater than 250 cm. From Fig. 7 
it can be seen that this particle can only be a 
proton. A meson of mass 200 with a range of 
800 cm would have a track curvature with only 
a 60 cm radius. This proton has an energy 
between 20 and 30 Mev. 

The examples that have been shown are 
typical of the 102 tracks that were observed to 
stop in the chamber. Nong of these 102 tracks 
could be identified as being caused by a particle 


ed area shows the variation of curvature 


having a mass in the neighborhood of, or less 
than, the meson mass (i.e., 200 times the electron 
mass). All of these tracks showed a positive 
curvature. While the majority of these 102 tracks 
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Fic.*7. Diagram summarizing resolution of mass deter- 
mination as_a function of foil thickness (horizontal shaded 
areas). 
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Fic. 8. Range=14.5 cm air; p=75 cm; Hp=2X10 gauss cm; Identity: Proton, Energy =3.1 Mev. 
Fic. 9. Range = 13.5-44 cm air; p= 100 cm; Hp=2.6X 10° gauss cm; Identity: Proton, Energy =3-5.8 


Mev. 
Fic. 10. Range = 47-92 cm air; p= 120 cm; Hp=3.2X 10 gauss cm; Identity: Proton, Energy =6-8.8 


Mev. 
Fic. 11. Range = 105-195 cm air; p= 200 cm; Hp=5.3X 10 gauss cm; Identity: Proton, Energy =9.5- 


13.5 Mev. 


could definitely be identified as caused by pro- 
tons, there were a few cases observed in which the 
mass of the particle, as calculated for its meas- 
ured range and magnetic rigidity, deviated some- 


what from the proton mass. In only two cas 
however, did the mass determination lead to 
value for the mass less than 1500 times # 
electron mass. One example of this is shown! 





RANGE-MOMENTUM MEASUREMENTS 


Fic. 12. Range = 190-380 cm air; p~200 cm; Hp~6.6 X 10° gauss cm; Identity: Proton, Energy = 13.2- 


19.5 Mev. 
Fic. 13. Range = 400-800 cm air; p>250 cm; Hp>6.6X 10° gauss cm; Identity: Proton, Energy = 20- 


30 Mev. 
Fic. 14. Range >14 cm air; p=45 cm; Hp=1.2X10 gauss cm; Mass<1000 m.; Identity: Proton if 


Pocat < 150 cm. 
Fic. 15. Range >750 cm air; Energy (proton) >28 Mev. 


Fig. 14. This picture shows a particle whose 
range is greater than 14 cm of air. Only this 
lower limit can begiven for the range as the 
particle runs out of the lighted portion of the 


chamber before stopping. The radius of curva- 
ture of this track is 45 cm, giving a magnetic 
rigidity for the particle of 1.210° gauss cm. 
Using these values of range and magnetic rigid- 
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Fic. 16. Range distribution of protons. 


ity, one obtains an upper limit for the mass of 
this particle of 1000 times the electron mass. 
While such a value for the mass appears offhand 
to be incompatable with a proton, consideration 
of the effect of multiple scattering by the gas in 
the cloud chamber on the path of a particle 
indicates that the occasional observation of a 
track with abnormally high curvature is to be 
expected.5* For the particle in Fig. 14 to be a 
proton, it is only necessary to assume that, in 
addition to the magnetic curvature of the track, 


5See following paper entitled “Cloud-chamber study 
of multiple scattering of protons in air.” 
*H. A. Bethe, Phys. Rev. 70, 821 (1946). 


LUEBKE, AND BALDWIN 


there is also a curvature due to multiple Scatter, 
ing, ‘this scattered curvature having a rading ‘ 
150 cm. Such a curvature due to scattering is ng 
unreasonable. A similar analysis of the fey 
tracks that gave mass values deviating a litt 
from the proton mass, showed in every Case thas 
the particle was most probably a proton, 

it was concluded that every measurable tray 
that was observed was caused by a proton, 

During the course of this work, severaj Dat 
ticles were observed that traversed the entip 
chamber, penetrating all of the foils. Such, 
track is shown in Fig. 15. While tracks such » 
this cannot be classed as measurable tragy 
because their range is not known, it seem 
reasonable, in view of the identity of all oth. 
tracks observed, to assume these particles to} 
protons also. If this is true, these protons hay 
energies in excess of 28 Mev.’ 

That such high energy protons are obserya 
is interesting in that it indicates that, in somed 
the photo-disintegrations induced by high enery 
quanta, it is possible for a single particle to cam 
off a large fraction of the excitation energy jj 
the nucleus, rather than this energy always being 
distributed among several particles. 

Figure 16 is a plot showing the distributioni 
energy of the protons observed in this wok 
Further work is now being carried out to examin 
the tail of this distribution and determine tk 
maximum energy of the photo-protons produce 
by 100-Mev x-rays. 


7 J. Smith, Phys. Rev, 71, 32 (1947) 
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Experimental evidence is presented for occasional large apparent curvatures of short range 
proton tracks in a cloud chamber caused by multiple scattering in the gas. The observed 
curvatures agree qualitatively with the predictions of Bethe. 





METHOD in common use for determining 

the mass of a particle in cloud-chamber 
studies involves the simultaneous measurement 
of the range of the particle and the curvature of 
its path in a magnetic field. It has been pointed 
out by Bethe"? that, in employing this method 
when dealing with particles having short ranges, 
multiple scattering of the particles by the atoms 
of the gas in the cloud chamber can cause large 
apparent curvatures of the tracks and thus in- 
troduce considerable error in the mass deter- 
minations. It is the purpose of this paper to 
present experimental evidence of the occurrence 
of these large apparent curvatures. 

The experimental arrangement used in this 
work is the same as described in the preceding 
paper.* The cloud chamber was filled to atmos- 
pheric pressure with air and water and alcohol 
vapors, and a series of pictures was taken without 
a magnetic field being applied to the chamber. 
A total of 1180 photographs was taken and on 
these pictures there were 129 heavily ionizing 
tracks similar in character to the photo-protons 
described in the preceding paper. In accordance 
with the results of the previous work, it seemed 
legitimate to assume these 129 tracks were all 
due to protons. 

Analysis of these 129 tracks gave the following 
results: 64 of the tracks were straight; 46 tracks 
showed visible curvatures approximating circular 
arcs, 21 of these curving in one direction and 25 
curving in the opposite direction; 19 of the tracks 
exhibited an S-shaped curvature. 

Of the 129 tracks, 18 were short ranged, stop- 
ping in the gas of the chamber before reaching 


"Presented in part at the September, 1946 American 
Physical Society Meeting in New York City. 

'H. A. Bethe, Phys. Rev. 69, 689 (1946). 

*H. A. Bethe, Phys. Rev. 70, 821 (1946). 
’ *G.S. Klaiber, E. A, Luebke, and G. C. Baldwin, Phys. 
Rev. (preceding article). 


the first foil. As predicted by Bethe,? these par- 
ticles which stopped in the gas of the chamber 
showed instances of large apparent curvature. 
As the ranges of these particles stopping in the 
gas is known quite accurately, the observed 
curvatures, 1/r, of the tracks can be compared 
with the root mean square curvatures, 1/r,. due 
to scattering, as calculated according to Bethe.? 
Making this comparison, one finds 9 tracks 
having curvatures less than the calculated cur- 
vature, i.e., 1/r<1/ree, and 6 tracks having 
curvatures greater than that calculated, i.e., 
1/r>1/ree. The remaining three tracks which 
terminated in the gas of the chamber exhibited 
S-curvature. 

Tracks stopping in the first .foil may have 
ranges between 15 cm and 44 cm, with corre- 
sponding values of r,. in the upper part of the 
chamber of 170 cm and 490 cm; particles which 
stop at the center of this foil, having a range of 
30 cm, should have an rz, of 350 cm. Of the 21 
tracks observed to end in this first foil, 6 showed 
S-curvature, 6 had r¢350 cm, and 9 had r>350 
cm. 

Several interesting cases of extreme curvature 
by multiple scattering are shown in the accom- 
panying photographs. Figure 1 shows a direct 
(upper) and stereo (lower) view of a track with 
range 13 cm of air and an observed radius of 
curvature of 60 cm. This curvature exceeds the 
calculated root mean square curvature by a 
factor 2.5. Figure 2 is an enlarged view of the 
same track. A straight line is drawn near the 
track for comparison. 

Another example of apparent curvature due to 
scattering is the track of shorter range presented 
in Fig. 3. It has a range of 10.5 cm air and an 
observed radius of curvature of 100 cm. This 
curvature is 1.2 times the calculated r.m.s. value. 
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Fic. 1. Track of a proton of 13-cm range and 60-cm radius of curvature taken in zero magnetic field. 
Fic. 2. Enlarged view of the track of Fig. 1 with straight line for comparison. 

Fic. 3. Track of a proton of 10.5-cm range and 100-cm radius of curvature taken in zero magnetic field. 
Fic. 4. Enlarged view of the shorter track of Fig. 3 with straight line for comparison. 





The enlarged view of this track together with a than 14 cm air and the observed radius of cur 


straight line for comparison appears in Fig. 4. vature is 80 cm. This exceeds the calculatei - 

The photographs of Figs. 5 and 6 present r.m.s. curvature by a factor greater than two. FA 
apparent curvature of the track of a particle that Figure 7 is an example of the “‘S’”’ curved track wn 
passes out of the illuminated portion of the cloud of a particle that stops in the first foil (rang - 


chamber. The range of this particle is.greater 12.5-42 cm air). Figure 8 shows this t rack with 


cu 








,- SCATTERING OF PROTONS IN AIR 
















d. Fic. 5. Track of a proton with range in excess of 14 cm, having an 80-cm radius of curvature in zero 
id magnetic field. 
. 1G. 6. Enlarged view of the track of Fig. 5 with straight line for comparison. 
Fic. 7. S-curved track of a proton of range between 12.5- and 32-cm air taken in zero magnetic field. 
The middle 5 cm of the visible track has a radius of 120 cm. 
Fic. 8. Enlarged view of the track of Fig. 7 with straight line for comparison. 


s of cur 

iIculated , : , —— 

a straight line drawn along it for comparison. An In view of this evidence for apparent curvature 
ed teak extremely short range track with large apparent of tracks due to scattering in the gas of the cloud 
il (range curvature is shown in Fig. 9. With range 3.5-cm chamber, it is of interest to consider the effect 
4c witht air and radius of curvature 12+3 cm, it has a_ on the range-momentum measurements reported 


curvature three times the calculated r.m.s. value. in the previous paper.* These effects are presented 
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Fic. 9. Heavily ionizing track of 3.5-cm range with 12+3- 
cm radius of curvature in zero magnetic field. 


in the curves of Fig. 10. The ordinate is radius 
of curvature and the abscissa represents the 
range in cm air equivalent to the extent of the 
five aluminum foils present in the cloud chamber. 
The extreme case is considered where the ob- 
served curvature (1/r.».) of a proton track is a 
combination of the magnetic curvature (1/rz) 
due to the 2650 gauss field used in the experiment 
plus three times the calculated curvature (1/fse) 
due to multiple scattering. In this case, the curves 
indicate that protons can be distinguished from 
mesons of mass 200 for ranges greater than about 
5-cm air equivalent in an experiment with mag- 
netic field 2650 gauss. 

The above examples of highly curved tracks 
are not sufficiently numerous for statistical 
analysis. However, they agree qualitatively with 
the predictions of Bethe and they indicate that 
the abnormally high curvature of several short 
range tracks which have appeared in recently 
published work** can be ascribed’ to the effects 


4 M. Schein, A. J. Hartzler, and G. S. Klaiber, Phys. Rev. 
70, 435 (1946). 
*D. J. Hughes, Phys. Rev. 69, 371 (1946). 
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Fic. 10. Radii of curvature of proton tracks of various 
ranges in the experiment of the preceding paper: rg, curve 
by the 2650 gauss magnetic field; r.., curved by multi 
scattering in air; fobs, curved by combination of magnetic 
field and extreme multiple scattering. These are compared 
with the radius of magnetic curvature of a meson track, 
shown in the lowest curve. — 


of multiple scattering in the gas of the cloud 
chamber. These curvatures, while not great 
enough to spoil the mass discrimination of the 
experiment with strong magnetic field described 
in the preceding paper, are sufficient to account 
satisfactorily for the few deviations from protonic 
mass encountered in that work. 

The authors wish to thank Professor H. A 
Bethe for his many helpful discussions on 
theoretical aspects, and Professor M. Schein for 
his interest in this work. 
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The gamma-ray yield from AF’ bombarded with 0.5 to 1.4 Mev protons has been investigated 


with a sensitive G.M. counter. By introduction of a “cutter” device which automatically 
stopped the measurements when the voltage deviated more than a certain amount from the 
desired value, a resolving power of 4 kev was obtained. Under this condition a number of 
previously known peaks have been resolved into groups of narrow lines. In the investigated 
region the proton width has been estimated to be less than 1 kev, and a rough estimate of the 


efficiency of the gamma-ray counter gave a radiation width of the order of 1 ev. 





INTRODUCTION 


HE reaction Al?’(p, y)Si®* was first investi- 
gated by Herb, Kerst, and McKibben." 
Using a thick target, they found a number of 
steps, which indicated the existence of resonance 
effects. Gentner? found two sharp resonances at 
425 and 560 kev. Plain, Herb, Hudson, and 
Warren? made a detailed investigation of the 
excitation curve up to 2.6 Mev, using a thin 
target. They found a multitude of resonances. 
Hole, Holtsmark, and Tangen** have found four 
weak resonances between 220 and 500 kev. 
Starting the work with the newly built Van de 
Graaff generator of the Institute for Theoretical 
Physics, Copenhagen, we intended to make a 
study of the excitation curve of this reaction, 
with the combined aim of matching the voltage 
scale with the results of the investigations 
mentioned above, and of finding new resonances 
with improved technique of detection. 


EXPERIMENTAL 
Generator 


The pressure-insulated Van de Graaff generator 
was built with a grant from the Carlsberg 


1R. G. Herb, D. W. Kerst, and J. L. McKibben, Phys. 
Rev. 51, 691 (1937). 

*W. Gentner, Zeits. f. Physik 107, 354 (1937). 

*G. P. Plain, R. G. Herb, C. M. Hudson, and R. E. 
Warren, Phys. Rev. 57, 187 (1940). ; 

*N. Hole, J. Holtsmark, and R. Tangen, Zeits. f. Physik 
118, 48 (194i), : 

*R. Tangen, Thesis, Trondheim, 1947 (to be published). 
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Foundation to Professor Bohr. It is of the 
Lauritsen® type and will be described in greater 
detail at another occasion. The plans were made 
in 1940, and the collaboration of Dr. T. Lauritsen 
who at that time worked at this institute was’ of 
invaluable help in the design and construction. 
The dimensions of the tank are 4.35X2.40 
meters, and the maximum pressure is 7 atmos- 
pheres. The two belts are 50 cm wide. They have 
independent motors and spray volt plants. There 
are 52 corona rings, and 26 sections in the accel- 
eration tube. The ion source is of the capillary 
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SPHERE 


Fic. 1. Rotating compensation voltmeter. 


* T. Lauritsen, C. C. Lauritsen, and W. A. Fowler, Phys. 
Rev. 59, 241 (1941). 
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Fic. 2. Target arrangement. 


type. The accelerated beam is analyzed mag- 
netically. 

A potential of approximately 1800 kv has been 
obtained with air at 6 atmospheres in the tank. 


Voltmeter 


The conventional rotating voltmeter measures 
the charging current of a disk, which is alter- 
nately exposed to and shielded from the electro- 
static field of the generator dome. The current is 
directly proportional to the voltage of the gener- 
ator, but its measurement depends on the re- 
liability of an amplifier or rectifier. 

This procedure has been changed into a com- 
pensation method using two plates P; and Pz» 
(Fig. 1). When the grounded wings W; and W, 
can be rotated without changing the potential of 
P,, which is grounded through a resistance, the 
generator voltage must be strictly proportional to 
the positive potential of P:, which is easily 
measured, being of the order of 1000 volts. 

If the compensation is not perfect, P; will have 
a small a.c. voltage. This voltage is amplified and 
indicated on a small neon lamp of the type used 
as tuning indicators in wireless sets. The amplified 
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signals from P; are also fed to the yo} 
stabilizer, which is of the McKibben type,” ang 
to the “‘cutter,’’ which blocks both gamma. 
counter and beam integrator when the Voltage of 
the generator deviates more than a prescribe 
amount from its proper value. 

The voltmeter is calibrated by means of the 
resonance in aluminum at 503 kev, whig, 
Tangen® has measured with high accuracy, both 
by comparison with the 440-kev resonance jy 
lithium, and by an independent calibration of his 
voltage scale. The expected linearity is ascertaingd 
by measuring the 503-kev resonance with both 
protons and the molecular ions H,* and H;+. 


Beam and target arrangement 


The ions were focused to a beam of 3 mm 
diameter and the current generally used was 1 to 
2 wA protons. The current integrator is of the 
type used by Tangen.® It consists of a recorder 
which counts the pulses in a neon lamp, dis 
charging a condenser constantly charged by the 
current. 

The target arrangement is seen in Fig. 2. S,and 
S: are stops of tungsten foil. A, B, and C are lead 
cylinders, which serve to shield the gamma-ray 
counter against x-rays from the generator. The 
bore of these cylinders is greater than the diame- 
ter of the stop S:, so that every ion coming 
through S, must strike the target. B is kept 
negative in relation to A and C, so that no 
secondary electrons from S: can pass to the target 
and vice versa. 

The lowest part of the target tube can be te 
moved for change of target. The target tube can 
be separated from the main vacuum system bya 
valve. The target was cooled by securing a good 
thermal contact between the target disk and the 
bottom of the lead box L. 

The lead box contained a single G.M. counter 
placed as near to the target as possible. A small 
tube is preferable in the detection of weak radia- 
tion,’ and a tube of 10-mm diameter and 40-mm 
length, was used. It was connected to a scale-of-32 


Preparation of targets 


Targets were prepared by evaporation d 
aluminum in vacum on disks of copper or silver, 


7D. B. Parkinson, R. G. Herb, E. J. Bernet, and J. L 
McKibben, Phys. Rev. 53, 642 (1938). 
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which had been given a fresh and clean surface 
by being turned carefully with a perfectly clean 
tool. 

As a multitude of targets of different thick- 
nesses had to be used, we have investigated the 
preparation of targets of known thicknesses in 
some detail. Because of the smallness of the 
resonance width, the resonances come out as ap- 
proximately rectangular curves when the stop- 
ping power of the target is much greater than the 
inhomogeneity in proton energy, and the breadth 
of this rectangle immediately gives the stopping 
power of the target. We found the target thick- 
ness to be proportional to the amount of alumi- 
num evaporated, so that targets of desired 
thicknesses could be produced with fairly good 
accuracy. 

In order to estimate cross section values we 
also determined directly the amount of aluminum 
in a 25-kev target, and found a specific stopping 
power of about 250 kev per mg/cm? at 0.5 Mev- 
proton energy, in accordance with Parkinson, 
Herb, Bellamy, and Hudson.® 

The aluminum had a specified purity of 99.998 
percent. The targets will of course have an oxide 
cover, but all measurements of peak heights have 
been made on 10-kev targets, which we suppose 
to consist for the most part of metallic aluminum. 
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MEASUREMENTS 


In determining the yield of the gamma-radia- 
tion over the energy interval examined, a great 
number of targets had to be used, since the 
targets were worn out after 2 to 3 hours of bom- 
bardment, even for a proton current of only 
1 to 2 pA. 

As a standard a target of about 1-kev stop- 
ping power at 500-kev proton energy was chosen. 
Targets of this type were used for all measure- 
ments on strong resonances (of intensities above 
0.5 in the curve), whereas regions of no, or weak, 
radiation were searched with targets of about 
4-kev stopping power. Finally a number of peaks 
throughout the spectrum were examined with 
targets thick enough (about 10 kev) to give 
saturation intensities. By these measurements 
the intensities obtained with the thinner targets, 
which could not be produced with exactly equal 
thicknesses, were brought to a common scale. 
The relative intensities of peaks are estimated to 
be reliable within 20 percent. 

In most measurements the cutter was adjusted 
to allow voltage fluctuations of +1.5 kv. With 
this voltage definition a 1-kev target gives about 
4 of the saturation intensity. 

The measurements on the 10-kev targets were 
also used for the establishment of the exact 














Fic. 3. Yield curve 
of aluminum bombarded 
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with protons. All peaks 
have been plotted with 
the intensities found with 
the standard targets of 
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about 1-kev stopping 
power (4 of saturation in- 
tensity), whereas the back- 
ground has been plotted 
with the intensities found 











with 4-kev targets. The 
natural counting rate has 
not been subtracted. 
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*D. B. Parkinson, R. G. Herb, J. C. Bellamy,”"and C. M. Hudson, Phys. Rev. 52, 75 (1937). 
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TABLE I. The first column gives the energy E in kev of 
the incident protons at resonance, while the second column 
contains the directly measured thick target yieid N in 
counts per microcoulomb. The last column gives the 
quantity wy from formula (1), where w generally is of the 
order 4, and where , except for the weaker resonances, 
deviates only slightly from the radiation width, when the 
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which are found in the curve of Herb et aj. 
resonance at 862 kev, which Herb ascribes to 
fluorine, has not been found, even though 
weak background radiation is present in that 
region. As no considerably stronger fluorine lines 


proton energy exceeds E~600 k 
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* Measured by Tangen*® and matched to our yield for the 503-kev 
resonance. 


resonance voltage. A simultaneous check of the 
voltmeter scale was made on the 503-kev reso- 
nance with H;*+, H;*, and H;* ions. The voltages 
given are supposed to be correct within 0.2 
percent on our voltage scale, which depends on 
Tangen’s value for the 503-kev resonance.® 

A survey of the results is given in Fig. 3 and in 
Table I, where the last column refers to the dis- 
cussion below. Only the range from 500 up to 
1400 kev is covered, since, owing to a break- 
down of the belts, the measurements above 1400 
kev are so far incomplete. 

The curve shows 31 prominent peaks, 16 of 


= exist in the examined voltage region, the curve 


can therefore not be influenced by fluorine cop. 
tamination. The complex peak at 1192 key seems 
to consist of three lines at 1189, 1192, and 1195 
kev, but nothing definitely can be said. 

The gradual rise of the background with jp. 
creasing voltage is mainly due to radiation from 
the silver disk. Copper gives considerable radia. 
tion above 1000 kev and has therefore not been 
used as target support at higher voltages. 


DISCUSSION 


According to the general scheme for nuclear 
reactions, the observed resonances correspond to 
stationary states of the compound nucleus §j 
formed and excited to about 12 Mev by the 
impact of a proton on a nucleus of Al*’, which is 
the only isotope present in natural aluminum. 
In the present case, where neutron escape is ex- 
cluded for energetic reasons, and where the 
emission of an a-particle is very improbable 
because of the large barrier, the final result of the 
collision will depend on a competition between 
the re-emission of a proton and radiation pro- 
cesses leading to a final capture of the proton. 

The yields of such resonances depend on 
the probabilities for these competing processes, 
specified by the so-called widths I’, and I, for 
proton emission and gamma-radiation, respect- 
ively, which are simply the probabilities meas- 
ured in energy units by multiplication by Planck's 
constant over 2x. According to the well-known 
formula of Breit and Wigner, the cross section for 
(p, vy) processes integrated over one of the narrow 
lines will be given by 

h? r,!; 
C=; v=, St 
4ME r,+l, 

where E is the resonance energy of the proton 
and M its mass, while w is a spin weight factor 
which, because of the high spin 5/2 of Al’, will 
be of the order }, although, for large angular 
momenta of the protons, it may deviate ap- 
preciably from this value. 

The thick target gamma-ray yield of a single 
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resonance, i.e., the fraction Y of the incident 
protons finally captured through the radiation 
processes, will therefore be given by the ex- 


pression 


= Swwy 
sMA 


=—, 2 
sAE @) 


where A is the mass number and s the specific 
stopping power of the target material, while so 
is an abbreviation for h?/4M? and equal to 2.5 
kev per mg/cm’. 

For the specific stopping power of aluminum 
we have used the values given by Herb e¢ al.,® 
which for 500 kev agree with our measurements. 

By comparing the yield from fluorine with the 
absolute yield given by Lauritsen e¢ al.,° and by 
van Allen and Smith,!° we find the efficiency of 
the counting arrangement to be 0.12 percent for 
6 Mev quanta. This corresponds approximately 
to 0.2 percent for the aluminum radiation. The 
estimate of the efficiency is rather rough and 
constitutes the main source of uncertainty in the 
values, given in Table I, of the quantity wy 
calculated by means of Eq. (2). 

From Eq. (1) it is seen that y will be of the 
same order of magnitude as the smallest of the 
two partial widths I’, and I,. In the low voltage 
region where, due to the influence of the barrier, 
I, is very small compared with I’, but increases 
rapidly with the proton energy, we should there- 
fore, according to Eq. (2), expect the resonances 
to rise gradually in strength with increasing 
energy E. 

In the high voltage region we may, of course, 
still find small values of I, corresponding to the 
impact of protons with high angular momenta, 
but the value of I’, for zero angular momentum 
should be large compared with I,. On entering 
this region the general increase in the resonance 
strength should therefore be expected to stop, 
since the radiation width is normally assumed to 
be approximately constant. 

Such a dependency on the proton energy is 
also indicated by the general trend of the yield 
curve (Fig. 3). The transition region between 
“low” and “high” energies appears to be situated 

* J. F. Streib, W. A. Fowler, and C. C. Lauritsen, Phys, 
Rev. 59, 253 (1941). 


J. A. van Allen and N. M. Smith, Jr., Phys. Rev. 59, 
501 (1941), 


in the vicinity of 600 kev. Below this energy the 
values given in the third column of Table I should 
thus represent wI',, while for higher energy they 
should represent wI',, except in cases where we 
have to do with comparatively large angular 
momenta of the incident protons. 

Since w will on the average be about 4}, one 
finds values of I’, of the order of magnitude of 
1 ev, in accordance with general theoretical esti- 
mates. For the three very strong levels at 986, 
1372 and 1379 kev it seems necessary, however, 
to assume excessively large values of T',, which 
might perhaps be caused by some specific charac- 
ter of the excitation states favorable for radiative 
transitions. 

As regards the order of magnitude of the proton 
width, we may assume that for zero angular 
momentum it will rise rapidly from about 1 ev 
for E~600 kev, to values of about 1 kev at 1400 
kev, as estimated from the penetrability of the 
barrier. Thus, for none of the resonances is the 
natural width, which is the sum of the partial 
widths, revealed by the measurements, but still 
the slope of the foot of the strong double-line at 
1372 and 1379 kev may give an indication of the 
actual width. In fact, by comparing the shape of 
the curve in this region with the foot of an ordi- 
nary dispersion line with the given area, one 
arrives at a total width of just about 1 kev. 

It need hardly be added that the “foot” might 
originate from one or more weak lines which have 
escaped detection by the present resolution, but 
there is reason to believe that the cutter arrange- 
ment is sufficiently reliable to prevent any ap- 
preciable influence from the central part-of the 
line at the distance in question, and this inter- 
pretation is also confirmed by the sharpness of 
the strong line at 986 kev for which, with the 
estimated much smaller proton width, no meas- 
urable foot should be expected. 

In the region from 600 to 1400 kev the mean 
distance between resonances is seen from Table | 
to be about 30 kev. The lines are within this 
region fairly evenly distributed, but there is a 
remarkable grouping and perhaps some regularity 
in the spacing of the levels of each group (e.g., the 
three lines at 757, 764 and 771 kev), although, it 
cannot be excluded that some of the coincidences 
are of a more accidental character. Still, especially 
in the case of the two very strong lines at 1372 











and 1379 kev, it is difficult to escape the impres- 
sion of the presence of a typical fine structure. 

In a forthcoming publication from this insti- 
tute it is planned to give a more detailed dis- 
cussion of the characteristics of the spectrum. 

The authors wish to express their thanks to 
Professor Niels Bohr for offering them the 
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for Danish-Norwegian Co-operation” and 

Norwegian “Foundation of 1919 for Scientific 
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in the work in Copenhagen. 

























been determined for several elements. 





INTRODUCTION 





HE scattering processes of slow neutrons 

are greatly complicated by interference 
phenomena due to the fact that the de Broglie 
wave-length is comparable with interatomic dis- 
tances. The general pattern of interference phe- 
nomena of slow neutrons is similar to that of 
x-rays, since both the wave-length and the scat- 
tering cross section of x-rays are comparable to 
those of neutrons. On the other hand there are 
considerable differences due to several factors. 
Among them is the fact that the scattering of 
X-rays varies regularly with atomic number while 
that of neutrons is a rather erratic property. 
Furthermore, in the case of neutrons the phase 
difference between scattered and incident wave 
may be either 0° or 180° as will be discussed in 
Section 1. For x-rays instead, it is always 180° 
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Various experiments involving interference of slow neutrons have been performed in order 
to determine the phase of the scattered neutron wave with respect to the primary neutron 
wave. Theoretically this phase change is very close to either 0° or 180°. The experiments 
show that with few exceptions the latter is the case. The evidence is based on the following 
types of measurements: (a) measurement of the intensities of Bragg reflection of various 
orders of many crystals, and comparison with the theoretical values of the form factor; (b) total 
scattering cross section of gas molecules for wave-lengths long compared with the molecular 
dimensions; and (c) determination of the limiting angles for total reflection of neutrons on 
various mirrors. The elements Ba, Be, C, Ca, Cu, F, Fe, Mg, N, Ni, O, Pb, S, and Zn were 
found to scatter neutrons with 180° phase difference; Li and probably Mn scatter with zero 
phase difference. The five elements I, Br, Cl, K, and Na behave alike and the phase with 
which they scatter is tentatively identified as 180°. Coherent scattering cross sections have 
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because x-ray energies are larger than most 
electronic resonance energies. Also, the absorp- 
tion properties of neutrons differ markedly from 
those of x-rays. 

The main purpose of this work was the investi- 
gation of various interference phenomena ia 
order to determine the phase change of the 
scattered neutron wave for a large number of 
elements. Section 1 contains a summary of the 
theoretical background of this work. Section? 
describes the measurements of the intensities of 
Bragg reflections of various orders and their 
interpretation. Section 3 is a discussion of some 
experiments on filtered neutrons. In Section 4, 
experiments on scattering of neutrons by gas 
molecules are presented. Section 5 describes 
measurements of the limiting angle for total 
reflection of neutrons. The general conclusions 
are discussed in Section 6. 
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1. Theoretical Considerations 


When a slow neutron is scattered by a nucleus, 
its de Broglie wave at some distance from the 
scattering nucleus may be written as the sum of 
a term exp(ikx) representing the primary wave 
and a term —al[exp(ikr)]/r representing the 
scattered wave. This is true for the case of slow 
neutrons because with very good approximation 
the scattering is spherically symmetrical. The 
coefficient of the latter term has been written 
with the minus sign for reasons of convenience 
that will be apparent later. If the constant a is 
positive, there is a phase change of 180° between 
scattered and incident wave, and if a is negative, 
the phase change is 0°. One proves in an ele- 
mentary way that the scattering cross section 
is related to the constant a by the equation 


o=4nrla|? (1) 


and also that the constant a is with very good 
approximation for slow neutrons a real number. 
Its imaginary part is very small and can usually 
be neglected except in case of extremely high 
absorption. The quantity a which has the dimen- 
sions of a length and the order of magnitude of 
10- cm shall be referred to as the “scattering 
length.” The main purpose of this paper is the 
experimental determination of the scattering 
length and in particular of its sign. 

For elements consisting of one isotope only 
and without nuclear spin, the magnitude of the 
scattering length can be immediately obtained 
from (1), so that only its sign needs to be deter- 
mined. Even in this simple case a small correction 
must be applied depending on whether the atom 
is free or bound. Because of the change in the 
reduced mass, the cross section of a free atom 
differs from that of the same atom bound in a 
crystal, by a factor [(A+1)/A }, where A is the 
atomic weight.! 

Since the scattering length is proportional to 
the square root of the cross section, the correc- 
tion factor will be (A+1)/A. It follows that 


a=[(A+1)/A Ja; (2) 


where a@ indicates the scattering length of the 
atom bound in a crystal and ay the scattering 


assy example see H. A. Bethe, Rev. Mod. Phys. 9, 71 
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length for the free atom. For elements that are 
mixtures of several isotopes with probabilities of 
occurrence pi, P2, «++, Pa, and scattering lengths 
1, G2, ++, @ the scattering length is the average 


A= P10, + podat +++ +Padn. (3) 


This is the magnitude that determines the inter- 
ference properties of the element. The scattering 
cross section is no longer given by (1) but by 


o =4r{ p,a;?+ pra? + - . *+Pnd,’}. (4) 


Actually in this case the relationship between 
a and a can be expressed by the inequality 


|a| <(o/4n)'. (5) 


The knowledge of the scattering cross section is 
insufficient, therefore, to determine even the 
magnitude of the scattering length. 

A similar situation obtains if the nucleus has 
a spin J different from zero. Here the spin vector 
of the scattering nucleus is oriented either 
parallel (total spin 7+ 4), or antiparallel (total 
spin J—}) to the spin vector of the neutron. 
Consequently the scattering length may take 
either of two values, @;4, or a;_;. The effective 
scattering length is the average of these two 
scattering lengths, each being weighted by the 
probability of occurrence of the corresponding 
spin orientation, namely 


I I+1 
62 —— 61-4 FT —0144- (6) 
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Fic. 1A. Positive scattering length. 
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Ata) 


a 


Pa 
Dil 


agi ones eine as 
Wiha a4 


Fic. 1B. Negative scattering length. 
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Fic. 2. Apparatus for measuring intensity of Bragg orders. 


The scattering length will obey (5), the equal 
sign corresponding to the case that a;,,; is equal 
to a;_, both in magnitude and sign. 

A simple geometrical interpretation of the 
scattering length can be obtained if the inter- 
action between the neutron and the nucleus is 
represented by a potential well. Let the eigen- 
function of the neutron for an s state of energy 
zero be ¥(r). In Fig. 1, rp(r) is plotted versus r. 
The curve becomes a straight line for r greater 
than the radius ro of the potential well. This 
straight line can be extended to intersect the r 
axis at the point P. One can prove that the 
scattering length a is the abscissa of the point P. 
Figure 1A corresponds to a case in which a is 
positive, and Fig. 1B to a case in which a is 
negative. By inspecting the two figures it is 
apparent that the probability is much higher for 
@ positive than for a negative, especially for 
heavy elements. This remark on the relative 
probability for a to be positive or negative is due 
to E. Teller and V. Weisskopf. 

A simple formula can be obtained in the case 
that the scattering is due to the effect of a single 
Breit-Wigner resonance level. It is, then 


a= —kXaIl',/[(W—R+1T ] (7) 


where W is the energy of the neutron, R is the 
resonance energy, I’, and I are the neutron and 
total half width at half maximum, and Xz is the 
de Broglie wave-length at energy R, divided 
by 2x. Since in most cases I’ is small compared 
with (W—R), it follows that a will be positive 
or negative depending on whether the energy of 
the neutron lies below or above the resonance 
energy. This picture is a great oversimplification 
since there are no cases in which the Breit- 
Wigner scattering of a single level is the dominant 
phenomenon, and complicated situations arise 
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due to interference of Breit-Wigner and potential 
scattering. 


2. Intensities of Various Orders of 
Bragg Reflection 


One can compare the changes in phase in the 
neutron scattering of different elements by meas. 
uring the intensity of the Bragg scattering fo, 
various orders and various crystalline planes of 
crystals containing at least two elements. 

The simplest case is that of a crystal in which 
the planes are equidistant and consist alternately 
of two kinds of atoms, as for example, the 1, 1, { 
planes of NaCl. These planes are equidistant and 
consist alternately of sodium and of chlorine. Ip 
the first order Bragg reflection, the optical path 
for reflection from sodium planes differs from 
the path for reflection from chlorine planes by 
4/2. Consequently if sodium and chlorine nuclej 
cause the same change in phase of the scattered 
neutron wave, their contributions will subtract 
and the order will have low intensity. If they 
scatter with opposite phase change their con- 
tributions will add and the order will have high 
intensity. 

The situation for the second order is reversed. 
Here the difference in optical path for reflection 
from the two kinds of planes is \. Consequently 
if sodium and chlorine scatter with the same 
change in phase the order will have high in- 
tensity, and if with opposite change in phase, 
low intensity. Similarly it follows that if the two 
kinds of nuclei scatter with the same change in 
phase the odd orders will be weak and the even 
orders will be strong. Conversely, if they scatter 
with opposite change in phase, the odd orders 
will be strong and the even orders weak. 

Analogously to x-ray scattering, there is super- 
imposed upon this effect a continuous decrease in 
intensity from order to order due to geometrical 
factors, to thermal agitation, and to imperfec- 
tions of the crystal. 

For more complicated cases the intensity of 
the various orders is determined as for x-rays 
by the form factor 


F= | 30; a; exp(2aind ;/d) | (8) 


where a; is the scattering length, d is the spacing 
of the lattice planes, is the order of the Bragg 
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reflection, and 4; is the perpendicular distance 
from the jth atom to the plane of reflection. 
Naturally the form factor is strongly dependent 
on the relative signs of the scattering lengths of 
the various atoms of the crystal. An analysis of 
the intensities for various orders and various 
planes of the crystal will often allow the determi- 
nation of the relative signs of the scattering 
lengths of the elements involved. In several cases 
it is possible to determine also the ratios of the 
values of the scattering lengths. 

A comparison of intensities for Bragg reflec- 
tions of various orders of many crystals have 
been made using monochromatic neutrons. The 
experimental arrangement is schematically repre- 
sented in Fig. 2. A beam of non-monochromatic 
thermal neutrons emerges from a 4” <4” hole 
containing a long collimator 4” wide by 1}” high 
in the thermal column of the Argonne heavy 
water pile. This beam falls on the 1, 0, 0 plane of 
a large CaF, crystal at a glancing angle of 
about 16°. Neutrons of energies satisfying the 
Bragg condition (A~1.50 angstroms for the first 
order) are reflected to form a beam that is used 
for investigating various crystal specimens. The 
beam contains a main component of neutrons of 
wave-length 1.50A plus a small fraction of neu- 
trons of energy four times as large due to the 
second order reflection on the CaF; crystal. The 
second order component is fairly unimportant 
because there is only a small number of neutrons 
of this relatively high energy in the Maxwell dis- 
tribution. We can, therefore, consider the beam 
reflected on the calcite to be approximately 
monochromatic. This monochromatic beam of 
neutrons falls on some plane of the second crystal, 
which is mounted on a rotating table. The 
neutrons undergo a second Bragg reflection, and 
are finally detected by an enriched BF; propor- 
tional counter. The counter is supported by an 
arm which rotates about the axis of the crystal 
table. Additional Cd slits not shown are used for 
more precise collimation of the neutron beam. 

A typical measurement was carried out as 
follows. Both counter and second crystal were 
set approximately at angles corresponding to a 
reflection of a given order. By rocking each 
separately the setting of maximum intensity 
was found. Also the inclination of the crystal was 
adjusted for maximum intensity. From the in- 
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tensity measured in these conditions was sub- 
tracted a background, measured by turning the 
crystal a small angle off the Bragg position. In 
the majority of cases the background was a 
small fraction of the total count. 

The results of these measurements are sum- 
marized in Table I. Columns 1 and 2 indicate 
the crystal and the plane. Column 3 is the form 
factor. Here the chemical symbol of the element 
has been used to represent its scattering length. 
The fourth column gives the absolute value of 
the form factor, using the values of the scattering 
length given in Table VII. The fifth column gives 
the measured intensities of the various orders in 
counts per minute. In the sixth column is the 
ratio of column 5 to column 4. If the values 
chosen for the scattering lengths are approxi- 
mately correct, one expects that the values in 
column 6, corresponding to any given plane, 
should show a regular decrease in intensity with 
increasing order, since dividing the intensity by 
the form factor should correct it for the irregular 
change of intensity from order to order. There- 
fore only the regular decrease should remain. 

There may be some doubt whether it is more 
appropriate to divide the intensity by the form 
factor or by its square. Theoretically one would 
expect that, for ideally perfect crystals, the form 
factor should be used, and for ideally imperfect 
crystals, its square. Actually, we have found that 
one obtains a much better fit by using the form 
factor. 

For the simplest crystals containing two ele- 
ments one can see immediately by inspecting the 
intensity data of column 5 whether the scattering 
lengths of the two elements have equal or 
opposite sign. For instance CaF, (100), NaCl 
(111), and PbS (111) show clearly superimposed 
on the general decrease of intensity with in- 
creasing order, an alternation of intensity with 
strong even orders and weak odd orders. Since, 
in these crystals, the planes consist alternately 
of the two kinds of atoms, we conclude that the 
scattering lengths of each pair have the same 
sign. The opposite is the case for LiF (111) where 
the even orders have intensity so low that we 
could not measure them, and the odd orders have 
a normal intensity. Therefore Li and F have 
scattering lengths of opposite sign. In some of 
the more complicated cases the analysis is also 
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TABLE I. Intensity of Bragg orders. 








Order Form factor Intensity 


Ca-—2F i 16300 
Ca+2 F { 20300 
Ca—2F 7 1287 


Na-Cl 2376 
Na+Cl 2750 


Pb-—S ‘ 7280 
Pb+S : 10700 
Pb—-S ‘ 808 
Pb+S : 750 


Pb+S > 19650 25800 
Pb+S ° 11420 15000 
Pb+S ‘ 2249 


Fe+ .62S ‘ 6893 
Fe—1.6 S ° 1354 
Fe—1.6 S o 726 


Mn—1.06S ; 7930 
Mn+ .06S ‘ 2560 
Mn+ .06S ‘ 670 
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12% 

1840 538 
10900 2730 
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quite straight forward as can be seen by com- 
parison of the observed intensities with the 
formulae for the form factors in column 3. In a 
few cases the interpretation is not unique. 

The following conclusions as to the signs of 
the scattering lengths can be drawn. The com- 
ponents of the following pairs have scattering 
lengths of the same sign: Na, Cl; Pb, S; Ca, F; 
Fe, S; Fe, O; K, Br; Mg, O; K, Cl; K, I. The 
components of the pair, Li, F have scattering 
lengths of opposite sign. The measurements of 
Mn&; do not allow a unique interpretation al- 
though the probable conclusion is that Mn and S 
also have scattering lengths of opposite sign. 

Using the evidence to be presented in Sections 
4 and 5 that carbon and oxygen have scattering 
lengths with the same sign, one can conclude 
from the measurements on calcite that Ca and O 
also have scattering lengths of the same sign. 

From the measurements on FeS; and Fes0,, 


one can conclude that S and O have the same 
sign; and from data on BaSQ,, one is led to 
assign also the same sign to Ba. The measure- 
ments on NaNO; indicate that N and O have 
the same sign. (See confirmatory evidence in 
Section 4.) The measurements of the (111) plane 
of NaNO; indicate although not quite con- 
clusively that Na has the same sign as the group 
NOs, and therefore as N and O. 

From these data it follows that C, O, Fe, Mg, 
Ba, Ca, S, F, Pb, and N all have the same sign. 
Na, K, Cl, Br, I all have like sign. The partial 
evidence just mentioned, according to which 
Na has the same sign as O, indicates that all 
this latter group should be included with the 
former. That this conclusion is correct is strongly 
supported by evidence to be presented in Sec- 
tion 6. Li and perhaps Mn have scattering 
lengths of the other sign. They represent together 
with hydrogen (see Section 6) the only exceptions 
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TaBLe I.—Continued. 

















Form factor Intensity 








BaSO. (100) 1 4Ba+4S+ 420 4.54 3606 794 
2 4 Ba+4S+14.6 O 13.19 4837 366 
3 4 Ba+4S+ 3.1 O 6.17 756 123 








BaSO, (011) 1 —3.4 Bat+3.8 S+4.70 1.24 1018 821 
2 1.75 Ba+3.1 S—960 3.61 1632 452 
3 42 Ba+2.0 S—7.90 3.93 651 166 
4 —2.46 Ba+ .75S—2.20 3.07 36 12 
101) 1 87 Ba—1.25S— 1470 _.56 737 1320 
BaSO. ( 2 1.13 Ba—2. 02 S— 4.560 2.46 1946 791 
3 3.77 Ba+ .98S+ 2.560 4.81 1477 307 
4 27 Ba— ‘0 S— 2120 1.47 on di 
5 1.04 Ba+2.12S—104 O 4.93 587 119 
100) 1 Mg+0O 1.17 10352 8800 
MgO ( 2 Mg+O 1.17 6258 5350 





MgO (111) 


1 
2 
3 





LiF (111) 


1 
2 
3 






KBr (111) 


ia 







KCl (111) 


1 
2 
3 
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NaNO; (211) 


1 
2 
3 





NaNO; (111) 


1 
2 
3 





KI (111) 


1 
2 
3 








Mg—O 05 764 15300 
Mg+0O 1.17 14175 12100 
Mg—O 05 132 2640 





Li-—F 1.19 10080 8470 
Li+F 1 ~~ _ 
Li-—F 1.19 300 252 








Br—K 21 1545 7360 
Br+K 91 2853 3140 


Ci-—K .78 2346 3010 
Ci+K 1.48 5620 3800 
CiI-—K 78 334 428 
Ci+K 1.48 160 108 


Na+N+ .880 1.96 29000 14800 
Na+N— .980 83 4171 5000 
Na+N+1.34 0 2.25 4148 1800 







NO;—Na 2.14 29400 13700 
NO;+Na 3.26 16790 5150 
NO;—Na 2.14 2243 1050 







K-—I 01 85 _ 
K+I 71 348 500 
K-I 01 24 _ 








found so far to the behavior of the majority of 
elements. 

This type of data allows only a comparison of 
the signs of scattering lengths of different ele- 
ments but not their absolute determination. The 
absolute determination of sign of scattering 
length will be discussed in Section 5. The over-all 
conclusions from Section 2 will be presented in 
Table VII. There, an attempt has been made to 
give the actual value of the scattering lengths 
for several elements. In calculating these values 
the attempt was made to obtain the best possible 
agreement of the observed intensities given in 
Table I, with the form factor after proper ac- 





count was taken of the natural decrease of 
intensity with increasing order. 

In x-ray analysis of crystals the Debye- 
Scherrer method of powder photography is ex- 
tensively used. For neutrons, however, the 
intensity of the Debye-Scherrer maxima is rather 
small, and in most cases the method will be 
impractical. For example, in the case of micro- 
crystalline graphite, we were able to detect only 
the maximum corresponding to the first re- 
flection on the 001 plane. 

A second possibility to study microcrystalline 
substances is to measure the total cross section 
of the substance for neutrons of wave-length 
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Fic. 3. Spectrum filtered through BeO. 


close to 2d where d is the longest lattice spacing. 
An analysis of this case for substances containing 
two elements has been made by Fermi and 
Sachs, and indicates that the total cross section 
in this region depends strongly on whether the 
scattering lengths of the two elements have equal 
or opposite sign. The practicability of this 
method is somewhat limited by the fact that 
the theory applies only when the microcrystalline 
grain of the substance is exceedingly fine. 


3. Spectrum of Filtered Neutrons 


In the study of interference in gas molecules, 
somewhat simpler results are obtained with 
neutrons of energy far below thermal energies 
and which have a wave-length long compared 
with interatomic distances. Such neutrons are 
present in very small percentage in the Maxwell 
distribution at room temperature. Consequently, 
one cannot isolate such low energy neutrons by 
Bragg reflection because in the reflected beam 
there is a contamination by high orders which 
is many times more intensive than the first order. 
A simple way to avoid this difficulty, without use 
of extremely low temperature moderators, is to 
remove the high energy part of the Maxwell 
distribution by passing the beam of neutrons 
through a filter of a microcrystalline material of 
low absorption,? such as graphite or BeO. This 


2 Anderson, Fermi, and Marshall, Phys. Rev. 70, 815 
(1946). 





E. FERMI AND L. MARSHALL 





filtering action of microcrystalline substances jg 
due to the fact that neutrons of wave-length 
longer than 2d, where d is the longest lattice 
spacing of the crystal, cannot satisfy the Bragg 
reflection condition for any of the many crystal- 
lites that they meet while crossing the filter, 
Consequently, they are scattered out of the beam 
only to a minimal extent. The production of very 
low energy neutrons by the filtering process has 
been previously reported? for the case of graphite, 
where the limiting wave-length is 6.7A. In the 
present experiments BeO filters have been used, 
for which the limiting wave-length is 4.4A. In 
order to get experimental proof of the correctness 
of the interpretation of the filtering process, we 
have investigated with the crystal spectrometer 
the spectrum of the neutrons filtered by a prism 
40 cm long and containing 100 g/cm?. The side 
dimensions of the prism were about 10X10 cm. 
The analysis was carried out with a crystal of 
celestite in the 001 plane. The results are plotted 
in Fig. 3, where the wave-lengths are plotted on 
the abscissae, and the observed intensities are 
the ordinates. The figure shows clearly that the 
spectrum of the filtered neutrons has negligible 
intensity for short wave-length, rises abruptly 
to a maximum value at wave-length of approki- 
mately 4.5A, and after that, decreases gradually 
with increasing wave-length. 

The decrease of intensity on the long wave- 
length side is somewhat more rapid than corre- 
sponds to the Maxwell distribution. This prob- 
ably is due to the fact that the neutrons emerging 
from the graphite column are not fully slowed 
down to the lowest energies of the Maxwell 
distribution by the collisions in graphite. 

The fact that microcrystalline BeO and Be 
metal are good filters for neutrons is an indication 
that the scattering lengths of Be for the two 
spin orientations of the neutron do not differ 
very considerably. Indeed if this were the case 
one would expect a strong incoherent scattering 
that should not vanish even for wave-lengths 
longer than 2d. For filtered neutrons we found a 
residual cross section for BeO of 0.7 X10- cm’. 
If all this residual cross section were due to the 
difference in the scattering lengths a; and a: for 
the two spin orientations, it could be calculated 
that a; and a; certainly have the same sign and 
that their ratio is about 2. Probably the actual 
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values are appreciably closer, because part of 
the observed residual cross section is certainly 
due to crystal imperfections and to thermal 


agitation. 
4. Scattering of Neutrons by Gas Molecules 


Some conclusions on the scattering of neutrons 
can be drawn from the study of the total cross 
section of gas molecules. For neutrons of wave- 
length short compared with interatomic distances 
the molecular scattering cross section is the sum 
of the individual scattering cross sections of the 
atoms in the molecule. When neutrons of longer 
wave-length are used, a number of complications 
arise, partly due to interference effects, and 
partly to the fact that the scattering now is due 
to atoms that no longer can be considered free. 

The theoretical calculation of the cross section 

has been carried out by Teller and Schwinger’ 
for the case of the He molecule for low rotational 
states. For heavier molecules, such as Nz or Os, 
even when low temperatures and slow neutrons 
are used, a fairly large number of rotational states 
is always excited so that a detailed calculation of 
the contributions of all rotational states becomes 
impractical. Simple results can be calculated by 
neglecting the neutron mass in comparison with 
the mass of the atoms in the molecule. In this 
approximation one obtains results identical with 
those of the classical theory of interference. The 
approximation is quite good in the case of x-rays, 
where the particle scattered is a photon of very 
small effective mass, and has been used currently 
by Debye and his co-workers. For scattering of 
neutrons from molecules such as Ne, Oz, etc., the 
approximation is not nearly as good. The cor- 
rections, however, have not been calculated, and 
the results will be compared with those of the 
approximate theory except in the case of H, 
where the theory of Teller and Schwinger has 
been used. 

The scattering cross sections of several mole- 
cules have been determined using very low 
energy neutrons obtained by filtration through a 
BeO filter. The spectrum of these neutrons is 
given in Fig. 3. Although these neutrons are not 
monochromatic they belong to a fairly narrow 
band of average wave-length 5.1A. For the 





* Teller and Schwinger, Phys. Rev. 52, 286 (1937). 


measurements at room temperature the gas 
under investigation was contained in a long 
aluminum tube of 2.5” id. and 365 cm long. 
Pressures up to 2 atmospheres were used. On H, 
we performed some measurements at liquid air 
temperature and in this case a tube 24” long was 
used. A copper coil was soldered around this tube 
and liquid air was circulated in it. The entire 
assembly was protected with rock wool insula- 
tion. The conditions were such that the ortho- 
para ratio was practically the same as at room 
temperature. 

Only in the cases of CF, and Hz measurements 
were performed also for neutrons of shorter 
wave-lengths obtained by Bragg reflection on a 
fluorite crystal. 

In Table II the total cross sections of BeO- 
filtered neutrons observed for various gases at 
room temperature are collected. In the third 
column of the table is given the sum of the total 
cross sections of the constituent atoms. One 
notices considerable differences between this sum 
and the observed cross sections. Except for the 
case of hydrogen which will be discussed later 
these differences are explained only to a small 
extent by the fact that the gas molecules are 
not at rest (so-called Doppler effect), and are 
due mostly to the interference of the waves 
scattered by the constituent atoms. In column 
4 the values of ¢ calculated from the classical 
interference theory (see above) are given. For di- 
atomic molecules containing two different atoms 
one obtains two different results, depending on 
the assumption that is made for the sign of the 
scattering lengths of the two atoms. In the cases 
of CO, and CF, a reasonably good agreement 
between theory and observation is obtained only 
if the scattering lengths of the component atoms 
are assumed to have the same sign. The same 
seems to be true for N,O, though the sensitivity 


TaBLE II. Cross sections for \=5.1A 
(gas at room temperature). 











Sum of the 
total cross 
sections of 
constituent ¢ calculated from classical! 
Molecule ¢ (observed) atoms interference theory 
CO: 24.5 13.0 24.8 or 4.1 for opposite phase 
N:O0 57.8 34 5S or 41 for opposite phase 
16.2 8.2 13.2 
or 13 1 3 7.5 for opposite phase 
‘ or or Oj te 
Hs 170 @ 4 

















674 E. FERMI AND L. MARSHALL 


TaBLE III. Cross section of CF, (room temperature). 











Velocity of 
neutrons 

»X108 (m/sec. ) o X10" 

881 4590 19.5 
1.099 3600 19.0 
1.316 . 3000 21.5 
1.492 2650 20.8 
5.1 775 41.5 








to a change in relative sign of the scattering 
length is not so pronounced here. 

The very large cross section observed for H: 
is strongly perturbed by the Doppler effect be- 
cause here the thermal agitation velocity of the 
molecules is larger than the velocity of the 
neutrons, so that the phenomenon can be better 
described by saying that the molecules hit the 
neutron rather than the opposite. This fact is 
brought out by the strong temperature depend- 
ence of the cross section which drops from 
170X10-* cm? to 8110-* cm? when the Hz 
is cooled from room temperature to liquid air 
temperature. This point will be discussed later. 

Table III gives the cross sections observed for 
CF, at room temperature for neutrons of various 
velocities. The table shows a rise by about a 
factor 2 when the wave-length increases from 
1.5 to 5.1A. Up to 1.5A the cross section is fairly 
close to 21X10-* cm?, namely to the sum of 
the cross sections of one atom of carbon and four 
of fluorine. The rise is due in large measure to 
interference. 

Tables IV and V give the results obtained on 
H: at room temperature and at liquid air tem- 
perature respectively. At room temperature one 
notices a steady rise in cross section with de- 
creasing velocity of the neutrons. ‘As previously 
indicated, this is to a large extent due to the 
thermal agitation velocity of the gas molecules. 


TABLE IV. Cross section of H; at room temperature. 





A correction factor, calculated on the assumption 
that the cross section is constant, is given jp 
the fourth column. The last column gives the 
cross section corrected by this factor. 

The data of Table V taken at liquid air tem. 
perature can be most easily compared with the 
Teller-Schwinger theory because at this tem. 
perature practically only rotational states 0 and j 
are involved. Column 3 of this table gives the 
observed cross section and column 4 gives the 
cross section corrected for the thermal agitation 
velocity of the gas molecules. The fifth column 
gives the theoretical values calculated with the 
Teller-Schwinger theory, assuming d)=—2,40 
<10-" and a,=0.54X10-" for the singlet and 
triplet state of the neutron proton system. The 
agreement is not too good especially for the 
lowest velocity, though the last value may be 
vitiated by the large Doppler correction. 

By performing a similar set of measurements 
on deuterium one might be able to draw some 
conclusions of the spin dependence of the cross 
section, and we hope to make such measurements 
in the future. 


5. Total Reflection of Neutrons of Mirrors 


The total reflection of neutrons on mirrors is 
theoretically expected to occur at very small 
glancing angles for substances with a positive 
scattering length. The index of refraction is 
given by 


n=1—)*Na/2r (9) 


where a is the scattering length taken with the 
proper sign, N is the density of atoms in the 
mirror, and d is the wave-length. Consequently 
is less than 1 (case of total reflection) if a is 
positive. The limiting glancing angle is then 


- TaBie V. Cross section of Hz: at 83°K (ortho-para ratio 


unchanged from room temperature). 











Factor for 
Velocity Maxwell 
4 X<108 (m/sec.) o X10% distribution o/factor 
-76 5212 50.9 1.048 48.6 
.93 4256 53.2 1.071 49.7 
1.11 3546 56.2 1.099 51.1 
46 2714 61.2 1.172 52.2 


1 

1.68 2347 64.4 1.231 52.3 
2.24 1765 68.0 1.407 48.3 
5.1 775 169.7 2.514 67.5 





absorption Teller- 
Velocity total X10*% and Doppler Schwinger 





»X108 (m/sec.) observed effect theory 
0.836 4730 44 43 _ 
1.202 3290 . a 49 _ 
1.638 2420 55 52 56 
1.884 2100 56 51 57 
5.1 775 81 52 80 
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‘Fermi and Zinn, Phys. Rev. 70, 103 (1946). 
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Fic. 4. Monochromatic total reflection on mirrors. 


given with very good approximation by 
69=[2(1—m) ]#=A(Na/z)!. (10) 


In most cases (w —1) is of the order of magnitude 
of 10-* and the limiting angle of the order of 10’. 
The total reflection of neutrons on mirrors has 
been previously observed‘ using non-monochro- 
matized thermal neutrons. In this case of course 
one cannot observe a limiting angle because 4 
[see (10) ] is proportional to \ so that neutrons of 
the various wave-lengths of the Maxwell dis- 








tribution drop out of the reflected beam gradually 
as the glancing angle is increased. In spite of this, 
the fact that strong reflection is observed on a 
substance can be construed as good evidence 
that the index of refraction for that substance is 
less than one, and that therefore a is positive. 
This is so because, although a substance which is 
not totally reflecting has a finite reflection coeffi- 
cient, it is so small that no very prominent re- 
flection could be observed in such cases. From 
the quoted experiments it could therefore be 
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Fic. 5. Typical results of measurements of a reflection on a Be mirror. 
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TABLE VI. Limiting angle for total reflection of 
neutrons of 1.873A. 








Limiting angle (minutes) 
bserved Caloula 





Mirror 
Be 12.0 11.1 
C (graphite) 10.5 8.4 
Fe 10.7 10.0 
Ni 11.5 11.8 
Zn 7.1 6.9 
Cu 9.5 9.5 








concluded that the following elements, Be, Cu, 
Zn, Ni, Fe, C, for which strong reflection has 
been observed, have positive scattering length. 

In order to test conclusively the theory of total 
reflection, the reflection experiments were re- 
peated using monochromatic neutrons. This 
made it possible to observe the sharp drop of 
the reflected intensity at the limiting angle and 
to measure this angle. 

The experimental arrangement is shown sche- 
matically in Fig. 4. A collimated beam of thermal 
neutrons emerging from the thermal column is 
made monochromatic by Bragg reflection on a 
fluorite crystal. It is further collimated by 
passing through two very narrow cadmium slits, 
S; and S;, 2 mm wide and about 3.5 meters 
apart. The narrow monochromatic beam emerg- 
ing from the second slit falls on the mirror which 
can be rotated by small angles. The reflected 
beam is detected by a BF; proportional counter 
about 3.5 meters away and with a cadmium 
entrance slit 2mm wide. By moving the counter 
across the beam, both direct and reflected beam 
are detected. 

Typical results of these measurements are 
plotted in Fig. 5 for the case of a reflection on a 
Be mirror. The abscissa is the position of the 
detector expressed in centimeters measured trans- 
versely to the direction of the beam. The ordinate 
is the neutron intensity in counts per minute. 
The curves whose experimental points are marked 
0, 2, 4, etc., correspond to the intensity plots 
found when the mirror is set at an angle of 0, 2, 4, 
etc., minutes to the direction of the neutron 
beam. When the mirror is parallel to the beam, 
there is observed a single maximum correspond- 
ing to that of the main beam. When the mirror 
is rotated by 2’ the reflected beam is barely 
resolved from the direct beam. At 4’ the two 
maxima are well separated. For larger angles the 
direct beam has not been drawn in the plot 
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TaBLe VII. Scattering lengths. 











tgs 
Ba 0.79 7.7 8 coll 
Be .89 8.0 6.1 M 
Br 56 3.9 <7 Cc 
. 67 48 48 A 
a “3CtaS Frog 
K 1 A, no spin and isotope 
effect considered 
= - 2” 7.2 M ” 
d " 41 A, no spin effect 
sidered sa 
Fe 82 8.1 9.2 M 
H —.39 48 21 From 
I 36 1.6 1.6 A, no spin effect con. 
K 35 1.5 1.5 A, no spin and isotope 
effect considered 
Li —.59 3.4 i Cc 
# 2 8 & fneee 
n =. 4 J » no spin 
N 87 8.3 8.3 A | effect 3 
d \ J »no 
sidered = 
Na 56 3.7 3.5 Cc 
Ni 1,09 14 13 M 
Oo 61 4.1 4.1 A 
Pb 48 2.9 10 © 
Ss .28 95S 1.1 Cc 
Zn 58 4.1 3.6 M 





———— 
—<—=} 





except for the largest angle, 13’, at which the 
reflected beam has completely disappeared and 
also the direct beam has disappeared because the 
mirror has turned by such a large angle that it 
completely cuts off the main beam. These points 
therefore correspond to the background, which 
actually was about 20 counts per minute. 

By inspecting the plot one sees that the in- 
tensity of the reflected beam drops sharply to 
background between the positions 10 and 13 
minutes. This drop is not quite sudden because 
the angular resolution was about two minutes. 
A study of the intensity shows the limiting angle 
to be 12 minutes. 

In Table VI the values of the limiting angles 
of several mirrors measured for neutrons of de 
Broglie wave-length 4\=1.873A are given. The 
experimental values given in the second column 
of the table should be compared with the theo- 
retical values calculated from formula (10) given 
in the third column. The agreement is excellent 
except in the case of carbon, for which a not too 
good graphite mirror was used. 

The absolute determination of the sign of the 
scattering length for the totally reflecting ele- 
ments, taken together with the results of Sections 
2 and 4, allows one to determine the actual sign 
for a large group of elements. 


6. Conclusions 


In Table VII all the data available at present 
on scattering length have been summarized. The 
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absolute sign given in the table has been ob- 

tained primarily from the data of Section 5 on 

total reflection, and indirectly, by combining 

these data with those of Section 2 on intensities 
of Bragg reflections, and Section 4, on molecular 
scattering. The absolute values have been ob- 
tained in various ways. 

For the elements having only one isotope or 
one strongly dominant isotope, and no nuclear 
spin, the absolute value of the scattering length 
can be calculated from formulae (5) and (2) 
provided the total scattering cross section is 
known. For example, this method has been 
applied in assigning the value 0.61 X 10-* cm for 
the scattering length of oxygen in a crystal. This 
differs by a factor 16/17 from the scattering 
length 0.57X10-" cm for a free oxygen atom. 
In some cases this same method has been applied 
to elements having several important isotopes, 
or a spin, or both. In such cases the absolute 
value is naturally much less reliable. The method 
has been used when no better way was available. 
In the last column of Table 7 the code letter A 
has been used to denote values obtained by this 
method. 

The results of Section 2 allow the scattering 
lengths of several elements to be related directly 
or indirectly to that of oxygen whose scattering 
length 0.61 <10-" has been used as a standard 
for these cases. 

For the five elements Na, K, Cl, Br and I, 
which are not conclusively related to the others 
(see Section 2), Cl has been taken as a standard 
with a scattering length of 1.1310-" cm. The 
sign has been taken positive partly on account 
of the evidence from Bragg reflection on NaNOs, 
as discussed in Section 2. Furthermore, both 
theoretical and experimental evidence indicate 
that cases of negative scattering length are very 
improbable. The probability that five elements 
all having the same sign should be negative is 
therefore exceedingly small. 

In other cases the value of the scattering length 
has been calculated from the value of the limiting 
angle using formula (10). Code letter M indicates 
these instances. 

In the third column of Table VII is given the 
expression 4ra*(A/A+1)*X10™. If there are no 
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important isotopes and no spin this should 
coincide with the scattering cross section. Other- 
wise it should in general be smaller, the moreso 
the greater the differences are in the scattering 
lengths for the various isotopes and spin ori- 
entations. 

The values in column 3 therefore should be 
compared with the experimental values of the 
scattering cross section given in column 4. 
Naturally the comparison is trivial for the cases 
designated by A, where formulae (2) and (5) 
have been used to calculate a. Comparison in 
the remaining cases shows that usually the data 
of the third and fourth columns are rather close. 
The few exceptions for which the data in column 
3 are less than those of column 4 can be explained 
as due to experimental inaccuracy. The fact that 
they are close seems to indicate that, in general, 
the spin and isotope dependence of the scattering 
length is not extremely pronounced. 

There are two notable exceptions. One is 
hydrogen, a case well known from the Teller- 
Schwinger theories. The other one is lead where 
perhaps some of the many isotopes may have a 
negative scattering length. 

The table shows that the scattering length is 
positive in the great majority of cases as dis- 
cussed in Section 1. 

This work was performed at the Argonne 
National Laboratory, and we greatly benefited 
by the help of the staff of the laboratory in 
providing and operating the neutron source, and 
in supplying many of the experimental facilities 
needed. Several of the crystals used were loaned 
to us by the Field Museum in Chicago through 
the courtesy of Mr. Changnon, by Dr. Howland 
of the Minerology Department of Northwestern 
University, and by Dr. Fisher of the Geology 
Department of the University of Chicago. Also 
some specimens were obtained from the Union 
Carbide and Carbon Research Laboratories, Inc., 
and the Eagle Picher Lead Company. Mr. 
Warren Nyer helped in some of the experiments. 

This document is based on work performed 
under Contract No. W-31-109-eng-38 for the 
Manhattan Project, and the information covered 
therein will appear in Division IV of the Man- 
hattan Project Technical Series as part of the 
contribution of the Argonne National Laboratory. 
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Neutron Scattering in Ortho- and Parahydrogen 


Morton HAMERMESH 
New York University, University Heights, New York, New York 


AND 


JULIAN SCHWINGER 
Harvard University, Cambridge, Massachusetts 


The cross section for the scattering of cold neutrons by parahydrogen is sensitive to the 
range of the forces in the triplet state. To facilitate the determination of the range from such 
measurements, expressions have been obtained for the dependence of the ortho and para cross 
sections on neutron energy and hydrogen temperature. The interpretation of transmission 
measurements is complicated by the existence of radiative capture, which is comparable in 
importance to scattering in parahydrogen. The different energy dependence of the two 
processes should permit an accurate determination of the individual cross sections. 





oo papers' have emphasized the importance of measurements of neutron scattering in 
gaseous ortho- and parahydrogen, as a means for obtaining accurate information concerning 
the range of nuclear forces. Preliminary experiments were performed by Alvarez and Pitzer,? but with 
insufficient precision for this purpose.* In anticipation of more accurate experiments, a letter to 
The Physical Review* presents curves showing the dependence of the ortho- and parahydrogen cross 
sections on the range of the nuclear forces and the free-proton cross section. These curves must be 
corrected for three reasons: 

1. The calculation of the triplet amplitude, a1, was based on the formula: a;= —(1+ar»)!/a, 
where a=[ME,)/h*]!, Eo is the binding energy of the deuteron, and 1 is the range of a rectangular 
potential well interaction. This expression is an approximate one, whose validity is restricted to 
ranges that are small in comparison with 1/a=4.35X10-" cm. This approximation may lead to 
errors of as much as 1 X 10~"* cm in the range deduced from the observed para scattering cross section. 

2. The curves are based on equations given by Schwinger‘ for the cross sections at a hydrogen gas 
temperature T= 20°K, and neutron energy E=kT. The constants in these equations were obtained 
by an approximate procedure and are in error by a few percent. ‘ 

3. The underlying theory has been based on a rigid-rotator model for the Hz molecule. This 
approximation introduces small but not negligible errors.** 

Accurate calculations of ¢para and gortho have been made using essentially the same procedure as for 
Dz,* but modified to include the effect of zero-point oscillations. The only change in the theory is the 
replacement of jz(kr.) by (jx(kr))w, the mean value of jz(kr) in the lowest vibrational state. In the 
calculation of this average we have performed an expansion in powers of y=r—r,/r, and dropped 
terms beyond g and (y’),. It should be noted that 7 must be included since the zero-point vibrations 
are anharmonic. (y*), is obtained directly from the harmonic oscillator model which gives 
(y")\w=4E:/hw where hw(=0.513 ev) is the excitation energy of the first vibrational state, and 
E;(=0.0147 ev) is the energy of the J=1 rotational level of Hz. A numerical value of 9 can be ob- 
tained by considering the dependence of the rotational energy levels on the vibration quantum 
number. One need only compare r,=0.7414 X10-* cm, the equilibrium separation, with ro =0.7506 


1 J. Schwinger and E. Teller, Phys. Rev. 52, 286 (1937); J. Schwinger, Phys. Rev. 58, 1004 (1940). 

2L. W. Alvarez and K. S. Pitzer, Phys. Rev. 58, 1003 (1940). ; 

* Recent experiments at Los Alamos have shown that ortho-para conversion during the course of the experiments is @ 
serious source of error. The data of these investigators, which differs considerably from that of Alvarez and Pitzer, will 
be published shortly. 

*C. S. Wu, L. J. Rainwater, W. W. Havens, Jr., and J. R. Dunning, Phys. Rev. 69, 236 (1946). 

‘); Schwinger, Phys. Rev. 58, 1004 (1940). 

** We are indebted to E. Teller for discussions of this point. 

5 M. Hamermesh and J. Schwinger, Phys. Rev. 69, 145 (1946). 
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Fic. 1. Dependence of parahydrogen cross section on the Fic. 3. Parahydrogen cross section asa function of neutron 
range of nuclear forces. (Gas temperature T=20°K, neu- energy, omitting amplitude factors (cf. Eq. (1)). (Gas tem- 
tron energy E=kT.) perature = 20°K.) 


10-8 cm, which is defined by 1/ro? = (1/r?)w1/r.2(1 —29 +3(y")a). This procedure gives J = 0.0337, 
(y?) = 0.0143. 

The cross sections for the various transitions have been expanded in powers of E/E,, where E is 
the neutron energy. The thermal averages of the cross sections for the various transitions are: 
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Fic. 2. Dependence of ortho- 
hydrogen cross section on the 
range of nuclear forces. (Gas 
temperature T=20°K, neutron 
energy E=kT.) 
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1 exp(—x*) 


Cee ——— ee ee @) 1429+ 0%) =| A(14- = om ta 


x 


+(45+5)|+ aaa 14494-10142) (=) 


1 7 33\exp(—x*) 15 
x{5(14+5+5 +1 ttt )e +) 
i at x? 4x4 x 2x? 4x* 8x6 
e111 Barta)? +2 (01 — a0)" 1Gin 


where G;--; is identical with Goo, except that the factor preceding (Z/E;)? is replaced by 
1024 


(1449+ 19140): 
(81) (225) 4 es 


644 E,\3 
ow1=—(1~09)*(—) [o.01982(1+1.739-0.15(9%)) 


Phas 1S /kT 
+0.02719(1 + 1.239 —0.664(y? r= (1 +) +0. oooass(—), {1+- "eae —) | (1) 
1 


@;,9 are the triplet and singlet amplitudes, T is the temperature of the gas in °K, x*=2E/kT and 
2 z 
(x) -—{ exp(—?*)dz. 
at 0 


For T=20°K and E=kT we obtain: 
Oo-0 = 6.444(3a,;+a)?, 73-1 =6.450[ (301+)? +2(ai1—ao)*], oo-1 = 1.753(a1—a)?, 


so that 
Tpara = 6.444(3a1+40)*, Fortho = 6.450(3a1-+-a0)? + 14.653 (a1 —ao)?. 


The dependence of a; on ro was computed for a rectangular potential well, using the expansion 
given by Kittel and Breit.* Then a» was calculated for several values of the free-proton cross section, 
ty. Lhe results are shown in Figs. 1 and 2. 

The parahydrogen scattering cross section is obtained from transmission measurements which 
determine only the total cross section of the molecule including the radiative capture cross section. 
From the measured value of the mean lifetime of slow neutrons in water,’ the capture cross section 
of the hydrogen molecule for 20° neutrons is calculated to be 2.5 X 10-* cm?, which is approximately 
equal to the parahydrogen scattering cross section corresponding to the range of nuclear forces 
customarily assumed in current theory, (ro=2.8X10-" cm). Transmission experiments at several 
neutron energies will enable a separate determination of the scattering and capture cross sections 
since they differ in their energy dependence. The capture cross section follows a 1/v law, whereas 
the scattering cross section varies with energy as shown in Fig. 3. Thus, for example, the capture 
cross section for 80° neutrons is one-half of that for 20° neutrons, whereas the scattering cross section 
at the higher neutron temperature is only 25 percent less than that at 20°. 


*c. bg and G. Breit, Phys. Rev. 56, 744 (1939). 
7J. H. Manley, L. J. Haworth, and E. A. Luebke, Phys. Rev. 61, 152 (1942). 
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The Double Focusing Beta-Ray Spectrometer 


FRANKLIN B. SHULL AND Davip M. DENNISON 
Randall Laboratory of Physics, University of Michigan, Ann Arbor, Michigan 
(Received February 10, 1947) 


The double focusing spectrometer recently proposed by 
Siegbahn and Svartholm has been analyzed and compared 
with the conventional semicircular spectrometer. A shaped 
magnetic field is required which possesses cylindrical 
symmetry as well as symmetry with respect to the median 
plane of the instrument. The magnetic field in the median 
plane is assumed to have the form 


H=H)—(r—a)aHo/a+(r—a)*BHo/a?’, 


where a and @ are experimentally disposable shaping 
constants. The analysis, which has been carried through 
the second order of approximation, shows that double 


focusing occurs when a=}. In general there will exist, for 
all types of spectrometers, second order defocusing terms 
which tend to lower the resolving power. However, these 
may be almost completely eliminated in the double 
focusing spectrometer by choosing B=}. A figure is 
presented showing the form of the focused image for a 
semicircular spectrometer and for several cases of the 
double focusing spectrometer. It becomes evident that 
when a=4$ and 8=}, the intensity of the image as well as 
the resolving power become many times greater than for 
the former type of instrument. 





1. INTRODUCTION 


T is customary in beta-ray spectroscopy to 

focus electrons with either the magnetic semi- 
circular spectrometer or the magnetic electron 
lens. The former requires a uniform field, which 
is created by an iron-core magnetic circuit with 
plane parallel pole faces. Large magnetic field 
strengths are feasible, and hence high energy 
electrons can easily be focused. The electron lens 
uses a magnetic circuit which is relatively free 
of iron. Much copper and power are necessary 
to produce large fields, and it is impractical to 
work with high energy electrons. In another 
respect, however, the electron lens seems su- 
perior. Assuming equal resolving power for both, 
the electron lens can focus a beam intensity as 
much as 15 to 20 times greater than the semi- 
circular spectrometer. 

Siegbahn and Svartholm! have proposed and 
constructed an electron focusing system which 
utilizes an inhomogeneous magnetic field pro- 
duced by an iron-core magnetic circuit with 
shaped pole faces. The magnetic field is shaped 
to fall off radially as 1/r*. The electron beam is 
focused after sweeping out an angle 2'r about 
the axis of symmetry. The focusing principle is 
basically a refinement of the semicircular system 
which permits a considerable increase of beam 
intensity without sacrificing resolving power. 
Such a spectrometer could compete more favor- 





'Siegbahn and Svartholm, Nature 157, 872 (1946). 
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ably with the electron lens in terms of focused 
beam intensity, while retaining the ability to 
focus high energy particles. 

It is proposed in the paper which follows to 
analyze this focusing system, evaluate its per- 
formance, and compare it to the conventional 
semicircular type. The equations of motion of 
the electrons will be solved by means of the usual 
perturbation methods, and the calculations will 
be carried through the second order of approxi- 
mation. 


2. THEORY 


It will be convenient to use the cylindrical 
coordinates r, 0, and z, and it will be assumed 
that the magnetic field is (a) independent of the 
azimuth angle @ and (b) symmetrical with respect 
to the median plane z=0. The vertical compo- 
nent of the field, H,, as measured in the median 
plane, is taken to have the form 


r—a r—a\? 
H,=Hy~(—)aie+(—) BHo—---, (1) 
a a 


where a is the radial coordinate of the center of 
the electron source (whose other coordinates are 
6=0, z=0), Ho is a constant magnetic field, and 
a and @ are arbitrary dimensionless constants. 
Clearly Eq. (1) represents a series development 
for H, through second order, and assumes that 
the equilibrium or zeroth order orbit will be a 
circle with the radius r=a. If, for points in the 





682 F. B. SHULL AND 
median plane, H, is assumed to fall off as 1/r*, 
the identifications a= and B= (n?+n)/2 can be 
made, but it should be emphasized that Ho, a, 
and 8 are experimentally disposable parameters 
in the present problem. Since curlH =0, it follows 
that the radial component H, of the magnetic 


field must be 


Hea stt{ -*+ a 2p(r— oe 


and since divH =0, it follows that a fourth term 
must be added to Eq. (1) in order to express H, 
at points off the median plane. Thus, 


r—a r—a\? 2? 
H,=Ho— (—)ai+(—) BH »——BHo. 
a a a? 


In any purely magnetic field, the motion of a 
charged particle must be such that its velocity, 
and consequently its mass, will remain constant. 
The equations of motion for the particle then 
reduce to the following expressions: 

di e 
——rb?=—+6H,, 
dt mc 


(2) 


“(7#)= Bs GH, —?fH,), 
mc 


(3) 


dz 


e 
= ——76H,. 
dt mc 


(4) 


As is well known, an integral of the motion may 
be obtained by multiplying (2) by 7, (3) by 86, 
and (4) by z, summing, and integrating with 
respect to the time. This integral is 


f+r¢+ 2% =y* =constant (5) 


where v is the initial velocity of the electron. 
Any three of the Eqs. (2), (3), (4), and (5) would 
be sufficient to describe the motion, and it will 
prove convenient in the present case to eliminate 
(3) and use only (2), (4), and (5). 

Assume that an electron originates in the 
source at the point r=a, z=0, @=0 and has the 
initial velocity components *=0, z=0, rd=v 
= —¢H,a/mc. Under these conditions, the solu- 


tion of the equations of motion becomes 
r=a, 2=0, 0=wt, 


where w=v/a. This circular orbit is hereafter 
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designated as the equilibrium or zeroth order 
orbit. In general, however, the initial conditions 
describing an electron orbit will be 


r=a+ir 
z= dz 
6=0 
t/v =¢r 
z/v = @¢: J 
where ¢, and g, are evidently angles which 
specify the direction of the initial velocity, ang 
ér and éz are coordinates of the initial position 
relative to the center of the source. 

It is proposed to solve Eqs. (2), (4), and (5) by 
standard perturbation methods, and to obtain 
solutions for r, 8, and z which are complete 
through second order. Let 


r=a+)Apityr’ps, 
6 = w+: +0%0s, 6) 
z= NMitr*Se, 


where pi, 6:, {1 and pe, 62, {2 are functions of the 
time, and \ is a parameter of smallness. It will 
be noticed that the equations of motion involye 
6 rather than 6, and consequently 6 has been 
specified above in place of @. The variables (6) 
are substituted directly into Eqs. (2), (4), and 
(5), and the three resulting equations are classi- 
fied according to powers of the parameter X. By 
equating the coefficients of \°, \', and 2? to zero, 
equations are obtained which describe the equi- 
librium orbit and the first- .and second-order 
corrections to the equilibrium orbit respectively. 

The zeroth-order equations yield the following 
relations, which determine the magnetic field 7, 
as a function of the momentum of the electron. 


w= —eH/me. (7) 


The first-order equations are 





v=wd, 


dpi 
A—+Aw?(1 — a) pi =0, 
dt 


df, 
A— + Aw*at =0, 
dt 


6:= —wp;/a. 


It is evident that the first-order solutions Apy 
6, and Af; will be sinusoidal, with circular 
frequencies w(1—a)*, w(1—a)!, and wat, respec: 
tively. 
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The second-order equations become after some simplification 


dp: w? Bow* 1 , 
Trae —a)ps] = | 08 —(ort- $1?) ——(o:?2+ | (9) 
dt a a 2a 


df Bw? 
| —+w%ats| = 2\*— pif, (10) 
dt a 
w w (617+ £1") 
762 =A*] —p1? ——p2 -—_——— |. (11) 
a? a 2a*w 


Both Eqs. (9) and (10) consist of a homogeneous part, to the left of the equality sign, and an inhomo- 
geneous part. The solution of the homogeneous part of each may be neglected, since it is identical 
to the solution of the corresponding first order Eq. (8), and thus contributes nothing new. The 
inhomogeneous parts which remain are now known functions of the time, and may be computed 
from the first order solutions. The equations can be integrated by elementary processes. 

The final expressions for r, z, and @, complete through second order, are 


br? (B—a) (48 —a—3) 6s? (3a8 —B—a?) , (48- 3a+1) 


se -o+|ir+—; Maa” aia os Guaticia’°" ii-an oe 


+—{ “es ér te) | inw(1— a) 
i —_ 's sinw = @ 
ia ”6=6huSi<@ ” =—CU= Ee 











|eosa(t —a)' 


+ aaa 26+1)( +) + 28- 2)( bs?-+— ~*~ 


(a +28) at,’ 
( bs?— 


(a+28— >/ a*¢,* 
——————[ 6r°— ———- 
4a(1—Sa) 


) cos2wa't 
12a(1—a) 


) cos2w(1 — a) #+ 


1—a a 


(a+2£8 —3) (a+28) ; 
—_——-¢,6r sin2w(1 — a) *¢-+——————¢, 4, sin2wa't, (12) 
6(1—a)! 2(1—Sa)a! 


28 / éréz 2a* ores 
)] coswa}t 
5a-—1 ~ Sa?—6a-+1 


1 gr (3a—1) 
of ~ eis) | sinwa't 
al 1—Sa 5Sa?—6a+1 
8 sinw[ (1 — a) at ¢rdz | 8 sinw[ (1 — a)*— at li ¢rbz -2] 
a—1—2(a—a*)! L(1—a)! ai a—1+2(a—a*)! L(i—a)? ai 























B cosw[ (1 —a)*+-a* iy drdz “ee B cosw[ (1 —a)'—at jiy drdz 1eres - 
a—1—2(a—a’)! La (a—a?*)! a—1+2(a—a’)! La ash ws) 
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(28B—a—a’*)/ér? 9? (a?—28) /é2?  ¢,” 
Pete 2)3eo 
4(1—a) a’? 1-—a 4(i—a)\a’®? a 
1 br br? (B—a) (48—a—3) 62? (3a8—B—a’) 

Seam 
(1—a)! 








@ a 3(1—a) 6(1—a)? a? (1—a)(1—Sa) 


sinw(1—a)*t-+——— 
1—a 





. (48 —3a+1) 1 gor (a+28 —3) $262 (a+28) 
¢r— ae 


** 9(1—a)(1—5a) .s a@48  « Ten oom 


(3a®—13a—28-+12)rar?—g? (Sa?+28) fas? ¢,? 
a | sinza (ta) 
M(i-a) le i—« 


8ai(1—5a) 
(3a?—13a—28+12) ¢,6r (Sa?+28) 9,62 
- cos2w(1—a)#¢-+-——_ cos2watt 
12(1—a)? a 4a(1—Sa) a 
¢r (a?—3a+28) g-dr (a?+28) ¢,d2 


ana 4(1—a)? a 4a(i—a) a 





|sin2.a' 


a? a 








(14) 


The integration constants have been chosen so that @, r, and z satisfy the initial conditions given 
earlier, namely, at ¢=0 


r=a+ér, @0=0, z=62, *f=awy,=v¢,, 2=dwye,=’¢e, P=P+r°F+2. 


It should be noted that in the orbital Eqs. (12), (13), and (14) those terms which do not contain 
ér, 62, y, or y, are the zeroth-order terms, while those which are linear and quadratic in these quanti- 
ties represent the first and second order terms respectively. 


3. THE SEMICIRCULAR SPECTROMETER 


When a=6=0, Eqs. (12), (13), and (14) describe the electron orbits in a semicircular spectrometer, 
and it will be fruitful to consider this case in some detail before proceeding further. Equations (12) 


and (13) become 


a 1 
r= a-+| ér — 3(e - o)| coswt+ ¢,(a+ dr) sinwt — — —a’¢g,*) cos2wt 
a 


ag,’ 


1 
—49,6r sin2wt-+—(ér?+a"¢,?) -—, 
4a 2 


z=d62+agwt. 


When the magnetic field is uniform, as it is here, there exists no restoring force in the z direction, 
and consequently there can be no z focusing. However, the first order sinusoidal terms in r show 
that r focusing will occur when wt=7, and that the image will be reversed in r with respect to the 
source. Let 6*, r*, z* be the coordinates of a focused electron, where 6* = x. Let ¢* be the time required 
for an electron to travel from the source to the focused position. In zeroth order 


6=wt, 0 =wt* =z, (17) 


while in first order 
or 
6 = wt —— sinwt+ 9, Coswt — ¢,, 
a 


F=xr=wi*—2¢,, wi*=2r+2¢,. 
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It will not be necessary for the present purpose to obtain the second-order correction to ¢*. Substi- 
tuting ¢* from (17) and (18) into the second- and first-order terms, respectively, of (15) and (16), 
r* and s* are found to be 






r*=a—br—a(o7+¢,), (19) 
2* = 62+ rag, t+ 2agr¢s. (20) 


It is evident that the second-order terms in Eq. (19) represent a departure from perfect focusing 
which is strongly dependent on ¢, and ¢,, and this will be discussed later. Since the orbits are all 
circular in form, it is clear that rigorous expressions for r* and z* may be derived from purely geo- 
metric considerations. When these exact expressions are developed through second order, they 
yield Eqs. (19) and (20), as they must. 
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4. THE DOUBLE FOCUSING SPECTROMETER 





When the first-order sinusoidal frequencies in Eqs. (12) and (13) are made precisely equal, a 
simultaneous focusing occurs in both the r and z directions. The necessary condition for this two- 
dimensional focusing is 






(14) 


w(l—a)ji=wal, a=}, 





In this case, (12) and (13) become 









































} given 
(26 —1) (88 —7) (1—88) wt 
rood ort - ee oa 
a 
contain (5—48) (1+48) _wt (3-46) (46 —1) 
juanti- +2 ae-+ ¢r5r+ oi] aa Tel (6r?+-2a?9,?) + (dz? +-2a*y,?) 
1 
+—[(48—5)(6r?—2a*¢y,?) — (1448) (dz? —2a*y,*) ] cos2twt 
meter, ide 
1s (12) 2! ' 
+-1e- 5) ¢-ér — (1+48)¢,62] sin2%wt, (21) 
48 168 , wt 48 _ wt 
s=| is iris ——anr, cos--+24| apt —(o.br+ oc) | sin— 
3a 3 J 28 3 24 
(15) r éréz 26f bréz 218 ; 
+28 —+2a¢,¢,|-—|-=— 2ae.e.| cos2!wi——[_g,6r-+ ¢,62] sin24wt. (22) 
L a 3La 3 

(16) 

f The first order sinusoidal terms in (21) and (22), namely those terms which are linear in dr, 4, ¢, 
— or yz, show that focusing will take place when wt/2!'=, and that the image will be reversed in r 
‘a re and inverted in 2. Let 6*, r*, 2* and ¢* be defined as before. Then, in zeroth approximation, 
quired O=wt, O* =wt* =2ie~254° 33’, (23) 

while in first order 

(17) dr = wt wt 

6=wt — 2'— sin—+2¢, cos——2¢,, 
 » * 24 
(24) 





O* =2r=wi*—4y,, wl* =2'r+4¢,. 





Substituting ¢* from Eqs. (23) and (24) into the second- and first-order terms, respectively, of 
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(21) and (22), r* and z* are found to be 
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As before, the second-order terms in (25) and 

(26) represent departures from perfect focusing 

in r and 32, respectively. The arbitrary constant 

8, which determines the second-order correction 

to the magnetic field, still remains as a disposable 

parameter and can be chosen so as to minimize 
. the second order defocusing. 

For the spectrometer of Siegbahn and Svart- 
holm, it is stated that the magnetic field was 
shaped to fall off as 1/r*. This is identical with 
Eq. (1) if a=4$ and B=}. When B=}, Eqs. (25) 
and (26) become 


, (6r?—627) 4 
r* =a — br +———_—_—-a(¢,?—,"), (27) 
6a 3 


éréz 
3* = —§s-+——. (28) 
a 


As in the semicircular spectrometer, the second 
order defocusing in r is strongly dependent upon 
¢- and g,, and this point will be discussed in 
greater detail later. 

With a different choice of the field parameter 
8, it is possible to eliminate ¢g, and ¢, entirely 
from the second order defocusing in r. Thus, 
when B=}, Eqs. (25) and (26) become 


(6r? — 62”) 


r* =a — br -+—___, (29) 
2a 
éréz 8 

s* = — 6s ++—— —-a¢,¢. (30) 
3a 3 


Comparing Eqs. (19) and (20), (27) and (28), 
and (29) and (30), the following conclusions may 
be drawn. It will be assumed throughout that 
the source is rectilinear in shape, with its long 
dimension parallel to the z axis. 

(a) For the semicircular spectrometer, the 
negative term a(g+ ¢,’) in Eq. (19) predicts 
that the center of gravity of the image will be 


a 


(328—12) 
2* = — §s-+-— —_+——"4 ¢- ¢. (26) 





3 a(¢r?— ¢,*), (25) 























a 
shifted to smaller values of r, and that the image 
will be correspondingly broadened. 

(b) For the double focusing spectromete 
where 6 =§, the term 4a(g,?— ¢,*)/3 in Eq, (27) 
may be either positive or negative. The j 
will be broadened toward both larger and smaller 
values of r. Its center of gravity may or may not 
shift; the magnitude and direction of the shift 
depend upon the relative maximum values of », 
and ¢,. The term. (ér*— 6z*)/6a in Eq. (27) pre. 
dicts that the image will be slightly concaye 
toward the origin of coordinates. 

(c) For the double focusing spectromete 
where 8 =4, Eq. (29) predicts a curvature of the 
image similar to that discussed under (b) above. 
However, the absence of yg, and gy, means that 
the image will be sharp, and its center of gravity 
is at r=a. 

(d) Terms of similar order of magnitude in- 
volving ay,y, appear in each of the 2* Eqs. (20) 
and (30), but not in (28). In addition, Eq. (20), 
for the semicircular spectrometer, contains the 
much stronger term rag,. Such terms mean that 
the image will be lengthened vertically, and this 
effect is especially prominent in the semicircular 
type. 

The conclusions listed under (a) through (@) 
above are illustrated graphically in Fig. 1. The 
rectilinear source has been assumed to have the 
dimensions 20 mm X2 mm, while the radius 
has been arbitrarily chosen as a=10 cm. Its 
further assumed that the angles ¢, and ¢, are 
limited by a system of baffles so that (9+ 9) 
=0.1 radian. Figure 1 illustrates the image 
formed in each of the three spectrometer types 
which have been discussed; only the upper hall 
of each image is shown, since the image i 
symmetrical with respect to the median plane 
The shaded fectangle represents a perfect image 
of the source, centered at r=a; it could also be 
considered as representing an exit slit, through 
which the focused beam passes into a Geiger 
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counter. The dotted contours represent the ex- 
treme boundaries of the focused image; it should 
be emphasized that the intensity distribution 
within the boundary is not uniform. Neverthe- 
less, it is evident that a large proportion of the 
focused particles in a semicircular spectrometer 
miss the exit slit entirely, and consequently the 
counting rate may be expected to be relatively 
low. On the other hand, the double focusing 
spectrometer, especially when the field parameter 
is =}, makes very efficient use of the focused 
beam, and the counting rate will be relatively 
high. Moreover, when 8 =}, the sharpness of the 
image leads to a considerably higher resolving 
power. If the exit slit is properly curved and 
elongated to fit the image, maximum counting 
rate and maximum resolving power will be 
attained simultaneously. 


5. CONCLUSION 


The analysis has shown that, through the use 
of a shaped magnetic field, the beta-ray spec- 
trometer proposed by Siegbahn and Svartholm 
may be made into a very effective instrument. 
The resolving power as well as the intensity of 
the focused beam are much higher than those of 
a semicircular spectrometer. Some technical diffi- 
culties may be encountered in the construction 
of the instrument, however, since it is necessary 
that the shape of the magnetic field shall remain 
independent of the magnitude of the field over 
the useful range of operation. 

With the shaped field where 8=} and where a 
curved exit slit is used, the resolving power is 
very high; in fact, the equations show that the 
image will almost precisely coincide with the slit. 
This nearly perfect predicted performance may 
be slightly illusory in the sense that third and 
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Fic. 1. The focused images for a semicircular spec- 
trometer and for double focusing spectrometers with 
B=} and }, respectively. 


higher order terms in the development of the 
solution will again lead to some defocusing. It is 
clear, however, that these will be small and of 
the order of magnitude of, say, 5z*/a* or of ag,?. 
For the spectrometer dimensions chosen in the 
present example, these terms would have extreme 
values of about 0.1 mm. In principle, through 
the use of a further shaping constant y for the 
magnetic field, it might be possible to cancel out 
much of the third order defocusing. This would 
appear to be too refined a process to be practical 
however. 
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The Variational Method for Asymptotic Neutron Densities* 


R. E. MARSHAK 
University of Rochester, Rochester, New York 


(Received January 24, 1947) 


Inhomogeneous integral equations connected with a certain class of neutron diffusion prob- 
lems are treated by means of a variational method. It is shown how to obtain quite accurate 
values of the asymptotic neutron density in the following three cases: (a) Milne’s problem for 
the plane, (b) infinite scattering medium surrounding black spherical core, (c) infinite scattering 
medium surrounding black spherical core with air gap. 





¢ i ‘HE variational method for solving eigenvalue problems associated with integral equations js 
well known. One starts with the integral equation: 


po(t) =c f dr’ooe)K (4, ©), (1) 
where c is the (lowest) eigenvalue to be determined, po(r) is, say, the neutron density as a function 


of the coordinates, and K(r, r’) is the kernel (assumed to be symmetric in r and r’); the integration 
is over the medium in question. It can then be shown that the functional : 


f aro()| o(t)—c f dr’ p(x’)K(r, r)| (2) 


is an extremum! for the actual solution po of Eq. (1); furthermore, it is evident that: 


f dt po*(r) 


"feats fate r’) | 





c 


Since the first variation of (2) vanishes, we may expect that the value predicted for the eigenvalue ¢ 
by the insertion of some trial function for the unknown pp will be closer to the correct value of ¢ 
than the trial function is to po. It can be proved that the approximate eigenvalue is always larger 
than the correct eigenvalue. 

The extension of the variational method to inhomogeneous integral equations is possible and leads 
to a simple procedure for obtaining fairly accurate solutions for a certain class of neutron diffusion 
problems. We shall first present briefly the method in general terms and then discuss specific appli- 
cations. Suppose the inhomogeneous integral equation is:? 


go(x) = f dx'K (x, x’)qo(x’) +f (2), (4) 


where go(x) is the unknown function, K(x, x’) is a positive symmetric kernel, and f(x) is bounded so 
that {| f(x) |dx exists. Suppose, further, that we can express [ fdxgo(x) ] in terms of go( ©), say 


ge(2o) =k f daego(x) f(x) +he 


* This work was done in the Montreal Laboratory of the National Research Council of Canada during the Spring of 
1944 and is now declassified. : 

1 This can be seen easily by writing the trial function p(r) as [po(r)+<(r) ],(e(r) is the correction function), and showing 
that the term of first order in « vanishes identically. 

? The method is presented for a function of one variable; the generalization to more than one variable, as in the cast 
of the homogeneous integral equations, is also possible. 
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with k:, &: known constants. Let us define the functional: 


fsa) | ace - favK(, wae] 


f axa(ws(e)| 


The reasoning then proceeds along the same lines as in the homogeneous case: (5) is an extremum® 
for the actual solution, go(x), of Eq. (4), and can be equated to [k1/(go( ©) —k2) ]. For an arbitrary 
trial function, g(x), the equation determining go() is: 





(5) 


f dxa(x){ a) - fave, va(e)| 


. iy o( 00 ) lbp 
f axa(a)fc) . 


We shall now show how several neutron diffusion problems requiring a knowledge of the asymptotic 
neutron density can be written in the form (4). The application of (6) will follow immediately. 

We consider first the well-known Milne problem, i.e., the case of a semi-infinite, isotropically 
scattering (but non-capturing) medium bounded by vacuum which sustains a constant current from 
infinity. A rigorous expression for the asymptotic neutron density is known for this problem;‘ this 
is not the case for the other two problems discussed below where at present only approximate solu- 
tions are available. The application of the variational method to the Milne problem will disclose 
the manner in which values of the asymptotic neutron density are arrived at and give an indication 
of the power of the method. 

The integral equation for the neutron density, po(x) (the boundary is taken at x =0 and the scat- 
tering mean free path is taken as the unit of length), in the Milne problem is: 


(6) 





pols) =3 f dx'E,(|x—x'|) pox’), (7) 


where E;(x) is the exponential integral function of order one.’ Since the asymptotic neutron density 
is linear, we write po(x)=x-+qo(x); for large x, go(x) approaches a constant which we denote by 
go(). Substitution for po(x) into (7) yields: 


cs) E;( ) 
a(x)=3 f dx’ Ey(|x—x’ | )go(x’) + 2 ° (8) 





2 


where E;(x) is the exponential integral function of order three. Equation (8) is of the form (4) but 
we still require an expression for f°dxqo(x)E3(x) in terms of go() before applying (6). 
To obtain such a relation,* we make use of a theorem proved by Davison,’ namely that if we have 


* Dr. B. Davison, in the following article, discusses the conditions under which (5) is a minimum, not only an extremum, 
for the actual solution. 

*E. Hopf, Mathematical Problems of Radiative Equilibrium, Cambridge Tracts No. 31 (1934). 

* The exponential integral function of order nm is: E,(x) = (i°[e~**}dt/t. 

*In the original report, the relation between go(#) and Jo dxgo(x)f(x) in all three problems considered, was derived 
by another method. It is more elegant to make use of Davison’s theorem. 
hs The proof of this theorem is simply as follows: Q(r) in (9) can be regarded as r times the neutron density in an infinite 
(isotropically) scattering medium due to a spherically symmetrical system of continuously distributed sources, the 
source density at distance r from the origin being F(r)/r. The tota yy of these sources will be /dr(F(r)/r) 
=4r f{.°rdr F(r) and Q() is equal to 3/4x times the output of all sources; this leads to (10). 
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an integral equation of the form: 
0) = [ar Oe ILBA(Ir—) —Eilr ++ FO) 
and Q(r) is its solution bounded at infinity, then 
Q(«)=3f “rdrF(r). 
0 
A comparison of Eqs. (8) and (9) yields: 


go( © ) -if xdxEa(x) +4 f xdx f dx' qo(x')Ex(x+<x’). (11) 
0 0 0 


Interchanging the order of integration in the second term on the right-hand-side (r.h.s.) of (11) an 
evaluating the resulting integral and the first term on the r.h.s., we get: 


aul) =3 f dgo(x)Ba(x) +8. (1) 
0 
Equation (12) is of the desired form and Eq. (6) becomes: 


ff avace ax) —3 f dx!q(e)Ex(|x—x'|)| 
0 0 


ty f ” o(0)—} 
i| f dxq(x)Ex() | ig 


Equation (13) is the final result, and from it we may calculate go( ©) for different choices of the trial 
function g(x). The procedure now is to assume a function for g(x), and choose the parameters in this 
function so that the left-hand side (I.h.s.) of (10) is an extremum; the resulting function is then 
used to calculate go(). In this process the trial function g(x) is determined except for a constant 
factor; this factor may be taken so that g(x) is asymptotic to the value obtained for go(). 

As has already been remarked, the whole point of the variational method is that the choice of 
q(x) may be quite rough and yet give a fairly accurate value for go( ©). Thus, in the present instance, 
even if we assume g(x) =constant (for which no variation is necessary in the left-hand side of (13)), 


| f “dxkx(s)] 
f as] 1-3 f dx! Ex(\x—2'|)] 


=17/24=0.7083, 





we get: 





qo( ©) = #4 








which is only 0.3 percent less than the correct value 0.7104.* The result can be considerably improvel 
by starting from a more elaborate trial function and actually using the extremum property of the 
functional. In fact, LeCaine has shown that by choosing a trial function which simulates mor 


closely the the actual g(x), namely: 
q(x) =const.[1—A E(x) +BE;(x) ] 


* Cf. E. Hopf, reference 4; the value 0.7083 should be less than the true value (which it is!) since the I.h.s. of (13) 
should give a minimum for the correct solution of the integral equation (8) (cf. reference 3), 
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and extremizing the left-hand side of (13) with respect to A and B, one arrives at a value of go( ©) 
which agrees with the correct value to six decimal places (0.7104457 as compared to the correct value 
0.7104461). 

The second application of the variational method considered here is an estimate of the asymptotic 
neutron density in an infinite, isotropically scattering (but non-capturing) medium supporting a 
constant flux from infinity and surrounding a perfectly absorbing spherical core. This problem was 
treated by the spherical harmonic method,’ and we shall compare the results of the variational 
method with those results. The integral equation for the neutron density in the medium surrounding 
the black core can easily be derived; we find: 


po(r) = rf dr’ po(r’) {Ex(|r—1|) — Ei ((r? —a*)#+-(r’?—a*))}, (15) 


where po(r) is r times the neutron density and a is the radius of the black core. 
We are interested in a solution of (15) which is asymptotically linear, and we write this solution as 


po(r) =r—Q(r). (16) 


The existence of this solution can be established in the same way as the existence of the solution of 
Eq. (7). On substituting (16) into (15), we find the integral equation for Q(r), namely: 


Q(r) = f dr’ Q(r”)K(r, 1’) +e(r) (17) 


K(r, 1’) =3{Ei(|r—r’|) —Ex((r?—a*)!+ (7’2—a*)))}, (17a) 

es(r) = 4a E,(r —a) + E;((r?—a*)*) —E;(r —a) J. (17b) 

The kernel K(r,r’) is obviously, symmetric; it is also positive, as can be seen by noting that 

|r—r’| <[(r?—a?)*+ (r’?—a*)*] and that E(x) is a monotonically decreasing function of x. Moreover, 

the function ¢3(7) is bounded and the integral /,*| €s(r) |dr exists. Hence, the variational method can 
be employed to estimate Q() provided Q(«) can be expressed in terms of /,°drQ(r)es(r). 

A relation between Q() and JS,*drQ(r)es(r) can be found by again making use of Davison’s 


theorem (cf. Eqs. (9) and (10)). For this purpose, we must define Q(r) and ¢;(7) in the interval 
0<r<a; we choose:!® 


on=4 f dr'Q(r'){Ex(\r—r'|)—Exlr+r’)}, (O<r<a), (18a) 


e;(r) =0, (O<r<a). (18b) 
Further, we write: 


AQ(r) =4 f dr’ Q(r’){Ex(|r—r’|) —Exlr+r’)}, (19a) 
0 


A:0(r) =4 f dr’ Q(r’){Ex{r—r'|) -Ex(r+r)}, (19b) 
0 


an0(r)=4 f dr'Q(r’) { Ex((r? —a*)#+-(r’?—a*)'!)-Ex(r+r’)}, (ra), (19c) 


=0, (0<r<a). (19d) 


*R. E. Marshak, Phys. Rev. 71, 688 (1947). 
_ Cf. Davison’s declassified Montreal report MT-232, “Influence of an Air Gap Surrounding a Small Black Sphere 
Upon the Linear Extrapolation Length of the Neutron Density in the Surrounding Medium.” 
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In view of the above definitions, Eq. (17) becomes: 
Q(r) = AQ(r) —AiO(r) — A20(r) +63(7r). (20) 
Equation (20) is also satisfied for 0< r<a. If we now use Eq. (9), we get: 


Q() = 3| ff rrsin =f rarcas0091 $y ff rans. . (21) 


The evaluation of the first term on the r.h.s. of (21) is straightforward; the result is 3a?/4, The 
second term on the r.h.s. can be evaluated by interchanging the order of integration to give 
—3S,*rdrQ(r). However, Q(r) for 0<r<a is defined by (18a); inserting (18a) and interchanging the 
order of integration a second time yields: 


=a ar ow) f rartB\r—r) Burt}. 
a 0 
The third term on the r.h.s. becomes, on interchanging the order of integration: 
=a ar'0e) f rant Bi((—a4 (72-049) —Br41)}. 
It is easy to show that: 
aff rdrBa(lr—r'l) +f rdrB*—0)+ (72-09) — f rire +r) | =e;(r’). (22) 
0 a 0 
Hence, Eq. (21) leads to the desired relation : 
Q(2) -— 3 f “drO(nalr). (23) 
Substituting (23) into (6) yields an equation corresponding to (13), namely: 


f “ara ow) -} J “ar Q")K(, 2) 











= ; (24) 
«© 2 a? O( 0 ) 
[ f 2a | — 
a 4 3 

Again, as in the plane Milne case, we choose the simplest trial function, i.e. Q(r)=const.; we get 

ai - p 

s|-——+ f arEy(r—<"))| 

3a? ee ae A 

Q( 2) ay ay ; (25) 


[i+2 f * f “ar E04 (7-0) | 


Values of [a—Q()]—i.e., the ‘extrapolated endpoint”—as a function of a computed on the basis 
of (25) are given in column 2 of Table I, while the values of the extrapolated endpoint predicted by 
the P;-approximation of spherical harmonic method? are given in column 3. It is seen how powerful 
the variational method is for the determination of the asymptotic neutron density. Of course, 4 
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slightly more accurate choice of trial function, TaBLe I. Extrapolation endpoint as function of a. 
which takes into account the exponential-type 
approach of Q(7) to Q(«), would lead to con- 








a Variational method Ps-approximation 





siderably more accurate values of the extrapo- y 0.347 0.334 
jated endpoint. — ge oe that 20 426 = 

1 solution o maximizes the ex- 5. .688 690 
the actua uM -— 690 








trapolated endpoint as is to be expected from 


Davison’s paper.’ 
Finally, we have applied the variational method developed above to the case of a black sphere 


(of radius a) surrounded by an “air gap’’ (of outer radius 6) which in turn is surrounded by an 
isotropically scattering (but non-capturing) medium supporting a constant flux from infinity. The 
integral equation for r times the neutron density turns out to be: 


oo(r) =4 f dr’ polr')R(r, 7’) (26) 


where: 
R(r, r’)= E\(|r-r’|) — E;((r? —b?)'+ (r’2 — 6b?) t) 


1 
+ — exp] —R+ (UR =r’) —4R—b9)I |, 
a R R 
with 

a=(r*—b*)14(r'2—b4), B= (r*—a*) + (r’2—a")), 


As in the two preceding cases, since capture is absent in the outside medium, the asymptotic solution 
is linear. Following the procedure for the black sphere without gap, we obtain the integral equation 


for Q(r), (po(r) =r—Q(r)), namely: 
Q(r) =4 f dr’ Or) Rr, 7’) alr), (27) 


where: 
és(r) = $[bE2(r —b) — Es(r —b) +Es((r? —a*)* — (b? —a*)!) — (b? —a*) 4 ((r? —a*)*— (6?—a?)#) ]. (27a) 


We find for Q(«) (cf. Eq. (23)): 


Q(~) -—-3f ara) és(r). (28) 


Substitution of (28) into Eq. (6) then permits an immediate application of the variational method 
to estimate Q(«). No numerical results have been obtained for this case. 

It isa pleasure to thank Professor J. Schwinger for a valuable conversation on variational methods, 
Dr. Davison’s cooperation is also appreciated. 
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The nature of the extremum of the functional 
S ace) eer fa@)K vay ]ax /| faceyoceras] 
associated with the inhomogeneous integral equation, 
go(x) =A f goly)K (x,y)dy+fole), 


is investigated. Two sets of sufficient conditions for the extremum reached for g(x)=go(x) to 
be a minimum (and neither a maximum, nor a “saddle point’’) are determined. 





1, INTRODUCTION 


R. R. E. MARSHAK in the article on 
“The Variational Method for Asymptotic 
Neutron Densities,’’ appearing in the same issue 
of this journal, solves certain inhomogeneous 
integral equations by use of the following theo- 
rem: 
Let go(x) be the (bounded) solution of the 
integral equation 


go(x) —d f K(x,y)go(y)dy—fo(x)=0, (1) 


in which the parameter X is positive and its 

value is known beforehand, the kernel K(x,y) is 

positive and symmetrical and the free term fo(x) 

is bounded, and the integral | fo(x) | dx exists. 
Then the functional 


f a(a)| ace) a f K(xy)a(o)dy fs 


i= » (2) 


f a(e)folxddr| 


in which g(x) is any bounded real well-behaved 
function, reaches its minimum for g(x) =qo(x), 
provided either 


af K(x,y)dy <1 (3) 





for any fixed x, or 
A<Ao (4) 


in which Xo is the limit, for 5, of the lowest 


eigenvalue of the equation 


o(x) =)’ f K(x,y) o(y)dy, (5) 


and provided K(x,y) satisfies certain general 
qualitative conditions, which will be formulated 
later. 

The purpose of the present note is to prove 
this theorem. 

(Notice that the conditions (3) and (4) are not 
equivalent, but represent two different types of 
cases in which one can use the functional (2); 
There are cases when \=Xo but (3) is satisfied, 
and there are cases when (4) is satisfied but (3) 
is violated.) 

We shall split the proof of the theorem into 
two parts. In the first part we prove the lemma: 

The functional (2) reaches its minimum for 
q(x) =qo(x), if X and K(x,y) are such that 


f ae) ac0)-r f K(x)a(0)dy fs 6 


is positive for any bounded, not identically 
vanishing, real, well-behaved g(x), and that 


f aerttarar= ff aoerfcerde 


0 0 


in which f(x) is any function satisfying the same 
assumptions as fo(x) and q:(x) is the solution of 
the Eq. (1) with fo(x) replaced by fi(x). 

In the second part, under some general quali- 
tative assumptions with respect to K(x,y) we 
prove that (6) is positive for any bounded, not 
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identically vanishing, real, well-behaved @(x), 
if either (3) or (4) is satisfied. 


2. PROOF OF THE LEMMA 


Let g(x) be some trial function which we 
substitute into the functional (2). Present it in 


the form 
q(x) =go(x) +€q:(x). (8) 


Let f:(x) be the value of the left-hand side of 
(1) when we replace there go(x) by g(x) (i.e. 
qi(x) is the solution of (1) with fo(x) replaced by 
fi(x)), and substitute (8) into (2). Then, after 
a slight rearrangement, and on using (7), we get 

1 € 2 


I=-— 1+ 
f go(x) fo(x)dx 








- al, @) 
f aerfeerar 


in which 


f fo(x)go(x)dx f fo(x)qu(x)dx 


A= . M , (10) 
ff s@aciae J A@acar 
0 0 








or, which is the same as, 


A= f fo(x)qo(x)dx 


f Cfa(x) —ufo(x) TLar() —ngo(x) Idx, (11) 


in which » is a constant determined from the 
condition 


J go(x) Lfi(x) — ufo(x) Jdx =0. 


Put gi(x) —ugo(x) = g(x), and notice that thereby 
fils) —ufolx) =4(2) af K(w.9)G(0)dy. 


Thus if (6) is necessarily positive then A is 
positive, and the expression (9) reaches its mini- 
mum for e=0, i.e., the functional (2) reaches its 
minimum for g(x) =go(x), which we wanted to 


prove. 
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Notice that if \ and K(x,y) are such that (6) 
can have either sign, depending on the actual 
form of G(x), then, assuming /(*fo(x)go(x)dx to 
be positive, A can have either sign depending on 
what trial function we used, and (9) will have 
neither a minimum nor a maximum, but a 
“saddle-point”’ for e=0, i.e., the functional (2) 
will have neither a minimum nor a maximum 
for g(x) =qgo(x), but a “‘saddle-point.”’ 


3. CASE WHEN 4<2A. 


Turn now to the determination of conditions 
under which (6) is positive for any not identically 
vanishing, real, well-behaved function G(x), and 
let us deal first with the case \ <p. 

Let us examine first the case of an equation 
with a finite range of integration, i.e., the case 
when K(x,y) and f(x) vanish for x>b, so that 
the homogeneous equation corresponding to our 
inhomogeneous has the form (5). Let Ao’, Ay’, 
etc., be the complete set of eigenvalues for the 
Eq. (5) and go(x), ¢:(x), etc., be the corre- 
sponding eigenfunctions. It is known that since 
K(x,y) is symmetrical, all ),,’ should be real and 
the ¢m(x) should be mutually orthogonal.! We 
shall assume further that K(x,y) is such that all 
Am’ are positive and ¢,,(x) form a complete set of 
orthogonal functions. Let Ao’ be smallest of all 
Am’ (so that lims../Ao’ =Ao). Once we have as- 
sumed that ¢,,(x) form a complete set of orthog- 
onal functions, we can write 


G(x) = aaya(x) (12) 
n=0 


and hence? 


’ 


» oe A, —A 
a(s)—r f K(xa)40)dy=E ox eal), 


n 





so that 
b b 
f ae) a00)--r f K(s4)dy fx 
0 0 
» dv ot 
= 2 an" x J Gn*(x)dx. (13) 


1See e.g., Whittaker and Watson, Modern Analysis 
(Cambridge University Press, Cambridge, England, 1927), 
pp. 225-227. 

* See reference 1, p. 228. 
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If X\<Ao, then for sufficiently large values of 
we shall also have \<Ao’, hence all the terms on 
the right-hand side of (13) will be positive. Then 
the left-hand side of (13) should be positive too, 
and passing to the limit with b> we see that 
(6) is positive too. 

If A>Ao, one can see from the same argument 
that (13) and consequently (6) will be negative 
if G(x) is such that the term in (12) corresponding 
to the lowest eigenvalue [or the lower eigen- 
values] strongly predominates over the terms 
corresponding to high eigenvalues, and positive 
in the opposite case, so that for A>Xpo, the 
quantity (6) can have either sign depending upon 
the form of G(x), and hence the functional (2) 
will have merely a saddle-point for g(x) =qo(x). 
The case of A=Xo should be considered sepa- 
rately, and this will be done in the next section. 

Insofar as the assumption (7) used in proving 
the lemma is concerned, it may be pointed out 
that if X\<Ao and if K(x,y) is such that, in the 
integral 


f ee ff) K (es) K (6459) fi 
XK(Sn,y)fily)dxds,---dsydy, (14) 


we can change the order of integration in any 
arbitrary manner, then (7) is automatically 
satisfied. For if \<Apo, then go(x) and qi(x) can 
be determined by means of Neumann’s pro- 
cedure, once fo(x) and f:(x) are known, i.e., 


go(x) =fol(x)-+X f K(x,y)foly)dy 


+f Jf KeK(soyfolo)drds +ete., (15) 


and a similar expression for g.(x), the series 
appearing in both these expansions being abso- 
lutely and uniformly convergent. And then from 
the possibility of interchanging the order of 
integration in (14), (7) follows at once. 


4. CASE WHEN 4=20 


In this case we can prove that (6) is positive 
. for any bounded, not identically vanishing, real, 
‘ well-behaved g(x), and consequently (2) reaches 
its minimum for g(x)=qo(x), provided (3) is 


satisfied, and provided \K(x,y) satisfies the fol. 
lowing additional assumptions:* 

(i) The homogeneous equation corresponding 
to (1) cannot have any bounded solutions (so 
that the bounded solution of the inhomogeneoys 
equation is necessarily unique) 

(ii) f(x) is bounded and the integral /¢*| f(x)|dy 
exists 

(iii) For any f(x) satisfying (ii) the series of 
the form (15) converges for any particular valye 
of x and its sum is uniformly bounded with 
respect to x; in other words, using Hopf’s termij- 
nology,‘ the Eq. (1) possesses a unique bounded 
N-solution for any f(x) satisfying (ii) 

(iv) If f:i(x) and f2(x) satisfy (ii), and g(x) is 
the solution of (1) with go(x) replaced by q,(x), 
then in the integrals 


is f KerKes me 


XK (Sn,y)q2(y)dydsn---dsdx, (14’) 


we can interchange the order of integration in 
any arbitrary manner. 

Notice that (7) follows at once from our equa- 
tions and the assumption (iv). Notice also the 
following corollaries of assumptions (i) and (iii). 

(A) If the free term is non-negative, then the 
corresponding bounded solution is also non- 
negative. (This follows since every term in the 
series of the form (15) is non-negative, conse- 
quently its sum is also non-negative. On the 
other hand this sum is obviously a solution, and 
according to (i) the solution is unique.) 

(B) Put for brevity 


Af(x) =) f K(x,y)f(y)dy; 


(16) 
AA"—!f(x) = A*f(x). 
Then 
limn+eA"f(x)=0 (x fixed). (17) 


*Since our purpose is merely to give a rigorous justi- 
fication of the method of solution of certain problems given 
in Marshak’s paper, we shall not attempt to prove our 
result under the most general assumptions ible, nor 
to make our assumptions mutually independent. But we 
shall assume such properties of K(x,y) as it does possess in 
the examples treated by Marshak, and shall assume all the 
required properties independently, irrespective of the fact 
that some of them may be corollaries of the others. __ 

*E. Hopf, Mathematical Problems of Radiative alib- 
rium (Cambridge Tracts on Mathematics and Mathe- 
matical Physics, 1934). 
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(This is merely the necessary condition for the 
convergence of the series (15) with fo(x) replaced 
by f(x).) 

It should be mentioned that the only way in 
which we shall use the assumption (iii) is through 
its corollaries (A) and (B). It should also be 
mentioned with respect to (B) that we shall not 
use this corollary in full, but merely a somewhat 
weaker and clearer property, namely: 

(B’) Either one of the consecutive terms in the 
series (15) with fo(x) replaced by f(x) becomes 
non-negative (or non-positive), or for any par- 
ticular value of x the consecutive terms of this 
series tend to zero (i.e., we shall need to use (17) 
only for the cases when none of the functions 
A*f(x) conserves its sign throughout the range of 
x-values). 

Having clarified our assumptions, turn now to 
the proof that (6) is necessarily positive. Put 


fix) =a(x)—-) f K(xy)G(v)dy, (18) 
and : 


N(f)= f Fix)a(x)dx, (19) 


and notice that f(x) =f1(x) —ufo(x) satisfies the 
assumption (ii) and hence, according to (3), so 
do all A*f(x). Using the assumption (iv), and 
noticing that G(x) is the bounded solution of (1) 
with fo(x) replaced by f(x) and Ag(x) is the 
bounded solution of (1) with fo(x) replaced by 
Af(x) we have 


N(f)= J . f*(x)dx+r J ; J "K(xy) 


X f(x) fly)dxdy+N(Af). (20) 


The first term on the right-hand side in (20) is 
obviously positive. Let us show that the absolute 
value of the second term is less than the first 





term. Use for this purpose Schwartz's inequality, 
which in two dimensions has the form 


f f F,?(x,y)dxdy - f f F,*(x,y)dxdy 
>| ff rccormteaddady. 


F,(x,y) and F2(x,y) are real. Taking here 
F,(x,y) =f(x)AK (x,y), Fa(x,y) =F(y)AK (x,y), 


interchanging in the integral {” / F;*(x,y)dxdy 
the variables x and y, and extracting the square 
root, we get 





d J i J “Keso)fte)Foddeds 


<X f f K (x,y) f?(x)dxdy. (21) 


On the other hand multiplying (3) by f(x) and 
integrating over all x, we see that the right-hand 
side of (21) is less than Jo*f*(x)dx. Then it 
follows from (20) that 


_ N(f)>N(Af), (22) 
and similarly, 


N(Af)>N(Af), ete. 


If one of the functions f(x), Af(x), A’f(x), etc., 
A™f(x), say, conserves its sign then N(A™f) is 
positive, and, according to (22), so is N(f). If 
none of the functions A*f(x) conserves its ‘sign, 
then, according to (B’), A*f(x) will tend to zero 
with increasing m, and N(A*f(x)) will also tend 
to zero with increasing m, and the formulae (22) 
then show again that N(f) is positive, i.e., the 
expression (6) is positive and the functional (2) 
reaches its minimum for g(x) = qo(x). 
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The True Height of an Ionospheric Layer 


J. A. PIERCE 
Cruft Laboratory, Harvard University, Cambridge, Massachusetts 


(Received November 4, 1946) 


A simple method has been developed for the analysis of ionospheric sweep-frequency records 
in terms of the scale height and true height of maximum of the layer of Chapman form which 
most closely fits the observed data. It is shown that the height of maximum determined by the 
method of Appleton, Booker, and Seaton is the height of a parabolic layer which is not coin- 

,cident with the parabola that most closely fits the Chapman distribution. The present method 
includes a correction for retardation in one or more layers of similar form which lie below the 





INTRODUCTION 


OOKER and Seaton! have published a 
discussion of a method for analyzing virtual 
height-frequency records in terms of the assump- 
tion of a parabolic distribution of densities, in 
order to determine the true height of the maxi- 
mum and the relative thickness of the dis- 
tribution curve. It is the purpose of the present 
paper to develop a similar treatment for iono- 
spheric data of the same type, which is based 
upon the assumption of a Chapman distribution. 
The assumed variation of the density of 
ionization with height in an ionospheric layer is 
shown by the solid curve of Fig. 1. The relative 
density is given by the equation 
N/No=([exp(a+1 —e*) }! 


where NV=ion density at a height 4, No=mazxi- 
mum ion density, e=base of the natural loga- 
rithms, a=(h4y—h)/H, hu=height at which 
N=Np, and H=scale height. The constants used 
for the curve shown are hy = 300 km and H=50 
km. Equation (1) is given by G. W. Pierce® and 
defines the Chapman distribution in ratio form, 
Its derivation is based on the concept of mono- 
chromatic ionizing radiation falling upon an 
atmosphere of constant temperature with a 
density that decreases exponentially with in- 
creasing height. The equation actually defines 
the rate of production of ions under these con- 
ditions and H is the height of the homogeneous 
atmosphere. As diffusion, among other factors, is 
not considered, we may not leap to the con- 
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'H. G. Booker and S. L. Seaton, Phys. Rev. 57, 87 (1940). 
2G. W. Pierce, Cruft Laboratory Lectures. 
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clusion that the total density of ionization at 
any moment is defined by this equation, but it 
seems reasonable to suppose that the sustained 
levels of ionization will be represented by a curve 
of this same general form. It will be assumed, for 
the time being, that the sustained density has 
exactly this form. It does not follow, however, 
that the scale height, H, for the density distri- 
bution is exactly the height of the homogeneous 


atmosphere. 
A parabolic distribution is defined by the 


equation, 

N hu—h\? 

oat es | , 
0 T 





where + is the half-thickness of the parabola 
(where V =0) and the other quantities have been 
defined under Eq. (1). It will be seen from Fig. 1 
that the Chapman distribution may be fitted 
quite closely by a parabola, when N/Np is 
greater than about 0.4. In comparing the para- 
bola with the Chapman curve, we find that if 
7=1.82H the parabola lies wholly within the 
other curve and is approximately tangent to it 
at the point where V/No=0.710. If r=2H, the 
parabola is coincident with the Chapman dis- 
tribution at the maximum and the curves cross 
at the point where N/N»)=0.206. The parabola 
shown is drawn for r=1.84H and is considered 
to be about the best over-all fit. 

Synthesis of a virtual height-frequency curve 
on the assumption of a parabolic distribution of 
density neglects the retardation due to passage 
of the rays through the lower part of the Chap- 
man distribution, a part that must certainly 
exist in the ionosphere. This retardation will be 
small, but not negligible, for frequencies near the 
critical frequency, and will be of considerable 
importance at the lower frequencies where the 
rays are reflected at lower heights. 


THE METHOD OF COMPUTATION 


For the calculation of the virtual heights we 
shall make the assumptions that propagation is 
in a region of no magnetic field (or that trans- 
mission is for the ordinary ray in a direction per- 
pendicular to the field) and that the group 
velocity of the wave is the product of the velocity 
of light and the index of refraction. As our com- 
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putation is to be made for the case of normal 
incidence upon the ionized region, the virtual 
height will be the integral of 1/u from h=0 to the 
height where »=0. Unfortunately a solution of 
this equation has been found for only a few dis- 
tributions of density of ionization which do not 
include that defined in Eq. (1). It is therefore 
necessary to adopt an approximate method. The 
most straightforward treatment appears to be 
the computation of the group velocity, or prefer- 
ably its reciprocal, for a large number of 
laminations each of which is assumed to have 
uniform density throughout its thickness. This 
computation must be made for each lamination 
at each of a fairly large number of frequencies in 
order to define a virtual height-frequency curve, 
but fortunately the formula for 1/4 may be put 
into a form suitable for a calculation of this sort 
so that a virtual height for a single frequency 
may be determined, for a case involving a 
hundred or more laminations, within twenty or 
thirty minutes after the density values for the 
various laminations have been found. 

The most convenient form for the equation is, 


1 -| (f/fo)? ] 
(f/fo)?—N/Nol 


where fo is the critical frequency or the frequency 
for which yu is zero when N = N». After tabulating 
values of N/N, for the various laminations in a 
vertical column, the values of 1/u can easily be 
worked out with a comptometer and set down in 
a parallel column which is preferably headed by 
the value of (f/fo)*. Since the value of 1/,y is 
defined by the ratio (N/No)/(f/fo)*, the com- 
putation may be shortened by preparing a table 
relating these two quantities. This table need not 
be carried beyond the point where 1/y is 2.5 or 
3.0 as the labor involved in carrying the table 
further would become greater than that saved by 
its use. 

Because 1/u goes to infinity at the point of 
reflection, special consideration must be given 
to the last few laminations. One method is to 
take thinner and thinner laminations as the 
point of reflection is approached, but this is quite 
laborious and only defers the question of how to 
treat the last layer of all. A rough method, which 
has no theoretical justification but which seems 
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RELATIVE FREQUENCY 


Fic. 2. Comparison between the laminar computation 
and an accurate method of calculating virtual heights. The 
data are for a parabolic distribution at a height of 300 km 
with a half-thickness of 92 km. 


accurate enough for practical use, has been 
developed by a study of the parabolic distri- 
bution for which virtual heights may be cal- 
culated precisely. This method is to compute the 
value of N/ No for a point eight-tenths of the way 
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from the top of the last complete lamination tp 
the point of reflection (where f?/fo? = .N/Ns) and 
to multiply the corresponding value of 1/x by the 
thickness of this last fractional lamination. Tha 
this quantity, when added to the sum of the 
tabulated values of 1/u for the complete lamina. 
tions, gives approximately the correct tota} 
virtual height is indicated by the fact that the 
treatment may alternatively be applied to the 
last one-and-a-fraction laminations instead of to 
the last fraction alone, with an alteration in the 
total virtual height which is generally not more 
than 0.02H. The errors of the method ar 
indicated in Fig. 2, where the solid curve is the 
relation between virtual height and relative fre. 
quency for a parabolic distribution of density 
(hu =300 km, r=92 km) calculated from the 
equation given by Appleton,’ and the circles are 
points computed as described above. It appears 
that the errors in these calculations, which were 
made for lamina of 1-km thickness, are small in 
comparison with the usual experimental error in 
measuring virtual heights. 


METHODS OF ANALYSIS 


For any law relating density of ionization to 
change in height it is necessary to calculate only 
one virtual height-frequency curve because the 


TABLE I. Values of (h’—hy)/H as a function of f/fo. 








h’ —hu h’ —hu 
H S/he H t/fo 


h’ —hu 
H 





— 2 0.22 — 1.92 0.44 
— 2.99 0.23 — 1.89 0.45 
— 2.83 0.24 — 1.86 0.46 
—2.73 0.25 — 1.84 0.47 
— 2.65 0.26 —1.81 0.48 
—2.58 0.27 — 1.78 0.49 
—2.52 0.28 —1.75 
—2.46 0.29 — 1.73 
—2.41 0.30 — 1.70 
—2.37 0.31 — 1.67 
—2.33 0.32 — 1.65 
—2.29 0.33 — 1.62 
—2.25 0.34 — 1,59 
—2.21 0.35 — 1.57 
—2.17 0.36 — 1.54 
—2.14 0.37 —1.51 
—2.10 0.38 — 1.49 
— 2.07 0.39 — 1.46 
— 2.04 0.40 — 1.43 
—2.01 0.41 —1.41 
— 1.98 0.42 — 1.38 
— 1.95 0.43 — 1.35 
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*E. V. Appleton, Proc. Roy. Soc. A162, 451 (1937). 
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Fic. 3. The relation i 
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increase of virtual over true height of reflection 
is proportional to the scale height, as is the dif- 
ference between true height of reflection and the 
height at which maximum density occurs. Thus 
the virtual height may be referred to the height 
of maximum density by the equation 


h'—hy 
i= hut | (4) 
H 


where h’ is the virtual height of reflection and the 
other quantities have been defined above. The 
quantity in brackets is a function of f/f» which 
Booker and Seaton! have called “y’’ when 
dealing with a parabolic region, but which we 
shall leave unnamed as it may assume any 
number of forms corresponding to various laws 
connecting the density of ionization with height. 
The relation between (h’ —hy)/H and f/fo which 
has been calculated for the Chapman distri- 
bution is given in the graph of Fig. 3 and in 
Table I. 

Appleton and Beynon‘ have noted that the 
virtual heights of reflection from a parabolic dis- 


*E. V. Appleton and W. J. G. Beynon, Proc. Phys. Soc. 
52, 518 (1940). 
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tribution fall in a straight line when plotted 
against 
f fotf 
fe fe-S 

The curve of Fig. 3 gives a line which is approxi- 
mately straight, in the region where f/f» is 
greater than about 0.4, when plotted against this 
function of f and fo. Because of the retardation 
in the lower part of the Chapman layer the 
straight part of the curve lies above the line for 
the equivalent parabola and has a slightly 
smaller slope, as shown in Fig. 4. This means that 
if the analysis for a parabolic distribution is 
applied to the virtual height-frequency record 
obtained from an ionospheric layer which actu- 
ally has the Chapman form, it will result in 
values for the true height of maximum and the 
half-thickness which are respectively too large 
and too small. Thus if the virtual height- 
frequency curve for a Chapman layer with 
hu = 300 km and H=50 km is constructed from 
Eq. (4) and Fig. 3, and then the Booker and 
Seaton form of the parabolic treatment applied 
to the curve, the resulting values of Ay and r 
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will be 307 km and 84 km, respectively. The 
parabola having these constants is shown by the 
dotted curve of Fig. 1. 

In practice the Booker and Seaton! method for 
the parabolic distribution may be converted into 
an analysis for the Chapman region by these rules: 

1. Determine hy, and 7, by the method given 
in reference 2. (Note: Subscripts have been used 
to differentiate between the values derived by 
the ‘parabolic’ and ‘‘Chapman”’ treatments.) 

2. Apply the following corrections: 


a. H=0.60r>, 
b. hue=hu,—0.14H, 
c. Te= 1.1075. 


TABLE II. Values of R/H as a function of f/fo. 
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, It may be-desirable, and it is equally easy, to 
analyze records directly in terms of the best fit 
to a Chapman distribution. The true height of 
maximum density, for example, may be found 
immediately as it is equal to the virtual height 
at the frequency for which f/f) =0.811. A number 
of other ratios of f to fo for which there are simple 
relations between hy, h’, and H are listed in 
Table I. If the virtual heights are read from the 
record at the frequencies corresponding to the 
tabulated values of f/fo, it is easy to determine 
both hy and H. (Note that corrections to be 
described below must be applied in case of the 
existence of a retarding layer below the reflecting 
layer.) In general it will be found that the values 
of H determined from the larger values of f/fi 
may, for daytime records, be greater than the 
values found from the middle frequency range. 
This is an indication of a deviation from the 
Chapman distribution of densities, a deviation 
which is presumably caused either by imperfect 
superposition of a number of density curves, each 
of which may have the Chapman form, or by 
upward diffusion of free electrons. 

By an approximate computation exactly 
similar to that described above, it is possible to 
determine the retardation of a ray that passes 
completely through an ionized layer having a 
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Fig. 5. The “retarda- 
tion” of a ray which 

normally through 
a layer having a Chap- 
man distribution. 


ie) 
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Chapman distribution. These data are exhibited 
in Table II and in Fig. 5, which is drawn to a 
scale of f/fo which is linear in log (f/fo—1). The 
quantity R/H is the “retardation,” or the con- 
tribution to the total virtual height which is due 
to the decrease in group velocity in passage 
through the retarding layer, divided by the scale 
height for that layer. The data given assume that 
the retarded ray passes completely through the 
layer, but the figures are reduced by only a few 
percent if the ray passes no further than a point 


1.02 
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fe 


which is three or four times the scale height 
beyond the maximum. 

Combining the data from Figs. 3 and 5, we 
find that, for a ray whose frequency is greater 
than the critical frequency for the lower layer 
and smaller than the critical frequency for the 


upper layer, 
h' —hyu R 
H Jr Hig 


where the subscripts F and E refer to the con- 


h'’=hur +H, 
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Fic, 7. An example of one method for analysing experi- 
mental virtual height-frequency data for true height and 
scale height. 


stants for the upper and lower layers, respec- 
tively. 

Figure 6 exhibits the variation of the virtual 
height with frequency for an hypothetical iono- 
sphere having the following characteristics: 
hur =300 km; Hr=50 km; huz=110 km; 
Hg=10 km. Here the ratio of the critical fre- 
quencies for the E and F regions is 0.3. It will be 
noted that, unless the retarding layer overlaps 
the reflecting layer significantly, the virtual 
height for the upper layer is independent of the 
true height of the retarding layer. 

It is probable that the best simple method of 
analyzing data for the F: region is to replot them 
to a. standard scale of f/f» and to find by inspec- 
tion'the best fit between these data and one of a 
family of virtual height curves, calculated for 
various values of H, which have been drawn on a 
transparent template. The height of maximum 
may be found from the distance from the point 
where the curves intersect on the template 
(where f/fo=0.811 in the case of no E retarda- 
tion) to the base line of the experimental data. 
This procedure, of course, yields the true height 
and scale height of the Chapman layer that most 
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closely corresponds to the actual distribution. A 
diagram of a rudimentary template of this sort 
is shown in Fig. 7, in which the circles represent 
points on an experimental curve supposedly seen 
through the transparent template, and in which 
the base line of the experimental data is shown 
dotted. In this case (00*24" on October 2, 1940) 
the F layer may be represented by the quantities: 
hu =39245 km; H=65+5 km; and N,,, 
= (2.46+0.01)10° ions/cm’. 

At times when E-layer ionization is in eyj- 
dence, the analysis may be made by preparing a 
number of templates of this sort from the data in 
Tables I and II, each template being for a 
specific ratio of E-layer critical frequency to that 
of the F region. When this has been done, it is 
simply necessary to select the proper template 
for each experimental curve before reading off 
the values of hy and H. This procedure implies 
that H is known for the E layer. 

In the presence of F; layer ionization the 
analysis is more difficult because the F; and F, 
layers overlap to such an extent that the virtual 
height for the F; layer cannot be obtained cor- 
rectly by estimating F; retardation in the way 
in which it can be done for the E layer. It is still 
possible, of course, to trace out the total virtual 
height by the laminar method, but this process 
has not been carried to the point where a sum- 
mary of the results can be expressed easily. A 
rough method for subtracting the effect of the F, 
region so that the constants for the F2 region 
can be obtained is evident, however, from the 
construction of Figs. 8 and 9. In Fig. 8 the 
diagram at the top represents a variation of 
density of ionization with height obtained by 
adding three Chapman distributions, as follows: 


E layer F, layer F; layer 
hy =110 km hy = 185 km hu = 300 km 
H= 10km H= 10km H= 50km 


Nuox™ 0.2 No Nox ™ 0.25 No Nuex= Neo 


At the bottom of the figure is the corresponding 
virtual height curve which was calculated by the 
laminar method. The dotted curve is the F; layer 
part of the figure as it would appear if the F; 
ionization were subtracted from the total. 
Exactly the same construction is followed in 
Fig. 9, with the important exception that here 
the height of the F; maximum has been taken as 
198 km. Under these conditions no subsidiary 
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maximum appears on the height-density curve 
and consequently there is no critical frequency 
for the F: layer. 

In each case it is apparent that the dotted 
virtual-height curve, representing the F:; layer 
in the absence of the F; layer, passes nearly 
midway between the minimum virtual height for 
the F: layer and the minimum virtual height for 
the F; layer. This clue enables us to obtain 
approximate values of the Chapman distribution 
constants for the F: layer. If the experimental 
virtual-height curve resembles the solid curve 
of Fig. 8 or Fig. 9, it is easy to sketch upon it a 
single smooth curve that will not differ from the 
curve which would be obtained from the F; layer 
alone by more than a few kilometers of virtual 
height. This curve may be compared with the 
templates described above in order to determine 
hy and H for the F: layer, while these constants 
for the F; layer can be found by the treatment 
which has been described for two layers. 

It is not within the scope of this article to 
discuss the presence of an E» layer or of other 
ionization between the E and F regions. Further 
study may yield a sufficiently accurate method 
of analysis for the retardation introduced by such 
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Fic. 8. An assumed distribution of electron densities in 
the ionosphere, and the corresponding virtual height- 
frequency diagram. The dotted curve would be obtained if 
the F, ionization were subtracted from the total. 
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ionization so that we can form some conclusions 
about its distribution. 

It should be noticed that the assumption of no 
magnetic field leads, in the case of the ordinary 
ray which alone is considered here, to values of 
the virtual height that are always too large 
except for the case of transverse propagation. 
These errors, become smaller at the higher fre- 
quencies. 

It may also be worth while to note that if the 
assumption of Chapman form for the various 
layers is at all justified, the construction of Fig. 8 
indicates that the scale height for the F, layer 
must be small compared with that for F2. It is 
evident that no F; maximum can occur unless the 
negative gradient of the F; ionization above the 
maximum is steeper than the positive gradient 
of the bottom of the F, layer. Since the curves 
are quite symmetrical near the maximum and 
have the same form, their slopes will have the 
same magnitude only if the scale heights for the 
two layers are approximately in the same ratio 
as their maximum densities of ionization. This 
argument, together with the known variations of 
the F; critical frequency, indicates that the F; 
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Fic. 9. Another assumed distribution of ionization with 
the corresponding virtual height-frequency curve. In this 
case the F; ionization is at such a height that no density 
maximum exists for that layer. 
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layer is very similar to the E in scale height, in 
maximum density of ionization, in the extent and 
symmetry of its diurnal variation, and perhaps 
in its behavior over longer periods. In case these 
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ideas are correct, we must look for an explanation 
of the seasonal appearance and disappearance of 
the F; critical frequency in the distribution of 
ionization with height in the F; layer. 
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Formulae for the calculation of transition probabilities in an atom with two electrons are 
derived for transitions from discrete s to both continuous and discrete p-states. 


I, INTRODUCTION 


T is of interest, especially on account of its 
application to astrophysics, to investigate the 
influence which the mutual interaction of the 
two electrons in helium has on the matrix ele- 
ments of the coordinates. The various perturba- 
tions, the coulomb interaction, the exchange, 
and the “polarization” all contribute to the 
energy. As a consequence, the chief quantum 
number, m, is no longer integral. These pertur- 
bations are taken together as the ‘‘ Rydberg 
correction,” 6, so that »=[mn]+6. We assume 
here that the results of the energy calculation 
are known! and wish to calculate the coordinate 
matrix elements for variable n.? In order to find 
the matrix element for an atom, 4 is taken from 
one of the calculated energies and the value of 
n is then substituted in our formulae. 

The Rydberg correction is particularly large 
for s-states. We should therefore expect the 
greatest deviations from the hydrogen matrix 
elements for the transitions sp and these are 
accordingly treated here. 


1. Transitions from discrete s to continuous 
p-states 


The orbital eigenfunction Y, after splitting off 
from the spin function, must be either symmetric 


1E. A. Hylleraas, Zeits. f. Physik 66, 453 (1930), 65, 
209 (1930). H. A. Bethe, Handbuch der Physik [1], Vol. 


24, Pp 347. 
#1 have to thank E. A. Hylleraas for this suggestion. 
See also reference 5. 


(para) or antisymmetric (ortho) in the electrons 
v= 2-4 W (71) Va(re) +V (re) W.(r;) }. 
The inner electron is in the ground state 


Yi=e 272, a=h*/moe*. 


Z is the nuclear charge in the field of which the 
inner electron moves. For the outer electron 
moving in the reduced field z, we will assume 
that the eigenfunction Y, may be expressed as 
the product of a radial function and a spherical 
harmonic function, ¥.=y¥(r)- Yin(@, ®). If this 
is substituted in the expressions for the matrix 
elements of the coordinates, exactly the same 
selection rules and the same coefficients, de- 
pending only on / and m, as in the hydrogen 
matrix elements result.? Only the integrals de- 
pending on the radius are altered. If we denote 
with capital letters N, L, M the quantum 
numbers corresponding to the state with the 
higher azimuthal quantum number, L =/-+1, and 
substitute for r the dimensionless quantity 
p=2sr/a, we obtain the expression R,,;""-a/2zs, 
which Gordon denoted by C. 
For the (sp)-transitions 


RY *2= (II aw+JIwnIn) K [Inn] (Land +I,?) "3 (1, 1) 
where 


Iaw= f Ynonip*dp; Inw= f vni'p*dp; 
0 4 


*W. Gordon, Ann. d. Physik 2, 1031 (1929). 
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-f vno?p'dp; =f e~20l*9*dp =2(2/Z)*; 
0 0 


v= f e~ 20/2) yi p*dp; r= f e~ 20/2) ,op*dp 
0 


2. The approximate functions and their nor- 
malization 


The difficulty lies in discovering suitable ap- 
proximate functions for the discrete spectrum. 
In his earlier calculation of the helium matrix 
elements, E. A. Hylleraas used the functions 
which the variation method yielded from the 
energy calculation. Lately,® Hylleraas has used 
complex integrals in order to extend the eigen- 
functions for integral m, the Laguerre poly- 
nomials, to arbitrary m. This method shows 
itself to be especially suitable for the calculation 
of the transitions to the continuum. 

The differential equation for y, if we neglect 
the polarization* and assume that the inner elec- 
tron effectively reduces the nuclear charge to ze, 
is 


py" +2py’— (1(/+-1)—p+p*/4n*)y=0, (2, 1) 


which is satisfied by 


¥(p) = pin-t¥(2ni)-+ ferme 


X(1—s)rgeedg, (2, 2) 


if the path of integration is taken so that at both 
ends 


(1 —¢) i+1+n¢l+1—ne—p (1-2) /2n (2, 3) 
has the same value. 

For the discrete s-states, following Hylleraas, 
we take the path from {= — © round ¢=0 back 
to {=—«, whereby (2,3) vanishes so long as 
p>0. If we let ¢=s/(1+s) and integrate by 
parts 


‘E. A. Hylleraas, Zeits. f. Physik 106, 395 (1937). 
ass. A. Hylleraas, Archiv f. Math. og Natv. 48, No. 4 

* We take account of the polarization approximately as 
it is included in n=[n]+6. 
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(0, +) 


ee (lf) /2n 


Ya(0) =(2mip)- if 


X(L—s)rgr dg 
(2,4) 


(0, +) 
=n(2nip)-* e—P(1—2)/(1+8)2ng—n—Igg 
=! 


Yn does not fulfill the boundary conditions at 
p=0, but can be used as an approximate function 
for the outer electron for the calculation of the 
matrix elements, as the deviation from the exact 
solution in the region p=0 contributes very 
little to the matrix integral J,y. 

The normalizing integral is easily found. It is 


Inn = 2n*{1—['(m) —(1/m) + (1/2n*) ] 
X (sinww)*/x?}. (2, 5) 


The expansion of ©’ (mn) =d?/dn? logI'(m) in recip- 
rocal powers of nm is’ W’(m)=(1/m)—(1/2n*) 
+ (1/6n*) —(1/30n5)+---- so that the coefficient 
of'(sin mr)? is small compared with unity. 

The interaction integrals for large, real 
quantum numbers are proportional to #~* and 
vanish at the boundary of the continuum. The 
eigenfunctions of the continuum are hydrogen 
eigenfunctions but in the reduced field, z. We 
can take them from Sommerfeld,* remembering 
that his p is our p/n 


II(n+/) ee 

sg IR 
II(n—1—1) 

_ Unt) — ol 2ei)- rf weeny 
~ W(n—I—1) july 


X(x—4)—*-Hle-rz/ndy, (2, 6) 


Sommerfeld gives for the normalizing integral 
- T 

J R?r*dr ae ellie IC (n+1+-1) |—*(an*/z); 
0 


k,=3s/|n\a. 


te, 1938) 3 20" Funktionentafeln (B. G. Teubner, Leip- 
zig, 1 

8A. nad aie 0 Atombau und Spektrallinien II (Fried- 
rich Vieweg und Sohn, Braunschweig, 1938), p. 121, Eq. 
(28). The eigenfunction ¥ is also given directly from our 
Eq. (2.2) by integrating by parts (2/+1) times over the 
closed path about ¢=0 and ¢=1. 
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The last factor appears when we normalize for 
the quantum number interval dn=1. 

The principal quantum number of the con- 
tinuum is purely imaginary and we put N =i, 
y= (2moE)~*hz/a, where E denotes the energy of 
the free electron. 

Hence, for the p-state 


Inn = 2v(v?+1)(1 —e7?*”) (2, 7) 
and with x=+4 


vn = (N?—1)Np-*(2ni)“ (1+”)*- 


(—1, 0) 


X NB 1+29)/2Ndn, (2, 8) 


3. The matrix integral J, 
The integral J,» is first integrated over p 


(0+) 
Inn =nN?(N?—1) (221)? f (1—¢)"-1g-! 


x rf (140) ¥29-¥-(q— no)—'dndt. (3, 1) 
(—1, 0) 


The integral with respect to 7» is equal to the 
negative residue at the pole 


no= —$+N(2¢—1)/2n. 


The transformation 2/=1—1/y together with 
the abbreviation 











+y I+ky 
f=n log —N log 
y 
=n log —2yv arc tan(ny/v); 
l—y 
k=n/N=n/iv (3, 2) 
gives 


Inn = 25n5(1 —1/N?)e***(22i)—! 
(1+) 
x efy*(1 —y?)-“"(1—k*y?)“*dy. (3, 3) 
0 
Integrating by parts 
Ian = —8n?(1—k*)-*{2(1—k?)(sinnx)/x 


+G2+2nGitGo} (3, 4) 
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where 


Go= —(sinnx)/z, 
(1+) (3, 5) 
G;=e'"*(2xi)- ‘f yay j= 0, A 2. 


This form of the integrand suggests the intro. 
duction of f as a new variable.® 

In order to calculate G; we deform the path of 
integration so that it passes through the regular 
point y= © ; in fact, we put y=re—**’5 along the 
path from 0 to © and y=ret**/® along the path 
back to 0. Then 


Tre 


1 y=r exp (ir/3) de! 
G;=— lim Rel ie-'n f yi—dy }. (3, 6) 
0 dy 


We expand f in powers of y and split up the 
coefficient of the lowest power, f= — $n(1—k?)¢. 
For the integrals G; we require y™ and y~ as 
power series in ¢. 


We let 
1 xr 
y lett 8=—= = a,(e?**/3(1 +k?)/5)?o@r 1)/3 
f p=0 
and , 
y-2e2in/3 = — = > by (e***/8(1 +82) /5) 72? (3, 7) 
r? p=0 


where dy, 5, are coefficients to be determined 
later. Then 


1 2p—1 
Gi=- 2 a,tt(——) ca +#)/5)" 
T Pp 3 


x [2n(1 —k2)/3]¢-2””? sin(—n—4+42p/3)x, 


1 
2=-) bd u()catersy 
x [2n(1 —k2)/3]?-2”)!8 sin(—n — 3+2p/3)r. 


. (3, 8) 
4. The calculation of the coefficients a,, b, 


The most convenient relation for the calcu- 
lation of a5, 5, is 


d 
3—y = (1—1/y*)(k?—1/y*) (4, 1) 
dg 


® Method of steepest descent. 
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TABLE I. Table of coefficients from the recursion formula (4, 4). 






















































, 5) 1- +41: +(a:)*- +(a2)*- +(a:)*- +(as)*- 
> 0 i 
| O.st111111  —0.22222222 
a =| —0.11111111 —0.14141414 —0.54545455 
~ eel 0 —0.35431235 —0.24941725 
ro. a, = | —0.01975309 0.13029440 0.00676854 0.16969697 
a =|  0.03209877 0.08669704 0.60214615 0.44129988 
a @ 0 0.15052357 0.21678623 0.48512744 0.09073510 
1 of - - 0.00463779 —0.06659891 —0.02767115 —0.35543989 —0.17129693 
| ay= | —0.01039911 —0.04400100 — 0.44506038 —0.70035860 — 0.80378893 —0.14795090 
lar n= 0 —0.06435739 —0.14151900 — 0.63408920 —0.54310902 —0.27641329 
" 1- +a: +. (a2)?- +(a2)* ° 7 (a2)* ° + (a2) ° 
bs = 1 2 
bh = 0.22222222 1.55555556 
coi © —0.72727273  —0.09090909 
by = | —0.22222222 —0.76923077 — 2.03418803 
, 6) bh =| —0.02716049 —0.71964085 —0.71874299  § —0.75151515 
by = 0 0.45194022 0.45084327 0.96558344 
by = 0.07654321 0.38762433 2.14170974 2.13151380 0.81838480 
the by = 0.00488601 0.34851639 0.65085535 1.95543443 0.91814678 
bio= 0 —0.23956751 —0.44040667 — 1.85259927 — 1.16922010 —0.29955557 
)¢. bu= | —0.02793128 —0.18833842 — 1.60090500 — 3.14552213 —4.50545115 — 1.79771769 
as 
i ther with | = 
which, toge (2p+3)ap=5bp_1— Dd (debp-g+2agdp-<); 
q=l 
, . (4, 2) ne (4, 4) 
= a " -- 
. q=0 _ bp=2apt+ 2 aap>. 
q=l 
7) give The next coefficients are given as polynomials in 
ao=1; ai=1; ag=(1—4k*)(4—k*)(1+k)2/7; 2 in Table I. The series (3, 8) for’G, and Gy are 
bo=1; bi=2; be=1+2a2 (4, 3) asymptotic series in which the terms decrease 
ed rapidly after the first. Only for very small », i.e., 
and the following coefficients from the recursion for high energies of the free electron, is a greater 
formulae number of terms necessary. 
5. The exchange integrals Jy and I, 
r, The integral Jy can easily be calculated, as the path of integration for y in (2, 8) is closed. We 
integrate for p, put 27=y—1 and m=2/Z and integrate by parts. 
Iy = —2°N*(N?—1)(2xt)-! (y+1)*¥-*(y—1)-*-4(2y +N) (y —yo)~'dy. (5, 1) 
(-1, +1) 
) The path of integration can be contracted round the pole yo = —N/m=—iv/m. Then 
Iy =44N4(N?—1)(N—m)*—*(N+m)-*-*(m —2)m5 (5, 2) 
or 
Iy = 4° y4(v? +1) (v?-+m?)-*(m — 2)m exp[ —2y» arc tan(m/v) }. (5, 3) 


a 
The integral J, is integrated first over p and then by parts 


I,=4n*m*(n?—m*)-!; [2n(1—m)(n+m)—!—1](sinnx)/r 





(0+) 
+4mn*(1—m)(n-+-m)-*(2ni)—* f s-"(1+st)—'ds}; t=(n—m)/(n+m). (5,4) 
=f 
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Expanding the denominator and integrating round a unit circle 





4m*n * 
I,=4n?m*(n+m)- ote-ieet- +4mn*(1—m)(n-+m)~ 2) t?-1/(b+1—n) t(sinnx)/x. (5,5) 
—n p=1 


This series converges rapidly only when ¢ is small compared with 1. As m=2/Z <1, generally n>m, 
In order to find an expansion in powers of n~! we transform the pole to 0 by putting st=y—1 ang 
deform the path of integration so that it runs from y=1—¢ in the negative imaginary half-plane to 
y=-— o and back in the positive imaginary half to 1—#. Then 


(0+) 1—t 
f s~"(1+st)—'ds = —1-4(cosne gf (1 —y)~"dy/y+2i sinnr [ (1 —»)-*dy/y). 
1 


—o 


The first integral is given by contracting both paths about the pole and has the value 27i. Syb. 
stituting y=x/(m-+m) in the second integral 


f (1—y)-"dy/y= f e*/x exp —x—n log(1—x/(n+m)) |dx. 


Expanding the second exponential function in powers of 1/m and integrating term by term, we have 
finally 


I,= —4m*n?(n? —m?*)—(sinnx)/x — 16m*n*(1 —m) (n? —m?)—*¢" 
X {cosaa + [(Et)w(2m) +e"C ](sinnr)/1r}, 
= —$m+(2m—1)/12n?+ (20m —4m*+ 6m? — 6m + 3)/90n*+ ---. (5, 6) 


In the case of helium, Z =2, m=2/Z =}, so that 
I,= —4n?(4n? —1)—(sinna) / — 16n*(4n? —1)-*t" 
X {cosna+ [(Ei)w(1) —e+ (e/40n') +--+ ](sinnr)/x}, (5,7) 


to an accuracy of over 0.01 percent for all m>1.5. 


Il. TRANSITIONS FROM DISCRETE s- TO DISCRETE /-STATES 


We wish to investigate the junction between the continuous and the discrete spectrum. We shall 
then be able to neglect the Rydberg correction for the discrete p-state in comparison with its much 
greater value for the s-state, so that N is a whole number. With this assumption E. A. Hylleraas! 
has calculated the matrix integrals. His series are, however, differently constructed so that an im- 
mediate comparison is not easily made. We shall therefore indicate shortly how our formulae for 
the discrete spectrum are to be altered. 

As approximate function for the discrete p-state we use (2, 2) in which, as N is integral, the path 
of integration can be taken round ¢{=0 and ¢=1. Integrating by parts three times gives (2, 8). 
The normalizing integral is 


Inn =2N(N?—1). (6, 1) 


At the boundary of the continuum No, vo, Eggs. (6, 1) and (2, 7) become identical. 
Equations (5,1) and (5, 2) for Iw are always valid because N>1>m. Our expression for the 
matrix integral J,y can be taken over unchanged so long as Nn. 
The matrix elements of the continuum therefore join continuously with those for the discrete 
spectrum. 
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When N <n the path of integration is taken along y =re***‘/*, This gives the same formulae (3, 8) 
for G;; instead of (1 — k*) (*?)/8.sinf —n—(j—2p)/3 ]e we now have 
(k? — 1) (-*)/8.sinf —n—(2j7+2p)/3]x. 


The asymptotic series (3, 8) for G; and G; become useless as N approaches n. We shall therefore, 
for the sake of completeness, also give expressions which can be used when N <n. 

We substitute y=(1+s)/(1—s) in Eq. (3,5) and expand the integral in powers of 
e=(N—n)/(N+n). Then 


G2t+Go=2ne(A+B)(sinnxr)/x; Gi=ne(A —B)(sinnr)/x 


ant Iw = 16n?N(1+) {(1+2)?+n(A +B) +n*(A —B) } (sinnr)/(n—N)x, 
where 
N o 
=— (N+n)"'So! +2 Sp-i°e?!/(N—n—p) +X So?t!(—€)?*!/(N—n+)), 
B=+(N+n)- SHE ‘Sputet!/(N— n— p)+E So?-"(—€)?-1/(N—n+)). 
For 


N- S(O re 
a+q ' 


we have the following recursion formulae 
(N+5) - (So?! + e Soot! ] = So? [(1+62)b+2Ne*], 
(N—a) -[ Sa? +? Sa41°]=S.°-[(1+¢*)a—2Ne*], 
So'—Si°=2S0°, Swe°=Sw°=0; Sy-= 


and 


In order to calculate transitions from discrete p- to continuous s-states the approximate function for 
the p-state must be changed, as also for transitions with higher /. The matrix integral J, can be 
calculated as the corresponding integrals G; can be reduced to our G,; and Gz by integration by parts. 
But the approximate function can then no longer be normalized. 

Investigation and also numerical calculations along these lines are in progress. 

I wish to thank Professor E. A. Hylleraas at Oslo, Norway, for his most valuable encouragement 
in this work. I should also like to thank him and Professor I. Waller at Uppsala, Sweden, for enabling 
me to carry it out during a critical time. | 
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The Mechanism of the Negative Point Corona at Atmospheric Pressure in 
Relation to the First Townsend Coefficient 


LEONARD B. LoEB 
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(Received February 11, 1947) 


On the basis of recent work by P. L. Morton and G. W. Johnson on ionization by electrons in 
divergent electrical fields, it is shown that at atmospheric pressure it is permissible to compute 
the ionization in electron avalanches by means of the Townsend expression e 
points of the order of 1X 10-* cm radius. The phenomena resembling a glow discharge, observed 
at atmospheric pressure with the negative points normally used, must be ascribed to ionization 
of the Morton type which can only be produced in the positive space charge fields built up at the 
point Jater in the course of pulse formation. The ionization at both positive and negative 
points must initially be of the Townsend type. This conclusion clarifies difficulties, and is 
consistent with the pulse build-up time of 0.4 microsecond just observed by W. N. English 
with a Synchroscope. An explanation is also given for the equality of starting potentials of 
pulses from positive and negative points in air, while in pure Nz and Hz, the negative pulses 
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HE work of Trichel,' Weissler,? and of Loeb, 

Kip, and Hudson® on the negative point 
corona in air and other gases at atmospheric 
pressure has established a close analog between 
its character and structure and that of the well- 
known glow discharge. This behavior, in contrast 
to the phenomena observed with a positive 
point corona under similar conditions, leads to in- 
consistencies in the earlier interpretations given.’ 
The difficulty was introduced through the use 
of the function e/* in calculating electron 
multiplication in the divergent fields about gaps, 
with @ representing Townsend’s first coefficient 
for ionization by electron impact, and x, the 
distance traversed in the field. The -function 
yields a maximum of ionization adjacent to the 
cathode surface.! On the publication of Morton’s* 
important study concerning the ionization by 
electrons in highly divergent fields, the difficulty 
became more acute. Morton showed that the 
Townsend function could not be applied to 
fields whose divergence was such that the field 
changed by more than 2 percent over an electron 
free path. Not only did Morton show that the 
ionization in such cases could reach values of an 
order of magnitude more than given by e/ =, 
but he further showed that under these condi- 


1G. W. Trichel, Phys. Rev. 54, 1078 (1938). 

2G. L. Weissler, Phys. Rev. 63, 96 (1943). 

3L. B. Loeb, A. F. Kip and G. G. Hudson, Phys. Rev. 
60, 714 (1941). 

*P. L. Morton, Phys. Rev. 70, 358 (1946). 
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tions the ionization near the negative point 
electrode was very small, becoming appreciable 
only at some distance from the point. Such a 
distribution of ionization is just the inverse of 
that given by e442, Thus the Morton process 
does not yield the conventional Townsend elec- 
tron avalanches. Morton’s findings are in ex 
cellent accord with the situation required to 
account for the observed mechanism near the 
cathode of a glow-like discharge. However, they 
introduce some difficulty when it is noted that 
the use of e/@* nicely accounts for the streamer 
and other phenomena at a positive point of 
moderate diameter,’ while the Morton conditions 
seem to apply to the negative corona from the 
same point, at fields of the same magnitude. 

A comprehensive extension of Morton’s initial 
studies over a pressure range from 0.01 mm to 
760 mm pressure in pure Hz and mercury free 
air has just been completed by G. W. Johnson‘ 
in the writer’s laboratory, and will be published 
shortly. On the basis of the new data the writer 
was able to make some calculations which enable 
the difficulty properly to be clarified. It appears 
that at atmospheric pressure for all gases studied, 
the geometrical field conditions for concentric 
cylinders are such that the electrons remain 
pretty well in equilibrium in the electrical fields 
down to wire diameters of the order of 1X10" 
























5 K, E. Fitzsimmons, Phys. Rev. 61, 175 (1942). 
*G. W. Johnson, Phys. Rev. 71, 278(A) (1947). 
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cm. With other geometrical configurations, such 
as hemispherical capped cylinders or confocal 

baloids,’? the value of the critical point 
radius might be slightly larger, but not materially 
so. Thus, the use of eS 42 in calculation for most 
corona points commonly used at atmospheric 
pressure is justified, unless space charge conditions 
alter the fields materially. 

This also justifies previous practice in com- 
puting the creation of the avalanche leading to 
streamer formation in positive point corona, 
which determines onset.’ In such calculations, 
theory and experiment have always been in 
satisfactory agreement.® ® 

In the case of the negative corona the calcula- 
tions for electron multiplication leading to evalu- 
ation of the onset threshold must extend beyond 
the first ionizing avalanches to which eS edz can 
apply; for, unlike the streamer, or the burst 
pulse, the negative corona pulse builds up more 
gradually.* In this instance, at the threshold, 
the first electrons leaving the negative point 
region move outward into lower fields producing 
e/edz electrons and positive ions as the static 
field requires. Photons produced by this ava- 
lanche, together with some positive ion impact 
on the negative point, then produce more ava- 
lanches, and a continuous or chain succession of 
such events builds up. During this period, how- 
ever, the relatively slowly moving positive ions 
are approaching the point, building up a notable 
positive ion space charge at the point which is 
only partially neutralized by the rapidly moving 
secondary electrons and their progeny leaving 
the point. In consequence, there is built up a 
very strong potential gradient at the cathode surface 
which is many times the geometrically fixed gradient 
starting the phenomenon. There is nothing in the 
corona environment to prevent the development 
of cathode gradients from becoming quite com- 
parable to those observed in glow discharges. 
When such gradients are achieved, in conformity 
with the Morton findings, the electrons then 


“1saa) H. Fisher and G. L. Weissler, Phys. Rev. 66, 102 
_* Recent studies of the Trichel negative corona pulses in 
air by means of a Synchroscope with high sweep rate by 
W. N. English in this laborat indicate that the un- 
resolved pulses of Trichel do indeed build up gradually, 
rising to a peak and declining over a time interval of 0.4 
microsecond, while streamers build up more rapidly. 
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TaBLE I. Data on the potential in volts for the beginning 
of bursts for positive and negative points (Vg+) and for 
the onset phenomena of Weissler (Vonser). 











Self-sustaining corona 
Gas Ver+ Vge- Vonest + ‘onset — 
He: 3500 2800 _— _- 
N: 4800 3700 — = 
Air hee reat _ = 5500 8000 
Air (English) $100+50 5000450 5550+50 8000-10,000 
Indefinite oscillo- 
graphically, indi- 
vidual pulses still 
resolved 








gain some twenty to a hundred electron volts of 
energy in ten to twenty electron free paths in 
the dark space. With ionization and excitation 
probabilities of the order of a few percent these 
electrons lose little energy. They thus do not 
ionize very much until they reach the weak 
field region beyond positive space charge. This 
mechanism yields the characteristic Crookes dark 
space and negative glow so often described.'** 
Using e/@4*, without the building up of a space 
charge and the Morton condition, as in the past, 
this observed condition was hard to understand.! 
Such a discharge will continue indefinitely with 
some fluctuations in magnitude until choked off 
by negative ion formation and space charge,’ or 
until alteration of the cathode surface by positive 
ion bombardment reduces secondary emission 
below the operating threshold.? The mechanism 
indicated at once removes the conflict between 
the positive corona, operating with e/@#, and 
the negative corona operating as a glow dis- 
charge, in that, initially, e/*¢* applies to both 
points but alters to Morton's condition as the 
discharge develops. The clarification, however, 
does more in that it makes possible an under- 
standing of sputtering** observed in negative 
coronas. With the positive space charge gradients 
the positive ions can strike the cathode with ten 
or more volts of energy instead of the estimated 
1 volt following on the Townsend calculation.* 
It had been observed by Weissler that in very 
pure Hz and N; the corona starting potentials 
were considerably lower for the negative than 
for the positive point. On the other hand, it has 
been generally asserted, without data being pub- 
lished, that the starting potentials for positive 
and negative points in air was about the same 
and higher than in pure Hey and Nz. This is 
illustrated by the data of Table I. These data 
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represent the beginning of bursts for the positive 
points, Vg+, and Trichel pulses, Vg—, for the 
negative points as shown by oscillograph for air, 
together with the corresponding data for the 
onset phenomena observed by Weissler in N:2 
and Hz». Weissler failed to specify his point 
diameter for the data given. He observed the 
onset visually and by current change for the 
self-sustaining coronas in He, Noe, and air. In this 
discussion the potentials for the beginning of 
pulses, Vg only is of interest. W. N. English 
kindly has measured the potentials for room air 
by using a 0.5 mm point. This, it appears, was 
also the point that Weissler used. The values in 
volts are: 

An explanation of this marked difference in 
behavior was indicated by Weissler,? but needs 
further clarification in the light of what has been 
stated above. There appear to be two factors 
present to effect such a change. First, in pure 
gases the secondary mechanisms responsible for 
onset of the negative corona are determined by 
secondary electron liberation by photoelectric 
effect from the cathode, by positive ion bombard- 
ment of the cathode, or by photoelectric ioniza- 
tion in the gas near the cathode. In the positive 
corona the fundamental secondary mechanism 
active is only photoelectric ionization in the gas 
leading to streamer, or burst pulse, production. 
Secondly, in air, electrons of low energy do not 
remain free but form negative ions of heat, of 
formation in the order of 0.2 volt. In pure H, 
and N; electrons are free. Thus, from the view- 
point of initiating electrons, the pure gases 
merely require adequate fields to initiate second- 
ary mechanisms near the point. The secondary 
liberation of electrons from a negative point in 
pure gases requires a relatively low field, since 
the photoelectric effect from the cathode metal is 
adequate in the relatively transparent gases to 
furnish such electrons at the start. With the 


*L. B. Loeb, Phys. Rev. 48, 684 (1935). See also N. E. 
Bradbury and G. Bloch, Phys. Rev. 48, 689 (1935). 
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discharge well developed, secondary electron 
liberation by positive ion impact on the cathode 
is also very effective. In the pure gases the photo. 
electric ionization in the gas is very meager, and 
it is not surprising to see that higher fields may 
be required for a discharge operating by this 
action only, in the case of the positive point jp 
pure H, and Nz. In air, on the other hand, the 
heterogeneous character of the gas facilitates 
photoelectric ionization in the gas, so that the 
positive threshold could be lowered relative to 
that in pure Ne and Hg, provided that Town. 
send’s first coefficient is not materially altered 
by Oz. On the other hand this circumstance lowers 
the efficiency of photoelectric ionization at the 
cathode by gas absorption of the photons. It can 
also raise the work function of the cathode be- 
cause of the action of oxygen on the surface, 
It is, however, hardly likely that at any time the 
photoelectric action in the gas can be the domi- 
nating element in negative point corona, though 
it, together with photoelectric action at the 
cathode, may help start the discharge in air, 
After space charges build up, positive ion bom- 
bardment doubtless predominates. Thus, in re- 
gard to secondary mechanism, the negative 
corona should usually have a lower threshold. 
However, neither of these discharges in air can 
start until the fields near the electrodes are 
sufficiently high so that a field strength to 
pressure ratio, X/p, exceeding 90, the value 
required to liberate electrons from the negative 
O, ions,® extends far enough out from the point 
to insure adequate electron ionization in the 
avalanches. This situation appears to be achieved 
in air at the starting potential for preonset 
streamers and burst pulses expected in positive 
corona.® Thus, it is this circumstance that prob- 
ably fixes the value for the onset of both the 
negative and positive corona pulses at about the 
same potential, under conditions where, without 
this difficulty, the negative starting point could 
well be lower. 
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A stack of N identical, parallel, and equidistant layers is used as a simple model of a layer 
structure. Irregularities in the relative displacements 5 parallel to the plane of the layers cor- 


respond to stacking disorder. The stacking disorder is described by means of a set of functions 
Wx (5) which measure the probability of finding a relative displacement 5 for layers M spacings 
apart. It is shown that x-ray diffraction data permit a direct determination of the Fourier 
coefficients of the functions Wy. Thus these functions can be synthesized and the nature of the 


stacking disorder directly deduced. 








INTRODUCTION 


RREGULARITIES in the relative displace- 

ment of the layers parallel to their own planes 
occur in many crystals of layer structure type. 
This stacking disorder gives rise to characteristic 
features in the x-ray and electron diffraction 
patterns of such crystals. The purpose of this 
paper is to show that the specific nature of the 
stacking disorder can be directly deduced from 
the observed diffraction effects. 

In order to demonstrate clearly the general 
principles upon which the method is based it is 
convenient to use a simple model of a layer 
structure. The extension of the method to layer 
structures of more general type can be made 
without serious difficulty. 


A SIMPLE MODEL OF A LAYER STRUCTURE 


Our model is a sequence of N identical, parallel, 
and equidistant layers. The periods within the 
layers will be denoted by a; and a, while a; is 
the period in the third principal direction (not 
necessarily normal to the plane of the layers). 
The individual layers in the stack will be 
numbered consecutively by means of an index 
L;, the layer L3=0 being chosen arbitrarily. The 
symbol 6z3; will be used to represent the dis- 
placement, parallel to the plane of the layers, of 
layer L; relative to layer L3=0. 

The position vector of any atom in the crystal 
has the form 


rj+Lia:+ Loa,+ Leas +523. 


tr; is the position vector of the jth atom within a 
unit parallelogram of a layer and L, Le, Ls are 
a set of three integers. 
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The model just described corresponds to an 
ordered layer structure if 5z3 varies periodically 
as the index L; is allowed to run through the 
sequence of integral value. We shall assume that 
the 5-values are not periodic, in which case the 
model corresponds to a layer structure with 
stacking disorder. 

We shall assume that the distribution of 
3-values throughout the stack is independent of 
the crystal surface. In order to study the con- 
sequences of this assumption we imagine the 
stack of layers to be periodically repeated in 
space. Our assumption implies that any sequence 
of N consecutive layers bodily removed from the 
periodic stack corresponds to a crystal specimen 
representative of the given disorder. Clearly there 
are N such different but statistically equivalent 
specimens. In discussing the disorder we shall 
therefore have to think in terms of the assembly 
of N specimens rather than in terms of a single 
specimen. 

The stacking disorder for the assembly can 
be described by means of a series of probability 
functions Wy(5), where M=0, 1, 2, ---, N—1. 
A function Wy represents the probability of 
finding a relative displacement 5 between two 
layers M spacings apart. 

The displacement 8 can be expressed in the 
form § = y:@1:+yea2 where 0<y;<1. Clearly Wo=0 
for all values 1, ye except yi:=32.=0. Since the 
functions Wy are even and periodic functions we 
may set 


Wun, 2) => Zz W*™ nye cos2xr(Hyyi+ His), (1) 
A, He 


where the Fourier coefficients conversely are 
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given by 


1 


1 
Win f f Wu cos2r(Hyy1 + Heys) dy,dy2. 
0 0 


THE X-RAY DIFFRACTION PATTERN 


This section is devoted to a theoretical study 
of the x-ray diffraction patterns of structures 
with stacking disorder. Since each individual 
layer is assumed to be a fully ordered two-dimen- 
sional lattice we shall obviously find that the 
scattering is restricted to directions with ‘Miller 
indices” (H,Hoz) when Hj, and H; (but not neces- 
sarily z) are integers. We shall show that the 
intensity of scattering in a direction (Hi Az) can 
be expressed in terms of the Fourier coefficients 
W™ He. Secondly we shall show that the in- 
dividual Fourier coefficients can be deduced 
from the experimental data and that hence the 
probability functions Wy can be synthesized. 

The mean intensity of the scattered x-rays 
from the assembly of N crystals discussed in the 
preceding section becomes 


I(s) =NA(s) > coss- Ma;(coss-5y)w. (2) 


In this expression s = 2r(k—ko) where kp and k 
are the wave vectors of incidence and of scat- 
tering. A(s) represents the intensity of scattering 
from a single layer. 

It is convenient to set 


8 = 24(Hibi+H2b2) + 24(p1b1 + pob2+2bs;) 


where Hj; and HA are integers and b;, be, bs the 
vector set reciprocal to a, a2, as. In this way the 
functions J#;H2(pipe2) represent the intensity 
distribution in the reciprocal lattice. 

The function A #;H2(p:, p2 2) for given z has a 
maximum at ~:=f:=0. The‘half width of this 
maximum depends upon the lateral dimensions 
of the layers. Unless the crystal dimensions 
parallel to the layers are of colloidal or subcol- 
loidal order of magnitude one will not be able to 
observe the variation of A #2 as function of pi 
and 2. For this reason we shall deal not with the 
intensity J#,H2(pipez) but with the experi- 
mentally significant quantity P#,#2(z) defined by 


PiyH:2(2) = f f TiyHo(pi, po, 2)dpid pr. 
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The expression for P#,#2 becomes 
PyHe(2z) = hd VBayH2(z) > W™ ay H2 cos2e Mz. (3) 


The quantities W“a,H2 are the Fourier coef. 
ficients which were introduced in the Preceding 
section. Ip is the intensity of the incident beam 
and 6V is the volume of a single crystal specimen, 
which for the sake of convenience is assumed to 
be so small that absorption and extinction can be 
neglected. The function Ba,#2 can be written jp 
the form 


Bayh, = V-'T | Pain|? 


where V=a,-a@X@s, ~#1H2 the structure factor 
and T = (1+cos*20)e*/2m*c*R? the Thomson scat- 
tering factor. 

Experimentally one does not measure directly 
the power distribution P#;#2(z) in the reciprocal 
lattice, but rather the power Q#1#2(qa) in terms 
of a parameter a which is suitable for the par- 
ticular experimental procedure which is used, 
The reduction of the observations Q#14:(a) to 
the reciprocal lattice distribution P#,,2(z) is, 
however, a trivial problem which involves the 
evaluation of the appropriate Lorentz factor. 

Since W™oo.=1 for any M the function Po(sz) 
will have maxima for integral values of z=H,, 
and the half width of these maxima is deter- 
mined by the total height of the stack of layers. 
In other words the stacking disorder has no 
effect on the reflections (00Hs). 

If the nature of the stacking disorder, i.e., each 
function Wy, is known, the calculation of the 
functions P#,#2(z) by means of Eq. (3) is a 
simple matter, and the nature of the stacking 
disorder can be found by a trial and error pro- 
cedure. Equation (3) suggests, however, a direct 
method for the determination of the specific 
character of the disorder. 

The distributions P#,H2(z) can be regarded as 
functions which are known from experiments, 
however, on a relative rather than absolute 
scale. The functions Bz,H2(z) are known if the 
structure of each layer has been determined. 

Hence, the quantities 


DayHe(2) = PHyHe/ Tod VBuy te 


are known functions except for a numerical scale 
factor. By the Fourier reciprocity theorem we 
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have 


3 
W" nye = f Duyn2(2) cos24Medz. 
-i 


Accordingly the probability functions Wy which 
characterize the stacking disorder can be syn- 
thesized. The unknown numerical scale factor 
can be determined from the requirement that Wz 
is everywhere positive and J{\f/Wudyidy.=1. 

In most instances it is found that the functions 
Wx approach a limit W., as M increases. It is 
therefore convenient to set 


Wu= W.t+Tmu 


where accordingly I'y represents fluctuations at 
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small values of M so that the Fourier coefficients 
I'™ ,H2 go towards zero with increasing M. 
One sees readily that 


Lu W™ mH cos2aMz = W?HyH2> mw cos2eMz 
+> uw [nye cos24Msz. 


Thus the functions W,, is generally to be asso- 
ciated with the sharp maxima, the functions I’ y 
with the diffuse maxima of the functions Pa, AH. 


If the function W,, is a constant, implying zero 
correlation between layers far apart, W#,4,=0 
unless H,=H,=0. Apart from the reflections 
(00H3) all maxima of P#,#2 will then be diffuse. 

The writer has used the method outlined above 
to determine the stacking disorder in several 
crystals of the layer structure type. 
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Localized states (Tamm levels), having energies distrib- 
uted in the “forbidden” range between the filled band and 
the conduction band, may exist at the surface of a semi- 
conductor. A condition of no net charge on the surface 
atoms may correspond to a partial filling of these states. 
If the density of surface levels is sufficiently high, there 
will be an appreciable double layer at the free surface of a 
semi-conductor formed from a net charge from electrons 
in surface states and a space charge of opposite sign, similar 
to that at a rectifying junction, extending into the semi- 
conductor. This double layer tends to make the work 


INTRODUCTION 


HE generally accepted view! of the nature 

of the rectifying contact between a metal 

and a semi-conductor is illustrated in Fig. 1 
which applies specifically to an excess semi- 
conductor. Figure 1a shows an energy level 
diagram of the metal and semi-conductor in 
equilibrium, but with the contact separated. 
The Fermi level is the same in both the metal and 
1See, for example, N. F. Mott and R. W. Gurney, 
onic Processes in Ionic Crystals (Oxford University 
Press, London, 1940), Chap. V. The theory of rectification 


is due in large part to W. Schottky. The most important 
of his papers is Zeits. f. Physik 118, 539 (1942). 


(Received February 13, 1947) 





function independent of the height of the Fermi level in 
the interior (which in turn depends on impurity content). 
If contact is made with a metal, the difference in work 
function between metal and semi-conductor is compen- 
sated by surface states charge, rather than by a space 
charge as is ordinarily assumed, so that the space charge 
layer is independent of the metal. Rectification charac- 
teristics are then independent of the metal. These ideas 
are used to explain results of Meyerhof and others on 
the relation between contact potential differences and 
rectification. 





the semi-conductor. As shown, the work function 
of the metal, x;, is greater than the work function 
of the semi-conductor, x2°, so that there is a 
contact difference in potential, x:—x2°. It is 
assumed that when the metal and semi-con- 
ductor are nearly joined, the potential distribu- 
tion is as shown in Fig. 1b. A double layer is 
formed such as to give a potential drop, go, from 
the metal to the interior of the semi-conductor 
equal to the contact potential difference. 

This double layer is assumed to consist of a 
space charge region in the semi-conductor, ex- 
tending to a depth of the order of 10-* to 10-‘ cm, 
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SEmi- 
CONDUCTOR 


(a) (b) 


Fic. 1. Energy level diagram for metal and semi-con- 
ductor in electrical and thermal equilibrium. 4; and yp: 
represent the Fermi levels in metal and semi-conductor, 
respectively, and ¢; and ¢2 the inner potentials. {; and {2 
are the chemical potentials and x; and x," the work func- 
tions. Figure la shows a wide gap with a contact potential 
difference x;—x,°. In 1b, the gap is very small. There is a 
space charge region in the semi-conductor which gives 
an electrostatic tential energy rise at the surface 
¢o0™x1—x:2". This is the usual picture in which no surface 
states on the semi-conductor are assumed. 


SEMI-CONDUCTOR 


and an induced charge on the metal surface. 
The space charge gives a rise in electrostatic 
potential energy at the surface of the semi- 
conductor. Electrons are depleted from the space 
charge region, giving a layer of high resistivity. 
If a potential is applied to the junction, most of 
the drop occurs across this barrier layer. If the 
potential of the semi-conductor is negative with 
respect to the metal, the electron energy levels 
in the semi-conductor are raised, and electrons 
may flow more easily over the potential hill into 
the metal. This is the direction of easy flow. 
On the other hand, if the semi-conductor is 
positive, the levels are lowered, increasing the 
height of the hill, and making it more difficult 
for electrons to travel from the semi-conductor 
to the metal. This is the direction of high 
resistance. 

According to this view, the equilibrium height 
of the potential hill, go, and therefore also the 
degree of rectification, depend on the work func- 
tion of the metal. For an excess semi-conductor, 
the larger the work function of the metal, the 
larger is the potential rise, and the larger is the 
reverse resistance of the barrier. If the work 
function of the metal is less than that of the 
semi-conductor there is a potential drop instead 
of a potential rise, and no rectification will occur. 
For a defect semi-conductor, in which the cur- 
rent is carried by holes, just the reverse is true: 
low metal work function gives high rectification. 

A number of investigations have been carried 


out which have verified these conclusions in some 
cases, and in other cases have not. H. Schwej 

as quoted by Schottky,’ has found a Correlation 
between the resistance of selenium rectifiers in 
the blocking direction and the work function of 
the metal. Selenium is a defect conductor, and 
high reverse resistance was found for low work 
function metals such as K, Na, Li, Ba, and low 
reverse resistance was found for such high work 
function metals as Ag, Au, Ni. The metal eles. 
trodes were put on by evaporation. 

W. H. Brattain,* working in this laboratory, 
has found a good correlation between degree of 
rectification and work function for metal con. 
tacts evaporated on cuprous oxide (a defect 
conductor) and on both N- and P-type silicon, 
Metals used, listed in order of decreasing degree 
of rectification, were Al, Ag, and Pt on cuprous 
oxide ; Pt, Be, Ag, Mg, and Al on N-type silicon: 
and Mg, Cd, Ag, and Pt on P-type silicon. He 
found that when contact is made to the semi- 
conductor by a metal junction in air, the rectif. 
cation is practically independent of the work 
function of the metal used. Results somewhat 
similar to those of Schweickert have been ob- 
tained by J. N. Shive* who studied the rectifica- 
tion characteristics of a number of contacts made 
by evaporation of various metals on selenium. 
Metal contacts, listed in order of decreasing 
degree of rectification, are Be, Zn, Pb, and Au. 

A. V. Joffe’ has studied contact potential 
differences and the resistances of the contacts 
formed from a large number of different semi- 
conductors and metals. Most gave very poor 
rectification characteristics. While there was 
some qualitative correlation of contact potential 
differences with contact resistance, quantitative 
agreement with theories of Schottky' and Davy- 
dov® was poor. 

2 W. Schottky, Physik Zeits. 41, 570 (1940). 

* Unpublished work done at the Bell Telephone Labora- 
tories in 1940 under the general direction of J. A. Becker, 
The author is indebted to Drs. Brattain and Shive for 
permission to quote their results. 

4 The designations N- and P-type refer to the direction 
of rectification. An N-type is an excess semi-conductor, 
and the direction of easy flow occurs when the sem 
conductor is negative relative to the metal. A P-type 
semi-conductor is a hole conductor, and the direction d 
easy flow is opposite. In both cases, the direction of easy 


flow is that in which the carrier moves from the semi- 


conductor to the metal. 
5A. V. Joffe, J. Phys. USSR 10, 49 (1946). 
* B. Davydov, J. Phys. USSR 4, 355 (1941). 











W. E. Meyerhof’ has recently been making an 
extensive study of the relation between contact 

tential difference and rectification for metal 
point contacts applied to Si and Ge. He has 
determined the potential rise, yo, from the 
rectification characteristics at different tempera- 
tures and has also measured the ordinary contact 
potential difference, or Volta potential, by a 
modified Kelvin bridge method. Preliminary 
results show little correlation between the two 
sets of measurements. In fact, he has found that 
go is practically independent of the work func- 
tion of the metal for metal point—Si rectifiers. 
The same metal may rectify with both P- and 
N-type silicon, in opposite directions. 

These negative results indicate that a closer 
analysis of the nature of the contact between a 
metal and semi-conductor is warranted. The 
main purpose of the present paper is to investi- 
gate the effect of electronic states on the surface 
of the semi-conductor on go. There has been 
considerable discussion in the literature con- 
cerning the possibility of surface states,* but 
little direct evidence as to their existence. 


IMPURITY LEVELS AND SURFACE STATES 


According to the modern theory of semi- 
conductors, there is an energy gap between the 
highest filled band of levels and the lowest state 
of the conduction band. In the language appro- 
priate to excess semi-conductors, conductivity 
results from electrons thermally excited to the 
conduction band from impurity levels. The im- 
purity levels, which are intermediate in energy 
between the filled band and conduction band, 
represent states in which the electron is localized 
around a foreign atom or other defect in the 
crystal lattice. In normal semi-conductors there 


™W. E. Meyerhof, “Contact potential difference in 
crystal rectifiers,” Technical Report No. 5, Univ. of 
Pennsylvania, BuShips contract NObs-34144, Aug. 10, 
1946; see Phys. Rev. 71, 727 (1947). 

*Surface states with energies in the “forbidden” band 
were studied first by I. Tamm, Physik. Zeits. Sowjetunion 
1, 733 (1932), and are often called Tamm states. Later 
theoretical work was done by R. H. Fowler, Proc. Roy. 
Soc. A141, 56 (1933), S. Rijanow, Zeits. f. Physik 89, 806 
(1934), A. W. Maue, Zeits. f. Physik 94, 717 (1935), E. T. 
Goodwin, Proc. Camb. Phil. Soc. 35, 205 (1939), W. G. 
Pollard, Phys. Rev. 56, 324 (1939), and W. Shockley, 
Phys. Rev. 56, 317 (1939). The last of these is the most 
critical, and shows under what conditions surface states 
are occupied in a normal crystal. 
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may be the order of 10-7 to 10-* impurity levels 
for each atom of the crystal. 

In addition to the impurity levels in the 
interior of the crystal, there may be localized 
states on the surface with energies in the 
“forbidden”’ region between the filled and con- 
duction bands. Shockley* and others have in- 
vestigated the conditions under which surface 
levels may be expected on an ideal crystal. His 
analysis, based primarily on a one-dimensional 
model, indicates that “in a plot of the energy 
spectrum versus interatomic distance the surface 
levels appear only at lattice constants so small 
that the boundary curves of the allowed energy 
bands have crossed.’ The number of such surface 
states is equal to the number of surface atoms, 
and, in a neutral crystal, the surface states are 
half filled. The conditions for surface states are 
fulfilled, for example, in a diamond-type lattice. 

Surface states may also result from surface 
imperfections, from foreign atoms on the surface, 
etc. On general grounds, there is good reason to 
suppose that the ratio of the number of surface 
levels to surface atoms may be much higher than 
the ratio of the number of impurity levels to 
atoms in the interior. 

The energy levels corresponding to the surface 
states may be discrete, or they may have a 
continuous distribution with all energies in the 
gap between the filled and conduction bands. 
There is little evidence of either an experimental 
or a theoretical nature on this point. Well- 
defined impurity states of the same nature, 
sparsely distributed over the surface, would all 
have about the same energy and form a discrete 
level. If the impurity atoms are densely dis- 
tributed over the surface, so that there is 
considerable interaction between them, a con- 
tinuous distribution in energy is to be éxpected. 
A continuous distribution is also expected for 
surface states on clean surfaces of the type 
discussed by Shockley, Goodwin, and others.* In 
the analysis to follow we assume that the surface 
levels are continuously distributed in energy. 
Many of the same general conclusions would 
follow, however, if discrete levels had been 
assumed. 

We will show that if there is a relatively high 
density of surface states, there may be a double 
layer on the free surface of a semi-conductor re- 
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sulting from a surface charge caused by electrons 
in surface states and a space charge of opposite 
sign extending into the crystal to a depth of 
10-* to 10-* cm. The double layer tends to make 
the work function independent of the impurities 
in the interior. If contact is made to a metal, 
the contact potential difference is compensated 
largely by a true surface charge rather than by 
space charge, so that the height of the barrier 
is largely independent of the metal. 


NATURE OF CONTACT 


Before discussing the effect of surface states, 
we will give a brief review of the nature of the 
contact between metal and semi-conductor. 

At a contact between any two electronic con- 
ductors there is an electric double layer which 
adjusts the potential of one relative to the other 
for equilibrium conditions. This potential differ- 
ence may be taken as a measure of the strength 
of the double layer. 

Some care must be taken in defining the 
absolute value of the potential. There is, of 
course, no way of measuring the potential inside 
a metal or semi-conductor. The potential just 
outside the surface can be measured, but the 
difference between the potential inside and out- 
side the surface depends on the double layer at 
the surface. Only differences in potential between 
different parts of a conductor can be observed 
directly. 

Nevertheless, it is desirable to consider the 
potential in the interior. The true time-average 
electrostatic potential is periodic, with the period 
of the lattice. The macroscopic potential, 6, may 
be defined as the space average of the actual 
potential over a region large compared with 
atomic dimensions. This definition is somewhat 
arbitrary ; another definition might yield a differ- 
ent value for the strength of the double layer at 
a contact or at a surface. However, only differ- 
ences in double layers have direct significance, 
and these are independent of the way the 
potential in the interior is defined. 

The general theory of equilibria’ shows that 
when two electronic conductors are in contact 

* An excellent discussion is given by R. H. Fowler and 
E. A. Guggenheim, Statistical Thermodynamics (Cambridge 
University Press, England, 1939), Chap. XI. The author is 


indebted to Dr. C. Herring for the method of presentation 
used above. 
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and in thermal equilibrium, the electrochemigy) 
potentials, 4; and ws, must be the same in both. 










#1 =e in equilibrium. (1) 





In the usual picture of an electronic condueto, 
# is just the Fermi level which determines the 
probability, p, that a state of energy & j, 
occupied. 

p=1/(1+exp[(E—x)/kT)). (2) 


The electrochemical potential in a region may 
be defined thermodynamically, in a way which 
is independent of any particular model, by the 
equation : 










re] 
u=—(U-TS), 
on 





where n is the number of electrons in the region, 
and U and S are the energy and entropy. The 
partial differentiation corresponds to a reversible 
change in which the volume and temperature ar 
held constant. The internal energy U includes 
an electrostatic potential energy, —e®, per elec. 
tron: The choice of the zero or reference level 
from which this potential energy is measured is 
arbitrary, and there is the same arbitrary choice 
in the definition of p. 

It is convenient to introduce a quantity which 
depends on the chemical constitution of the 
material, and on the electron density, but noton 
the electrostatic potential. This quantity, f, 
which is called the chemical potential, to dis 
tinguish it from the electrochemical potential, 
is defined by 



















f=u"u-¢, (3) 





where 





g=—eb 





is the electrostatic energy per electron. 

Let ¢, be the potential energy of an electron 
just outside the surface of a conductor. This 
potential is to be evaluated at a point where the 
image potential is negligible, but at a distance 
from the surface small compared with macro 
scopic dimensions. The work function of the 
surface is then: 
















Xs = Ge—H- (4) 





This is the energy required to take an electron 
from the conductor and place it at rest ata 
point, defined as above, just outside the surface. 
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Fic. 2. Diagram to illustrate double layers at a contact 
between two conductors in equilibrium. There is a double 
layer of strength ga—¢o at the surface of conductor 1, 


layer of strength ¢e2— ¢2, in the opposite direction 
om of csndeseat 2, and a double Ger of strength 
ga-ga= —(x1—x2) formed by surface charges on the 
two conductors. The total strength, y2— ¢u, is independent 
of the surface double layers. 


The work function depends on the double 
layer at the surface. By replacing u by {+¢, 
according to Eq. (3), the work function may be 


expressed in the form: 
Xs=Ge—o—f, (S) 


which shows explicitly that the work function is 
the sum of the energy required to take an elec- 
tron through the surface double layer, ¢.—¢, 
and a body term, —f, independent of the surface. 

It should be emphasized that the separation of 
the work function into a surface term and a body 
term depends on the precise way the potential 
energy, ¢, in the interior of the conductor is 
defined. We have suggested above that the 
electrostatic potential @ be defined as space 
average of the actual potential. Another defini- 
tion, might yield, for example, a smaller value 
for ¢.— ¢ and a correspondingly larger value for 
—{. Differences in double layers, such as might 
yield different work functions for different crystal 
faces of the same material, are not affected by 
the arbitrary definition of the inner potential. 
Since the double layer involves a potential 
difference, the value is independent of the 
reference level, or zero of potential, which is also 
arbitrary. 

If two electronic conductors are in contact, 
and in thermal equilibrium, the electrochemical 
potentials must be the same in both. Thus, 


fitgi=Seten, 


or 
g2— gi =Si-—$e. (6) 
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The strength of the double layer at the contact, 
¢1—¢2, is equal to the difference between the 
chemical potentials. Since the chemical potentials 
depend only on the internal constitutions of the 
two conductors, the double layer is independent 
of the work functions of the two surfaces before 
they are brought into contact. 

The situation when there is a small gap be- 
tween the two conductors in equilibrium is 
illustrated in Fig. 2. There is a double layer of 
strength ¢.:—¢ at the surface of conductor 1, 
a double layer of strength ¢.2— ¢e, in the opposite 
direction, at the surface of conductor 2, and a 
double layer of strength ¢2—¢31=—(xi—x:), 
formed by surface charges on the two con- 
ductors. The total strength of the double layer is: 


(Yer— $1) +(¢e2— Per) om (Ge2— $2) = ¢2— ¢1. (7) 


The surface double layers drop out, as they 
should according to Eq. (6). 

If the surfaces are in very intimate contact, 
so that the electron space charges of the two 
surfaces overlap, it is not possible to divide the 
total double layer into different parts which 
depend on the work functions of the individual 
surfaces. All that can be said is that the strength 
of the total double layer is g:—;. The work 
functions of the surfaces before they are brought 
into contact play no role.?® 


USUAL PICTURE OF RECTIFYING CONTACT 


If one of the conductors is a semi-conductor, 
and the other is a metal, it has generally been 
assumed that as they are brought into contact, 
no surface charge forms on the semi-conductor." 
Instead, there is a space charge near the surface 
which may extend to a depth of the order 10-* 
to 10-* cm. If the gap is small compared this 
distance, but is still large enough so that there is 


10H. Y. Fan, Phys. Rev. 61, 365 (1942); 62, 388 (1942), 
has attempted to calculate the electron space charge 
distribution in the double layer at the contact between two 
metals, and at the contact between a metal and a semi- 
conductor. Although his calculations are based on a highly 
idealized model, they may serve to give a rough picture of 
the charge density in the contact ion. He does not 
consider the effect of surface states on the semi-conductor. 

"N. F. Mott, Proc. Camb. Phil. Soc. 34, 568 (1938) 

ives a discussion of the contact between a metal and an 
insulator or a semi-conductor. A brief analysis is given of 
the effect of a surface charge resulting from electrons in 
surface states on the contact between a metal and an 
insulator. Mott assumes a discrete energy level for the 
surface states. 
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METAL 


SEMI-CONDUCTOR 


Fic. 3. Energy level diagram for metal semi-conductor 
contact illustrating notation used in text. The Fermi 
levels are wu; and yw: and the work functions x; and x2. 
The lowest state of the conduction band is denoted by C 
and the highest level of the filled band by F. If the surface 
states are filled to an energy ¢o below the conduction band 
there is no net charge on the surface atoms. 


no appreciable overlap, most of the energy drop 
x1—x2° will occur in the semi-conductor rather 
than in the gap. This situation is illustrated in 
Fig. 1b. The space charge raises the energy at 
the surface of the semi-conductor by an amount 
x1—xz2°. This is the usual picture which indicates 
that the height of the potential barrier should be 
equal to the difference in work functions, and is 
that which is discussed in the introduction. 

This picture is deficient in two respects: 

(a) The contact may be so intimate that a 
division of the double layer into one at the 
surface of the metal, one at the surface of the 
semi-conductor, and one caused by space charge 
may not be possible. There will be a double 
layer at the immediate interface, and another 
due to space charge, but it is not possible to say 
how the total strength is divided between them. 
A cruder way of stating this is that the work 
functions of the two surfaces are modified by the 
contact. However, if the contact is intimate, the 
separate work functions have no meaning. 

(b) There may be electronic states localized 
on the surface of the semi-conductor, so that the 
surface atoms can become charged. An appreci- 
able field may then exist in the gap. Part of the 
total drop from metal to semi-conductor will 
occur across the gap and part across the space 
charge region. If the density of surface states is 
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sufficiently high, part of the double-layer at the 
free surface of a semi-conductor may be formed 
by a surface states charge compensated by a 
space charge. The height of the potential barrie 
of a rectifier formed from the material will they 
be determined in part by the normal space charge 
of the free semi-conductor and in part by th 
work function of the metal. 









FREE SURFACE OF SEMI-CONDUCTOR 








We will first consider the free surface of a sem). 
conductor, and then discuss the rectifying cop. 
tact. The notation to be used is illustrated jp 
the energy level diagram of Fig. 3. The lowes 
state of the conduction band and the highest 
state of the filled band of the semi-conductor are 
indicated, with an energy gap €,. It is assumed 
that the distribution of surface states is such 
that the surface states give no net charge if 
the states are filled to an energy € below the 
conduction band. Since the Fermi level cuts 
the surface above the level determined by ¢, the 
surface as shown will be negatively charged, this 
charge resulting from electrons in states between 
€9 and the Fermi leve!. The picture applies tp 
an excess semi-conductor.'* In the body of the 






































o— 








_ Fic. 4. Schematic diagram showing how the potential 
rise ¢o is determined from the density of surface levels 
and the contact potential difference (see text). 






We refer to the net charge of the surface states a 
“surface states” charge to distinguish it from the 
charge, including space charge, in the surface layer, which 
we call “surface charge.” j 

18 The case of an excess semi-conductor seems to be easitt 
to visualize than that of a defect semi-conductor. Al 
results derived for one case, of course, apply to the other 
with obvious changes in signs of the charges. 
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semi-conductor, the Fermi level is an energy ¢ 
below the conduction band. The space charge 
region extends for an approximate distance / into 
the semi-conductor, giving a potential energy 
rise go at the surface. An energy ¢, is required 
to remove an electron from the lowest state of 
the conduction band near the surface to a point 
just outside the semi-conductor. The work func- 
tion, x2, depends on ¢, as well as on ¢o: 


x2=¢st got. (8) 


The amount and extent of the space charge 
inside the free surface is determined by the 
density of surface levels. For zero external field, 
the positive space charge raises the potential at 
the surface by an amount just sufficient to give 
a compensating negative surface charge. The 
larger the positive space charge region, the larger 
is go, and the smaller is the negative surface 
states charge. For some ¢o the two will be equal 
in magnitude. This is the equilibrium value. 

These relations are indicated in a schematic 
way in Fig. 4, which shows the variation of i, 
the total space charge per unit area, of o,, the 
surface states charge, and of the surface charge, 
o,+0;, With go. The other curves on the diagram 
will be referred to later in connection with the 
discussion of a metal contact. For a uniform 
space charge in the boundary layer, o; is pro- 
portional to the square root of go (see the 
appendix). The plot of o, has a steep slope corre- 
sponding to a fairly high density of surface states. 
It passes through zero (point £ in Fig. 4) when 


go=eo—f, (9) 


because, according to the definition of ¢ the 
surface states will then be filled up to the level 
corresponding to zero surface states charge. The 
point A represents a value of go for which 
¢,+0;=0, which corresponds to a neutral surface. 
It can be seen that if the density of surface states 
is high, this value of go will be close to that given 
by Eq. (9). 

For the limiting case of a very high density 
of surface states, the line ¢, becomes vertical, 
points A and E coincide, and Eq. (9) must be 
satisfied ; hence, according to Eq. (8): 


X2= Get €0. (10) 


The work function is determined entirely by the 
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Fic. 5. Energy level diagram showing contact potential 
difference between two ends of a semi-conductor which 
changes from N-type to ey along its length. (a) High 
density of surface states, small contact potential difference. 
(b) No surface states, large contact potential difference. 


surface, and is independent of the position of the 
Fermi level in the interior. 

On the other hand, if the density of surface 
states is small, ¢,+0; is approximately equal to o; 
and the condition for a neutral surface leads to 
a small value for go. In the limiting case of 
vanishing surface states charge, gy=0, and the 
work function is 


xx=o+f, (11) 


which of course does depend on the position of 
the Fermi level in the interior of the semi- 
conductor. 

Some semi-conductors can be made either 
excess or defect, or N- or P-type,‘ depending 
on the nature and -concentration of impurities. 
The distribution of impurities in a single sample 
may be such as to make one end N-type and the 
other end P-type. Silicon and germanium are 
examples of materials which may behave in this 
way. For an N-type conductor, ¢ is small; for a 
P-type conductor ¢ is almost equal to the energy 
gap ¢,, which is the order of one electron volt 
for these materials. 

If the density of surface states is low, a large 
difference in work function between N- and 













P-type is to be expected. The difference will be 
almost equal to the energy gap e,. For a high 
density of surface states, the difference in work 
function will be small. Figures 5a and b show 
schematic energy level diagrams of these two 
limiting cases. Both show the variation in poten- 
tial along, and the work functions at the two 
ends of a sample which changes from N- to P- 
type from left to right along its length. Figure 5a 
is that corresponding to a high density of surface 
states. There is a space charge region at each 
end. The potential at the left is raised so as to 
bring the energy level corresponding to a neutral 
surface close to the Fermi level. The potential 
at the right is similarly lowered. The difference 
in work functions, xp—xv, is small. Figure 5b is 
the case corresponding to a small density of 
surface states. There are no space charge regions 
at the ends, and the difference in work functions 
is large. 

A difference in work function can be detected 
experimentally as a difference in contact poten- 
tial. Meyerhof’ has in this way measured the 
differences in work functions of various samples 
of N- and P-type silicon and has investigated the 
effect of surface treatment on these differences. 
He finds that the difference between N- and 
P-type samples is about 0.25 ev, which is con- 
siderably smaller than the energy gap of about 
1.1 ev. It is quite possible that there is a high 
density of surface states, either on the pure 
material or resulting from surface films or im- 
purities, which is sufficient to account for this 
discrepancy. 

We have so far discussed the density of surface 
states in a purely qualitative way. It is of interest 
to make an estimate of the density required to 
produce an appreciable space charge layer at 
the free surface of a semi-conductor. The density 
depends on the density of charged centers, or 
donors, in the boundary layer. If there are N 
such centers per unit volume, and if the thickness 
of the boundary layer (assumed uniform) is /, 
we have for the total space charge per unit area :' 


o;=eNl. (12) 
The potential energy rise at the surface is 
go = (2n/x)e NP, (13) 


where «x is the dielectric constant. 
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Let the number of surface states in the er 
interval de be nde/éo. We assume, for simplicity 
that m is a constant, which is of the order of 
the total number of surface levels per unit area 
with energies in the gap. The change in surface 
states charge density corresponding to a poten. 
tial energy rise go at the surface is then 










Ao, =en¢o/ €o. (14) 






The condition for an appreciable space charge 
layer is that go be of the order of magnitude of & 
when Ag, =o;. This requires that be the order 
of Ni. Taking, for example, x~15, and N~19", 
Eq. (13) gives /~10-5, so that 


n~ NI~10" states/cm?. (15) 













This density corresponds to about one surface 
state per thousand surface atoms. Other things 
being fixed, the limiting density of surface states 
is proportional to the square root of the density 
of charged centers in the boundary layer. 








METAL SEMI-CONDUCTOR CONTACT 






As a metal surface approaches that of a semi- 
conductor, the electrostatic field in the gap 
increases. There is a surface charge oy on the 
metal surface, and a charge of equal magnitude 
and opposite sign divided between a surface 
charge, o,, and a space charge, o; (the latter has 
been defined as the total space charge per unit 
area of surface). Thus 


om = —(0.+0;). (16) 


In the case of the metal, the surface charge causes 
only a very slight change in work function. 
However, the charge on the semi-conductor may 
cause appreciable changes in its work function, 
so that its value x2 will differ appreciably from 
the zero field value x2°. These changes can be 
analyzed with the aid of Figs. 3 and 4. 

According to Fig. 3, the field in the gap 
between metal and semi-conductor is (x1— x2)/¢a. 
This is produced by the surface charge which 
gives a field 4r(¢,+0;). Equation (8) expresses 
X2 as a function of go. Equating the two ex 
pressions for the field leads to: 
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=04(¢0) +;(¢0). (17) 
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The left-hand side of the equation is plotted in 
Fig. 4 as the line “L,”’ and is a straight line of 
negative slope. The solution to (17) is repre- 
sented by the point B, which is the intersection 
of “L” with the line representing o,+¢;. The 
line “L’”’ crosses the axis at the point D, which 


corresponds to x2=X1, OF 
$o=X1— e—$. (18) 


It can be seen that if the density of surface 
levels is high, so that o, and o,+0; have steep 
slopes, the point B will lie close to the point A 
which gives the value of go for the free surface. 
In this case, 

gow eo—f, (19) 


and is practically independent of the work 
function of the metal. This is a possible explana- 
tion of the results of Meyerhof, who found that 
the value of go for metal points on silicon does 
not depend very much on the metal used. 

In the limiting case of vanishing surface states 
charge density, ¢,=0, and the solution of (17) 
is given by the point C of Fig. 4. The value of 
go is then close to that given by the intersection 
of the line L with the horizontal axis (point D). 
In this case, 

go~X1— Ge—$, (20) 


and the usual picture of the contact between 
metal and semi-conductor applies. 

As the metal is brought closer to the semi- 
conductor, the slope of the line “L’’ increases 
with decreasing ‘‘a’’ and one might expect a 
transition from a solution near A to one near D 
with a resulting value of go given by Eq. (20). 
However, even if there is only one surface state 
per hundred surface atoms, the slope of the line 
¢, is so great that the solution still lies near 
point A even when a is reduced to atomic dis- 
tances, say 3A. This may be verified by making 
a rough calculation based on Fig. 4. For the 
value of go to be largely determined by the semi- 
conductor, the value of gp for the solution B 
must lie closer to A then to D. The necessary 
and sufficient condition for this is that the slope 
of the line o, be much greater, in absolute value, 
than the slope of “‘L.”” We have 


slope of o,=m,=en/€o, 


slope of L= —m,= —1/4mea. (21) 





SURFACE STATES AND RECTIFICATION 









The condition that m,/my, be large is 


m,/m,=4nre’an/eo>1, (22) 
or 
n> €o/4re*a. (23) 


Setting for example, ¢)~ 10-" erg and a~3 X 10-8 
cm, this condition requires that 


n>>10", (24) 


Thus if there is appreciably more than about one 
surface state per hundred surface atoms, the 
metal work function will have little effect on go. 
This is about an order of magnitude larger than 
the surface density required for the existence of 
an appreciable boundary layer at the free surface." 

The theory is worked out in detail for the 
special case of a uniform Schottky exhaustion 
layer! in the appendix. 


CONCLUSIONS 


If the density of surface levels with energies 
in the “forbidden” band is sufficiently high 
(>~10"/cm?), there will be a double layer at 
the free surface of a semi-conductor formed from 
a surface states charge and a space charge of 
opposite sign. The space charge region is similar 
to that which exists at a rectifying contact. This 
double layer tends to make the work function 
independent of the height of the Fermi level in 
the interior, and so independent of the impurity 
content. 

The total strength of the double layer at a 
rectifying junction between a metal and semi- 
conductor is fixed by the difference in chemical 
potentials, and so depends on the body proper- 
ties of the metal and semi-conductor, and is 
independent of the work functions of the surfaces 
before they are brought into contact. The double 
layer consists of the following parts: 

4S. Benzer, Phys. Rev. 71, 141 (1947), has recently 
reported on the current-voltage characteristic observed 
when contact is made between two pieces of the same 
(homogeneous) germanium crystal. He states that “‘the 
characteristic observed for both polarities is more the 
order of the back resistance observed when either piece of 
the crystal is contacted with a metal; in both directions 
the negative resistance at high voltage appears, which is 
typical of the back characteristics of metal-Ge contacts 
using these alloys.” The presence of a space charge layer 
at the surface of each piece which is little modified by 
contact would result in a characteristic similar to that of 
two rectifiers in series opposition. Such boundary layers 
are to be expected if the density of surface states is the 
order of that given by Eq. (24) and if the contact is not 


too intimate. Benzer’s results are thus indirect evidence 
for surface states on germanium crystals. 
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(1) A double layer of atomic dimensions at 
the metal surface. 

(2) A double layer of atomic dimensions at 
the semi-conductor surface. ; 

(3) A double layer formed from surface 
charges on the metal and semi-conductor, both 
of atomic dimensions. 

(4) A double layer formed from a surface 
charge of atomic dimensions and a space charge 
extending to a depth of 10-* to 10-* cm into the 
semi-conductor. 

The strengths of the double layers may be 
estimated in different cases as follows: 

(a) If the density of surface levels is suffi- 
ciently high (> ~10'*/cm?) the double layer (4) 
will be independent of the metal, and will be 
the same as that for the free surface of the semi- 
conductor. The rectification properties will then 
be largely independent of the work function of 
the metal. The difference in contact potentials 
is compensated by (3). 

(b) If the density of surface levels is small 
(<~10'*/cm?), the double layer (3) will be 
small, and (4) will be determined approximately 
by the difference in work functions. 

(c) If the contact between the metal and semi- 
conductor is very intimate, it may not be possible 
to distinguish between the double layers (1), 
(2), and (3). The metal will tend to broaden 
the surface levels, but if the broadening is small 
compared with the energy gap, conclusion (a) 
will still be valid. 

(d) If the broadening of the surface levels by 
the metal is large, no conclusions about the 
space charge can be drawn from measurements 
of contact potential differences. 

It is believed that all of these cases may be 
realized. 

The author is indebted to various members of 
the technical staff of the Bell Telephone Labora- 
tories, particularly to W. Shockley, C. Herring, 
and W. H. Brattain, for the benefit of numerous 
discussions concerning the subject matter of this 
paper, and for helpful comments on the prepara- 
tion of the manuscript. 


APPENDIX 
Theory for Schottky Exhaustion Layer 


The detailed calculation of gy» requires a 
knowledge of the distribution of space charge in 


the semi-conductor and its dependence on O. 
The dependence of surface charge on ¢p is algo 
required. We will carry through the calculation 
explicitly only for the case of a Schottky exhauys. 
tion layer! with a uniform density of charge, 
This case is particularly simple and brings oyt 
the essential features of the problem. For sim. 
plicity we also assume that the surface states are 
uniformly distributed in energy. 

Let eN be the positive charge density in the 
barrier region, assumed constant. Let / be the 
thickness of the charged layer. Let nde/€, be the 
number of surface levels per unit area with 
energies in the range de. We assume that n is q 
constant, independent of energy. For other nota- 
tion see Fig. 3. 

We have the following relations: 

(a) The total space charge per unit area is: 


o, = Nel. (12) 


(b) The potential energy at the surface of the 


semi-conductor is 
2re2N 
¢o= #, (13) 
K 
where «x is the dielectric constant. 
(c) The surface states charge is 


go, = —en(€o— po—f)/€0. (15) 
(d) The total surface charge is 
oi+¢.= Nel —en(eo— go—$)/€0. (25) 


(e) The work function of the semi-conductor 

is: 
2re?=N 
xX2= Geth+ P. (26) 
kK 

Equating o:+o, with (x1—x2)/4mea, we get 

Eq. (17) from which go or the thickness of the 

barrier layer, /, can be determined. With the 


notation: 


n 
bb=—(1—$/¢0), (27a) 
N 


1= k€o/2me7n, (27b) 
Le? = «(x1 — @ —$)/20e*N, (27c) 

Eq. (17) reduces to: 
P/ly+1—lo = (le? —P?) /2«a, (28) 


which is a simple quadratic equation for /. 
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Whether the semi-conductor or the metal pre- 
dominates in determining /, and thus the barrier 
height, depends on the relative magnitudes of 
the terms on the left and right sides of Eq. (28). 
The terms on the right side will be negligible if 


1,<K2«a (29) 
and 
lo>1,?/2 «xa. (30) 
The first condition requires that 
n> €o/27e*a, (31) 


and the second that 
n>>(€0/41e*a)(x1— gs —$)/(€o—-$). (31a) 
These conditions are essentially equivalent to 
that given by Eq. (23), and lead to the conclusion 
that »>>10"* in order that the metal have little 
influence on the space charge region. In this 
case of high density of surface states, the equa- 
tion for the layer thickness reduces to that for 
the free surface of the semi-conductor : 


P/l44+1—1)=0. (32) 
Equation (32) may, of course, be used to 
estimate the thickness of the space change region 


at the free surface of the semi-conductor, regard- 
less of the density of surface states. The limiting 
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case of high density corresponds to 


Lo>l,, (33a) 
or 
n®*> Ne? /2re*(eg—f). (33b) 
This is essentially the requirement, stated above 
Eq. (15), that m be large compared with Ni. 
The approximate solution of Eq. (32) is then 


P lo). 


The product Jol; is independent of n, and is just 
the square of the thickness of the Schottky layer 
for a barrier height 

¢o= €o—f. (34) 
This is the condition that the Fermi level cross 
the surface near the energy corresponding to 
zero surface states charge. 

In the limiting case of a vanishingly small 
density of surface states, the thickness of the 
barrier layer is determined by setting the right- 
hand side of Eq. (28) to zero. This gives 


$o=Xi— ¢s—f. (35) 
The height of the barrier is equal to the difference 


in work functions. 
Equation (28) may be used for intermediate 


cases. 
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University of Pennsylvania, Philadelphia, Pennsylvania 


The rectifying portion of a crystal rectifier is the contact 
between a small point made of metal such as tungsten, and 
a semiconductor such as silicon or germanium containing 
suitable impurities. The potential energy of a conduction 
electron near the contact determines the rectifying action 
of the crystal rectifier. The most important feature of this 
potential energy, as far as the present paper is concerned, 
is the height of the potential barrier, which the electrons 
have to overcome when they pass from the metal to the 
semiconductor or vice versa. The height of the barrier is 
called here contact potential difference (c.p.d.), because 
theoretically it is equal to the difference in the work func- 


* Condensation of a dissertation in Physics presented to 
the Graduate School of the University of Pennsylvania in 

rtial fulfillment of the requirements for the degree of 

doctor of Philosophy. Dissertation published with limited 
distribution under ‘‘University of Pennsylvania, BuShips 
Contract NObs-34144, Technical Report No. 5, August 10, 
1946.” Copies available from University of Pennsylvania, 
Department of Physics, Philadelphia 4, Pennsylvania. 


(Received February 16, 1947) 


tions of the substances in contact. The c.p.d. has been 
measured using both n- and p-type silicon and different 
metallic contacts. (The c.p.d. can be obtained from the 
variation of the contact resistance with temperature.) The 
work function differences (w.f.d.) between the same sub- 
stances were obtained independently by a parallel plate 
condenser method (Kelvin method). The results showed no 
correlation between the c.p.d. and the w.f.d. The c.p.d. 
is practically independent of the kind of metal used and 
also of the structure of the silicon surface. These results are 
in contradiction to the present theoretical model of the 
silicon crystal rectifier. 


** This work was done under Contracts OE Msr-388 and 
NObs-34144 between the Trustees of the University of 
Pennsylvania and the Office of Scientific Research and 
Development and the Navy Department, Bureau of Ships, 
respectively, which assume no responsibility for the ac- 
curacy of the statements contained herein, 

*** Now at the University of Illinois. 
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Fic. 1. Energy level diagrams for metal and 
semiconductor in equilibrium. 


INTRODUCTION 
HE essential constituents of a crystal recti- 
fier are a metal and a semiconductor in 
contact. Under suitable conditions these contacts 
have rectifying properties. 

In general one distinguishes large-area-contact 
rectifiers, such as copper oxide and selenium 
rectifiers, from point-contact rectifiers, such as 
silicon and germanium rectifiers, which will be 
described here. In both kinds of rectifier the semi- 
conductor is in the form of a slab, making contact 
on both sides to a metal. Although each side will 
rectify in the opposite direction from the other, 
a net rectification is obtained in the case of the 
large-area-contact rectifiers by making the in- 
trinsic rectification property of the two contacts 
different. In the case of point-contact rectifiers 
the same effect is obtained by making the area 
of one contact very much smaller than the area 
of the other.! Any voltage drop across the rec- 

! Actually the ratio of the large to the small area must 


be appreciably larger than the back to front resistance 
ratio of the rectifier. 


MEYERHOF 


tifier is then concentrated across the small area 
(high resistance) contact; in considering the 
theory of the point-contact rectifier one May 
justifiably neglect the large area contact.! 


THEORY OF THE SEMICONDUCTOR. 
METAL CONTACT 


When a metal and a semiconductor are brought 
into contact, the electron distributions adjust 
themselves in such a way as to make the average 
thermodynamic potential? of the conduction elec. 
trons in the solids everywhere the same. This 
adjustment is always accompanied by a transfer 
of electrons from the solid with the lower work 
function® to that with the higher work function, 
The transferred electrons stay in the neighbor. 
hood of the contact surface and, together with 
the positive charges they leave behind, form a 
double layer which just compensates for the 
original difference in work functions. 

The potential (resulting from a superposition 
of the original lattice potential and of the poten- 
tial caused by the double layer) in which the 
conduction electrons move has been investigated 
carefully by Fan.‘ In the case of silicon and 
germanium, though, one may make certain ap- 
proximations resulting from the fact that the ac- 
tivation energy of the electrons (or holes) from 
the impurity levels into the conduction band is 
only of the order of a few hundredths of an 
electron volt® so that practically all the impurities 
are ionized at room temperature. The exact shape 
of the potential near the contact then depends 
uniquely on the distribution of impurities in the 
semiconductor.*-* For example, a uniform impur- 
ity distribution produces a parabolic potential.*’ 

In the present considerations, it is necessary 
only to note that the potential jump which a 
conduction electron experiences in traversing the 


2See R. H. Fowler, Statistical Mechanics (University 
Press, Cambridge, 1936), ny XI. ‘ 

* Work function is defined here as the difference in po- 
tential of an electron just outside the surface of a solid 
and of an electron with an average thermodynamic potential 
energy. 

4H. Y. Fan, Phys. Rev. 62, 388 (1942). 

5 K. Lark-Horovitz, private communication; G. L. Pear- 
son and W. Shockley Bull. Am. Phys. Soc. [6] 21, 9 (1946); 
B. Serin, University of Pennsylvania dissertation, Phys. 
Rev. (to be published). : : 

*N. F. Mott and R. W. Gurney, Electronic Processes m 
Ionic Crystals (Clarendon Press, Oxford, 1940). 

7H. A. Bethe, NDRC Div. 14 report, MIT Rad. Lab. 
43-12, Nov. 23, 1942 (Copies obtainable from Publications 
Board, Department of Commerce, Washington, D. C.). 

® B. Serin, Phys. Rev. 69, 357 (1946). 
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contact surface between a metal and a semi- 
conductor is, to a first approximation, equal to 
the difference in work functions of the two solids. 
If ¢u is the work function of the metal and ¢5 
the work function of the semiconductor, the work 
function difference, ¢, is defined here by 


g=os—om, (1) 


when the solids are in contact. When they are 
separated the same symbols with the subscript 
zero will be used. Figure 1 shows the various 
possible potential diagrams for these two situa- 
tions.* Since the activation energies of the im- 
purity states are only a few hundredths of an 
electron volt, as was mentioned above, the posi- 
tion of the average thermodynamic potential 
energy level may be taken with sufficient accu- 
racy at the top and bottom of the forbidden zone 
for n- and p-type semiconductors respectively.* ® 

The theory of the silicon or germanium point- 
contact rectifier was given by Bethe,’? who 
pointed out that in the case of these substances 
the number of conduction electrons colliding with 
the potential barrier is appreciably larger than 
the number of electrons scattered by the lattice 
and the impurities in the barrier region. Therefore 
the diffusion theory, developed by Schottky,’° 
Mott,® and others and generally applied to large- 
area-contact rectifiers, is not valid, but a simple 
diode theory will be more satisfactory for calcu- 
lating the current flow across the contact. If n is 
the density of impurities in the semiconductor, 
v is the thermal velocity of the electrons per- 
pendicular to the contact surface, e is the absolute 
value of the electronic charge, A is the effective 
contact area, k is Boltzmann’s constant, T is the 
absolute temperature and m is the effective mass 
of the electron, then the current i through the 
contact is 


i= joA -exp(—e]¢|/kT)-[expeV/(kT)—1], (2) 


where 
jo= nev =4ne(2kT/(xm))* 


for m small compared to the number of atoms per 
cm*. V is the voltage drop across the contact 
barrier." The quantity ¢ is called here the con- 


* See also F. Seitz and S. Pasternack, NDRC Div. 14 
report, University of Pennsylvania. D1—102, June 10, 1942. 

” W. Schottky, Zeits. f. Physik 118, 539 (1942). 

"' The voltage actually applied to the rectifier is V+iR,, 
where R, is the spreading resistance of the “point” contact. 





SILICON CRYSTAL RECTIFIERS 





729 


tact potential difference (abbreviated below by 
c.p.d.) and is defined by Eq. (1). It is the purpose 
of this paper to present some results of an experi- 
mental investigation of ¢. 

Actually it is necessary to take into account 
the effective lowering of the potential barrier at 
the contact because of the image force’ and the 
tunneling of the electrons (or holes) through the 
barrier. Courant” has shown that under those 
conditions Eq. (2) must be written: 


t= joA -exp(—6e| o| /kT) 
-Lexp(6eV/kT)—exp(—(1—@)eV/RT)]. (3) 


6 is a numerical factor, depending on (| ¢| — V), 
which essentially lowers the barrier height, as 
shown in Fig. 2. Figure 3 is a plot of 6] 9| 
calculated by Courant for a typical contact." 

It is now necessary to make several remarks 
about the rectifier characteristic, Eq. (3). First, 
considering Fig. 1 it is easily seen that an n-type 
contact will rectify only if g¢<0 and a p-type 
contact only if ¢>0. If these conditions are not 
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Fic. 2. Effective potential barrier for electrons or holes 
in a rectifying contact. 
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o , 
Fic. 3. Potential and observed potential for image force 
and tunnel effect (Courant). 


2 E, Courant, NDRC Div. 14 report, Cornell University, 
May 17, 1943; Phys. Rev. 69, 684 (1946). 

3 m=9.1<10-** g, n=5.93X10"* cm™, T=290°K, di- 
electric constant = 10. 
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Fic. 4a. Apparatus used to measure c.p.d. 
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Fic. 4b. Crystal rectifier assembly. 


fulfilled, there will be no effective barrier for the 
electrons (or holes) traversing the contact. The 
barrier resistance is then independent of the 
voltage V™ and equal to zero. Second, in general 
the c.p.d. will not be uniform over the entire con- 
tact area. This so-called multi-contact effect has 
been examined by Sachs and others."* The essen- 


TYPICAL PLOT OF 
LOW LEVEL RESISTANCE VERSUS +. 


Oly] *.276 volts 


@|p| +.29 volts. 


2a 26 28 30 32 12403 (2)! 

Fic. 5. Typical plot of low level resistance versus 1/T. 
B—Be doped Si, polished, treated 2 hrs. at 975° C; 0.007” 
W point, ground; 4 mil. spring deflection, crystal tapped. 


4 R. G. Sachs, NDRC Div. 14 report No. 168, June 15. 
1943; V. A. Johnson, R. N. Smith, and H. J. Yearian, 
NDRC Div. 14 report, Purdue University, August 14, 1943. 
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Fic. 6. Variation of observed c.p.d. Si— W as a function 
of W point area. 


tial result is easily seen by considering the contact 
as made up of two spots with different values of 
g and A. These spots will act like two re. 
sistances in parallel and most of the current will 
flow through the spot of low resistance, i.e., low 
|y|, if the spot areas are not too different. This 
means of course that the lowest absolute value of 
the c.p.d. will determine the effective exponent 
e| | /(RT) in Eq. (3). 


MEASUREMENT OF CONTACT POTENTIAL 
DIFFERENCE (c.p.d.) 


In the present work it was proposed to com- 
pare the c.p.d. with the actual work function 
difference of a semiconductor and a metal not in 
contact. The c.p.d. was determined by a previ- 
ously described method,'* which consists in meas- 
uring the zero voltage resistance of the contact" 
as a function of the temperature. Equation (3) 
yields for this resistance: 


av kT 
R= —) -(- ) exp(ee| el /A7). mn 
di/veo \joAe 


In general the linear variation of R with T can 
be neglected, especially since theoretically jy is 
proportional to T°. 

Figure 4a gives a diagrammatic sketch of the 
apparatus. The semiconductor-metal contact was 
mounted in a conventional temperature stable 
cartridge!’ shown in Fig. 4b. The oven was heated 


18 This was kindly pointed out to the author by B, Seria. 

16 W. E. Stephens, B. Serin, and W. E. Meyerhof, Phys 
Rev. 69, 42, 244 (1946). 

17 The spreading resistance of the contact was subtracted 
from the measured resistance whenever it was not neg 
compared to the barrier resistance R. 

18 See W. E. Stephens, Electronics 19, 112 (July 1946) 
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volt: 


Fic. 7. Distribution curve of observed c.p.d. Si—W. 


to a temperature of about 130°C in a few minutes 
and then allowed to cool. 

About 5 millivolts were applied to the crystal 
in both directions, and the average current was 
used to calculate the resistance of the contact, 
thus eliminating the effect of any thermal e.m.f. 
The cooling of the oven was slow enough (ap- 
prox. $ hr. from 130°C to 50°C) to permit a 
practically simultaneous measurement of tem- 
perature and resistance. Figure 5 shows a typical 
semi-log plot of R versus 1/T. The high tempera- 
ture end of the curve indicates that the semi- 
conductor had not reached the temperature of 
the crystal base at the time of measurement. 


RESULTS OF CONTACT POTENTIAL 
DIFFERENCE MEASUREMENTS 


In order to study the variation of the c.p.d. 
with the nature of the metallic component of the 
contact, a particular piece of p-type silicon'® was 
mounted in a cartridge and different wires were 
pressed against it with various pressures. The 
wires were previously sharpened on Arkansas 
stone to a point of less than 10~° cm radius so 
that different contact pressures produced differ- 
ent contact areas. Figure 6 shows the variation of 
the c.p.d. silicon-tungsten with contact area, 
measured microscopically. 

For most metals used, the c.p.d. was roughly 
independent of area, although there was a trend 
towards an increase in c.p.d. with decreasing area 
especially noticeable for areas less than 0.5 X 10-* 
cm*, This might be caused by a multi-contact 
effect which raises the probability of finding a 
larger minimum (i.e., effective) c.p.d. with small 
contact areas than with large contact areas. 

The unavoidable electrical non-uniformity of 


99.95 percent pure Dupont silicon, 0.002 percent 
boron, 0.02 percent beryllium; dry polished on 0000 emery 
paper; heat treated 2 hours at 975°C (see ref. 18). 
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Silicon ground with No. 600 carborandum. 





Silicon polished on 0000 emery paper. 








int pressed on 
polished silicon. 


Tungsten point pressed on Tungsten 


ground silicon. 


Fic. 8. Photomicrographs of silicon surfaces and tungsten 
points (1000 magnification). 


the silicon surface itself is brought out by Fig. 7 
which shows the distribution curve of approxi- 
mately 50 measurements of the c.p.d. silicon- 
tungsten using the same piece of silicon (2 mm? 
area). For comparison, Fig. 8 gives an idea of the 
geometrical non-uniformity of different silicon 
surfaces and the corresponding tungsten points. 

Table I summarizes the c.p.d. using various 
metals. The first column gives the measured 
c.p.d. 6g and the second column the c.p.d. ¢ cor- 
rected by means of Courant’s curve (Fig. 3). 
Each mean value of the c.p.d. (heavy type) is the 
average of at least 8 measurements. The extreme 
values of the c.p.d. measurements are also given. 

Table II presents some comparison measure- 
ments on a piece of p-type and a piece of n-type 















TaBLeE I, Contact potential differences of various 
metals with respect to p-type silicon. 
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TABLE II. Comparison of contact potential differences 
p- and n-type silicon with various metals, d 















=—=—=—: 
one. p-type peel .— - wetee ee c ~ p-type ee 
Metal 6¢ (volts) @ (volts) = : cae ec «atom Orrected » (volts) 
t .165, 0.20 0.305, 0. 
Al 0.00, 0.01, 0.09 0.00, 0.07,0.21 Ta 045 0:20 a3 ons 
Au 0.21, 0.23, 0.28 0.36, 0.39,044 ww 0.16, 0.18, 0.20 0.30, 0.325, 0,35 
Cu 0.15, 0.22, 0.36 0.29, 0.38, 0.53 c.p.d. n-type silicon—metal ° 
Fe 0.14, 0.17, 0.22 0.28, 0.32, 0.38 Pt —0.16, —0.16 —0.30, —0,30 
Mo 0.17, 0.26, 0.39 0.32, 0.42,0.57 Ta —0.09, —0.095 —0.22, —0.225 
Ni 0.16, 0.31, 0.33 0.30, 0.48, 0.50 W —0,.07, —0.09, —0.10 —0.19, —0,22 —0.23 
Ni-alloy”® 0.21, 0.26, 0.41 0.36, 0.42, 0.59 ’ 
Pt 0.10, 0.17, 0.20 0.22, 0.32, 0.35 
— 0.17, 0.19, 0.37 0.32, 0.34, 0.55 
a 0.10, 0.18, 0.26 0.22, 0.33, 0.42 microwave power (i.e., “burnt out’’). 
W 0.15, 0.23, 0.32 0.29, 0.39, 0.49 power ( ). Only two 








silicon treated as nearly alike as possible,” but 
differently from the piece of p-type silicon used 
in Table I. Only two or three measurements of 
the c.p.d. were made with each metal. 

A comparison of the c.p.d. for the p-type sili- 
con with that in Table I shows that the surface 
structure of the silicon does not seem to affect 
the c.p.d. (the silicon of Table II was ground and 
that of Table I polished—see Fig. 8 for difference 
in structure). 

‘ Both Table I and Table II indicate also that 
the c.p.d. is, within the electrical non-uniformity 
of the surfaces, independent of the nature of the 
metal for the p-type as well as for the n-type 
contacts. This is quite contrary to theoretical 
expectations considering the large variety of 
metals used, which have unoutgassed work func- 
tions between about 3 and 5 volts.¥ In fact the 
c.p.d. values for p-type silicon all lie between 0.6 
and 0.0 volts. Actually the only substance giving 
sometimes the value ¢=0 is aluminum, which 
has one of the lowest work functions (even un- 
outgassed) and which therefore should be ex- 
pected to give one of the highest values of the 
c.p.d. according to Eg. (1). 

Before considering further experiments it is of 
interest to present Table III which gives the 
change in c.p.d. when a crystal rectifier is me- 
chanically tapped and when it is overloaded with 


2096.75 percent Ni, 1.00 percent Al, 2.00 percent Mo, 
0.25 percent C, work hardened wire. 

2190 percent Pt, 10 percent Ru. 

2 For impurities in the p-type silicon see ref. 19; the 
impurities in the n-type piece were not known. Both pieces 
were ground with No. 600 carborundum and water, etched 
together in hot sodium hydroxide, washed and dried. 

% See for example, A. L. Hughes and L. A. DuBridge, 
Photoelectric Phenomena (McGraw-Hill, New York, 1933), 



























series of measurements were made. The last two 
columns in Table III show the barrier resistang 
at room temperature and the relative contag 
area compared to the initial value and calculated 
from Eq. (4) assuming jo constant. It can be seep 
that tapping increases the c.p.d. and the re 
sistance. Burnout decreases the c.p.d., the re 
sistance™ and the effective contact area. 


MEASUREMENT OF WORK FUNCTION 
DIFFERENCE (w.f.d.) 

Since the above mentioned experiments showed 
no correlation of the c.p.d. between silicon and 
various metals, and the work functions of the 
metals, it was thought advisable to measure the 
work function difference (abbreviated below by 
w.f.d.) directly, under various conditions. Al. 
though these measurements eventually did not 
suggest any new phenomena which would make 
the c.p.d. roughly independent of the silicon sur- 
face and the nature of the metal, they demon- 
strate very clearly the disagreement of the c.pd. 
measurements with the presently assumed model 
of the contact. 






























TaBLeE III. Change of contact potential difference ofa 
p-type silicon-tungsten rectifier with tapping and power 
overloading. 
























Relative 
c.p.d. p-type silicon— contact 
t en area 
Treatment Observed Corrected Resistance A/Am 
crystal rectifier 6g (volts) (volts) R (k-ohm) (%) 
Initial values 0.21, 0.17 0.36, 0.32 4.8, 2.4 100, 100 
After tapping 0.28, 0.20 0.44, 0.35 37, 17 200, # 
After 7 watt RF 0.18, 0.14 0.33, 0.28 6.4, 3.8 23, & 
peak power™ 
After 70 watt RF 0.08, 0.10 0.20, 0.23 0.4, 0.9 6, 5 





peak power 















* For a more detailed study of resistance as a function 
of burnout see A. W. Lawson et al., NDRC Div. 14 repot 
No. 113, November 1, 1942. _ 

*6 10 cm microwave microsecond pulses were used in this 
experiment; the power given is available peak power. 
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IV. Variation of work function difference between 
mms p-type silicon and various abraded metals. 
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TABLE VI. Data for the determination of the absolute 
work function of a silicon surface. 








———— 
w.t.d. p-type silicon 





—metal Change Change in Change in 
imin. SOhrs. in@om dom caused ¢om caused 
after after in 50 by adsorbed by adsorbed 

abrasion abrasion hrs.in oxygen water vapor 
of metal of metal air (Dubois) (Dubois) 
Metal ¢o (volts) go (volts) (volts) (volts) (volts) 
Al 0.75 0.21 
Al 0.69 0.26 
Au —0.31 -—045 0.14 0.14 —0.14, —0.21 
Cu —0.37 —0.42 0.05 0.15,0.16 —0.03, —0.06 
Cut —0.34 —0.40 0.06 
Fe —0.18 —0.41 0.23 0.22,0.34 —0.56, —0.56 
Mo —0.17 —0.42 0.25 0.18 —0.16, —0.18 
—0.21 -—0.38 0.17 
Ni —0.09 —0.18 0.11 0.15 —0.16 
—0.08 —0.12 0.04 
Ta 0.13 —0.04 
Pt —0.60 —0.61 0.01 0.04 —0.02 


Ww -0.19 —0.33 











The method of measurement of the w.f.d. was 
essentially a parallel plate condenser method 
(Kelvin method). The two substances under in- 
vestigation form the two plates of a variable con- 
denser. The field established by the w.f.d. (see 
Fig. 1) is brought to zero by an externally ap- 
plied voltage (equal to the w.f.d.). The details of 
the method are described elsewhere.” In opposi- 
tion to the c.p.d. measurement, this method 
gives practically the average w.f.d. between two 
surfaces. 

The surfaces under study were usually abraded 
with 000 emery cloth in order to produce a sur- 
face treatment similar to grinding the metal 
points in the c.p.d. experiments. Also an appa- 


TaBLE V. Variation of work functions of tungsten and 
silicon after abrasion. 








Increase in work function in 
10 min in gas after abrasion 





in vacuum 
Tungsten Silicon 
Gas Initial impurities (volts) (volts) 
Air — 0.33, 0.41 0.10, 0.12 p-type 
0.22, 0.26 n-type 
Oxygen — 0.42,0.46 0.11, 0.15 p-type 
= and — 0.46,0.50 0.16,0.18 “ 
ater 
vapor 
Nitrogen 0.1% O2,0.4% A 0.14, 0.16 0.10, 0.12 
Hydrogen 0.1% Ox 0.08,0.12 0.08, 0.08 
Helium 1.6% Ns, 0.2% hydro- 0.06,0.08 0.04 
carbons 
Argon 0.5% Nz, 0.04,0.06 0.03, 0.05 


0,001% O2, Hz, H:0 








** All metals guaranteed pure except the starred metals 
which are commercial stock. The purity of the metal does 
not influence the work function to any appreciable extent 
in these experiments. 

7 W. E. Meyerhof and P. H. Miller, Jr., Rev. Sci. Inst. 
17, 15 (1946). 





Work 





7 function 
Work function of silicon 
’ - of oxide w.f.d. silicon— surface 
Oxide coating (Lammermann metal (Table IV) oo8 
Metal in air and Lange) ge (volts) (volts) 
Al AlzOs; 3.9% 0.21, 0.26 4.11, 4.16 
Cu Cuz0 4.47 —0.42, —0.40 4.05, 4.07 
Fe FeO or FesO, 4.54 or 4.58 —0.41 4.13 or 4.17 
Pb* PbO 4.06 0.13 4.19 








ratus was constructed which permitted this abra- 
sion in a vacuum of about 5X10-* mm Hg. 


RESULTS OF WORK FUNCTION DIFFERENCE 
MEASUREMENTS 


It is recalled that the w.f.d., go, is defined by 
an equation similar to Eq. (1): 


Yo= dos - dom, 


where ¢os is the work function of the silicon and 
¢gom is the work function of the metal, as meas- 
ured in the present experiment. 

Table IV summarizes the experiments in air 
(on the same metals as used in Table I). A piece 
of polished p-type silicon'® was taken as standard 
surface. Each metal was abraded with 000 emery 
cloth and the w.f.d. measured approximately 1 
min. and 50 hrs. after abrasion. The results on 
Mo and Ni in the above table give an idea of the 
reproducibility of the results. The fourth column 
in Table IV gives the change in the metal work 
function during 50 hours exposure in air. For 
comparison the last two columns show some re- 
sults of Dubois® on the influence of oxygen and 
water vapor on outgassed work functions. It can 
be inferred from these results that the increase in 
work function of the metals after abrasion in air 
is mostly caused by (pure or activated) adsorp- 
tion of oxygen and that the humidity in the air 
does not play any significant role. 

This is also brought out by Table V which gives 
the change in work function of tungsten and sili- 
con (each measured against a standard surface) 
abraded in a vacuum of about 5X10-° mm Hg 
and then subjected to various gases at atmos- 
pheric pressure for 10 minutes. The numbers 


28 Compare with 3.95 volts obtained by H. J. Spanner, 
Ann. d. Physik 75, 608 (1925). 

*® Not shown in Table IV. 

* E. Dubois, Ann. de physique 14, 627 (1930). 
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Fic. 9. Work function of tungsten and tantalum after 
abrasion in partial vacuum. 


given represent averages of two experiments. The 
gases in which oxygen was an impurity were 
passed over red hot copper, so that the percent- 
age oxygen present was certainly less than shown. 
Table V would indicate then that nitrogen and 
hydrogen can also affect the work function of 
solids. 

It is known that most metals oxidize in air. 
Lammermann and Lange® have measured the 
work functions of several oxides and their values 
allow a calibration of the standard silicon surface 
used in Table IV, if we assume that the metals 
in the present experiment were coated with the 
oxides shown below, 50 hrs. after abrasion in air. 
From Table VI a work function between 4.1 and 
4.2 volts may be obtained for the standard silicon 
surface. 

Figures 9 and 10 show the variation of work 
function of tungsten, tantalum, n- and p-type 
silicon after abrasion in a vacuum and then sub- 
jecting the surfaces to air at atmospheric pres- 
sure within 4 min. The measured w.f.d. was 
calibrated by using the standard silicon surface 
and calling its work function 4.2 volts. The in- 


%1 Since nitrogen is generally an impurity in argon this 
might invalidate some of the conclusions which O. Klein 
and E. Lange, Zeits. f. Elektrochemie 44, 542 (1938), drew 
from a change of the work function with time in an atmos- 
phere of argon. 

# Values given by O. Klein and E. Lange (reference 31). 
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Fic. 10. Work function of p- and n-type silicon, 
after abrasion in partial vacuum. 


accuracy of this value will evidently not affect 
the relative scales of Figs. 9 and 10. These figures 
show that for each solid the work function in- 
creases on exposure to air. 


DISCUSSION OF EXPERIMENTAL RESULTS 


Several comments on these measurements are 
now in order. First, assuming the effect of abra- 
sion to be the same on the surfaces of both n- 
and p-type silicon, the difference in work fune- 
tions of the two solids should be equal to the 
width of the forbidden region in silicon (1.12 
volts*). Figure 10 shows that it is only about 0.3 
volts. Second, since the work functions of the 
solids examined change very markedly on ex- 
posure to air, one might assume that on forming 
a contact between a metal and silicon, there oc- 
curs a change in the oxide layers in such a way 
as to make the results of Tables | and II con- 
sistent with those of Table IV and Figs. 9 and 10. 

Holm* has found that in the case of metal- 
metal contacts the plastic deformation of the 
metal parts actually in contact causes a breaking 
of the oxidation layers. Let us assume a similar 
simple mechanism in the formation of crystal 
rectifier contacts. Since the silicon is too hard to 


% F, Seitz, NDRC Div. 14 report No. 113, November 1, 
1942. 

*“R. Holm, Die technische Physik der elektrischen Kon- 
takte (Julius Springer, Berlin, 1941; ‘ton by J. W. 
Edwards, Ann Arbor, 1944), pp. 63-139. 





—_ - +. —_— = Ps mms fF © = fs - oe pee 


no fF «€F oo -« FA 6D A. 





affect 


pures 
Nn in- 








be deformed, we may assume that its oxidation 
layer is conserved during contact formation, so 
that a crystal rectifier contact is perhaps roughly 
a silicon ‘‘oxide”’-metal contact. Table VII com- 
pares the w.f.d. (assuming this kind of contact 
and assuming a silicon “oxide’’-metal ‘‘oxide”’ 
contact) with the c.p.d. (Table II) for the same 
metals. (The work functions of the p- and n-type 
silicon were measured after treating the surfaces 
similarly to those of Table II* and were found 
to be 4.13 and 3.95 volts respectively.) It is easily 
seen that the w.f.d. measurements cannot be 
made consistent with the c.p.d. measurements. 
In this connection it should be noted that the 
multi-contact theory requires for a p-type con- 
tact ¢=0, if go<0, and for an n-type contact 
¢=0, if gor 0. 

The general nature of the results is in essential 
disagreement with some experiments on large- 
area contact rectifiers performed by Joffe.* Joffe 
finds that there is agreement between the sign of 
the w.f.d. and the finite or zero value of the corre- 
sponding c.p.d. It might not be justified, though, 
to compare those results (on semiconducting 
oxides mostly) with the present experiments (on 
silicon). 


CONCLUSION 


Assuming the correctness of the experimental 
results, it seems desirable to reinvestigate the 
physical picture of the silicon-metal contact. Es- 
sentially it is necessary to construct a model 
which makes the c.p.d. roughly independent of 
the surface structure of the silicon and roughly 
independent of the work function of the metal. 
It is possible to construct such a model by as- 
suming the surface of the silicon to be practically 


* A. V. Joffe, J. Phys. USSR 10, 49 (1946). 
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TaBLE VII. Comparison of values for contact potential 
difference and work function difference. 








w.f.d. for icular model 





of contact 
; c.p.d. Silicon “‘oxide”- Silicon “‘oxide’’- 
Nature of (Table II) metal metal “oxide” 
contact ¢ (volts) go (volts) ge (volts) 
p-Si—Pt 0.3 — —0.7 
p-Si—Ta 0.3 0.4 —0.1 
p-Si—W 0.3 0.1 —0.4 
n-Si—Pt —0.3 —_ —0.9 
n-Si—Ta —0,.2 0.2 —0.3 
n-Si—W —0.2 —0.1 —0.6 








“metallic” ** so that the potential barrier would 
be located inside the silicon surface. Any metal 
making contact to such a surface would introduce 
no further discontinuities in the potential, just as 
in the case of metal-metal contacts, so that the 
c.p.d. would be independent of the work function 
of the metal. 

The thickness of the “metallic” layer would 
not have to be more than a few atomic distances 
in order to satisfy the electronic conditions, so 
that this model essentially assumes the existence 
of particular surface states*’ for the conduction 
electrons. Calculations in this general direction 
might throw further light on the mechanism of 
silicon point-contact rectifiers.* 
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** This was kindly pointed out to the author by G. H. 
Wannier. 

87See for example, F. Seitz, The Theory of Solids, 
(McGraw-Hill Book Company, New York, 1940), sec. 70. 

* Note added in proof: Since this paper was sent in for 

ublication, J. Bardeen has kindly informed the author of 
Fis calculations on surface states, which appear in Phys. 
Rev. 71, 717 (1947). 
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The Scattering of Neutrons of Energy Between 
12 Mev and 13 Mev by Protons 


J. S. LAUGHLIN AND P. G. KRUGER 
University of Illinois, Urbana, Illinois 
April 7, 1947 


HE angular distribution of the recoil protons in the 
neutron-proton scattering process has been investi- 
gated in a high pressure (25 atmospheres) cloud chamber. 
The recoil protons were produced in methane gas by a well- 
collimated homogeneous beam of neutrons one-half inch in 
diameter taken in the forward direction from the D(d, )He* 
reaction. A thin (0.5-Mev) deuterium gas target was em- 
ployed and a magnetic field deflected the deuteron beam 
from the forward direction after passing through the target. 
This prevented neutrons from sources other than the deu- 
terium target from entering the cloud chamber. The neu- 
tron spectrum of this arrangement was examined with the 
cloud chamber and showed a strong peak between 12.0 
Mev and 13.0 Mev, superimposed on a low intensity con- 
tinuous spectrum originating from the thin-foil windows on 
the entrance and exit of the target. 

A total of 2000 recoil protons have been measured, of 
which 1636 were produced by neutrons between 12-Mev 
and 13-Mev energy and satisfied selection criteria. The 
data are given in Table I. 

The aximuthal correction factors given were determined 
from the experimental data and are in good agreement with 
the correction factors calculated from the geometry and 
the illuminated height of the chamber. For each scattering 
angle interval the data were plotted in 20° azimuthal angu- 
lar intervals. The shape of this distribution was in good 
agreement with the calculated shape based on purely geo- 


TABLE I. Experimental data for recoil protons produced by 
neutrons of energy between 12 Mev and 13 Mev. 








No. of 
recoil pro- 
tons per 


Azimuthal Length Solid 
unit solid 


Angular No. of correc- _correc- angle 
interval recoil tion tion interval angle 
c.m. system protons factors factors X1/2e X2r/10 


452+10.5 
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Fic. 1. Angular distribution of recoil protons per unit solid angle 
in the center of mass system. 


metrical considerations. The number of recoil protons 
scattered into those azimuthal angle intervals in which all 
tracks were observable determined the. nymber which 
would have been measured if all the tracks in all the azi. 
muthal angle intervals could have been measured. The 
ratio of this number to the actual observed number yielded 
the azimuthal correction factors. 

Since the neutron beam is collimated, all recoil protons 
originate from a 4-inch diameter cylinder through the center 
of the chamber, and since the light beam height is known, 
a calculation of the expected azimuthal correction factor js 
quite straightforward. 

A length correction factor is necessary since for each 
different scattering angle interval the recoil protons havea 
different energy and length. Since these lengths and the 
chamber geometry were well known, these factors could 
be calculated from simple geometrical considerations. 

The graph in Fig. 1 shows the number of recoil protons 
per unit solid angle per scattering angle interval in the 
center of mass system plotted versus the scattering angle in 
the center of mass system. Probable statistical errors are 
indicated in each interval. A ratio of the differential cross 
sections for neutron scattering in the backward direction to 
scattering in the perpendicular direction in the center of 
mass system of either unity (spherically symmetric scatter- 
ing) or slightly greater than unity can be consistent with 
these data. This possible small asymmetry is in the opposite 
direction to that found by Amaldi! and co-workers and by 
Champion and Powell.* Later results by Powell® apparently 
show spherically symmetric scattering. A complete discus- 
sion of these data and the analysis will be published later. 

1 E. Amaldi, D. Bocciarelli, B. Ferretti, and G. C. Trabacchi, Netur- 
wiss. 30, 582 (1942). 


F. C. Champion and C. F. Powell, Proc. Roy. Soc. A183, 64 (1944). 
§ Reference 38, G. Wentzel, Rev. Mod. Phys. 19, 10 (1947). 





Erratum: The K X-Ray Absorption 
Edge of Silicon 


[Phys. Rev. 71, 406 (1947)] 
Vota P. BARTON, Goucher College, Baltimore, Maryland 
AND 
GeorGE A. LINDSAY, Universily of Michigan, Ann Arbor, Michigan 


T was erroneously stated in the above paper that in 
Table I the wave-lengths were corrected only for normal 
refraction. Actually the terms for the anomalous effect had 
also been applied for quartz and the four micas, which are 
the only ones needing such a correction. The anomalous 
terms increased the observed wave-lengths by 0.6 x.u. for 
quartz and by 3.C x.u. for the micas. 
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Gerenkov Radiation Counter* 


R. H. Dicke 


Palmer Physiess Laboratory, Princeton University, 
Princeton, New Jersey 


April 23, 1947 


ETTING! has proposed a high speed fast-particle 
counter using, to actuate a photo-multiplier tube, 
Cerenkov radiation* produced as a charged particle trav- 
erses a plastic rod. The advantages of such a system for fast 
counting, timing, and velocity measurement are evident. 
Furry has recently proposed a measurement of meson mass 
using a modification of the velocity measurement feature.’ 
A device similar to Getting’s Fig. 3 was tried at Princeton 
last summer. The apparatus consisted of a Geiger counter, 
Lucite flared rod, lens system, and photo-multiplier, all in 
a vertical line. A fast charged particle would actuate the 
Geiger counter, then produce in the Lucite Cerenkov radia- 
tion which would in turn be detected by the photo- 
multiplier and associated amplifiers. The apparatus was 
run for 4500 hours between August 1, 1946, and October 18, 
1946, with amplifier discrimination set at various levels; 
coincidence counts definitely attributable to Cerenkov ra- 
diation were not obtained. 

Recently another attempt to demonstrate operation of 
the device was made using electrons generated in lead by 
x-rays from a 20-Mev betatron. In Fig. 1 is shown the 
arrangement of components. The x-rays strike a }-inch 
piece of lead; electrons and positrons generated therein 
pass down a collimating hole in a lead plug and then into 
the Lucite rod. Cerenkov radiation from the rod is focused 
onto the 1P28 photo-multiplier. The latter operates at 
liquid nitrogen temperature in an evacuated envelope not 
shown. Photo-multiplier pulses are fed from a cathode 
follower into a long transmission line, thence into a Los 
Alamos Model 501 amplifier. Pulses from the amplifier are 
fed into a scaler with a calibrated pulse height discrimina- 
tor. The betatron, run at reduced intensity, was monitored 
by a Geiger counter. 

Figure 2 shows the result of one set of measurements. 
Curve A represents the number of pulses greater than a 
certain height plotted against this height. Each point was 
taken for a definite Geiger-monitored x-ray flux. Curve B 
is obtained from A by subtracting dark pulses observed 
with the betatron off. Curve C is obtained from B by sub- 
tracting the x-ray background counts measured by inter- 
posing a thin piece of blackened cardboard between lens 
and photo-tube (Fig. 1). Curve C thus represents the dis- 
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Fic, 1, Cerenkov radiation counter. 
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Fic. 2. A—total pulses. B—total pulses minus dark pulses. 
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tribution of photon pulses. As a check on “pile up” a run 
was made with betatron intensity reduced by a factor of 
2.5; the resulting photon distribution is represented by the 
circles in Fig. 2. 

Curve C is in agreement with the energy distribution of 
electrons produced in lead by bremsstrahlung from 20-Mev 
electrons. Comparison of curve C with the distribution of 
pulses from single photons (unfortunately not taken at the 
same time) indicates probably no more than 5 photo- 
electrons in the largest pulses. This is about } of the com- 
puted? maximum (for photoelectric efficiency see RCA 
Tube Manual). 

In conclusion it can be stated that photon pulses are 
obtained from the counter shown in Fig. 1. The data are in 
reasonable agreement with the hypothesis that the photon 
counts are caused by Cerenkov radiation. Considering the 
geometry of the counter it seems unlikely that fluorescence 
is an important source of photons. However, fluorescence 
cannot be definitely ruled out as a contributing factor. 

The author wishes to express his appreciation to Colonel 
Carl Dutton and Miss Kathleen George of the Picatinny 
Arsenal for their hospitality and helpful assistance in mak- 
ing available the Kerst betatron at Picatinny. He is in- 
debted to Dr. W. C. Elmore for aid with electronic prob- 
lems, and to Dr. D. W. Kerst for advice and many helpful 
suggestions. 

* This work was supported in part by Navy Contract N6OR1-105 
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11. A. Getting, Phys. Rev. 71, 123 (1947). 

?P. A. Cerenkov, Comptes Rendus U.S.S.R. 8, 451 (1934). I. Frank 
and I. Tamm, AL = endus U.S.S.R. 14, 109 (1937). G. B. Collins 
we a hys. Rev. 54, 499 (1938). 

. Furry, Bull. Am. Phys. Soc. 22, No. 2, p. 27 (1947), to be 
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Nuclear Magnetic Resonances at 
Low Temperatures* 


F. Bitter, N. L. Atpert, H. L. Poss, C. G. Lenr, 
AND S. T. Lin** 


Research Laboratory of Electronics, Massachusetts Institute 
of Technology, Cambridge, Massachusetts 


April 14, 1947 


E have recently completed a preliminary set of ex- 
periments on nuclear magnetic resonances at tem- 
peratures down to 4°K, using the 94-inch gap of the M.I.T. 
cyclotron magnet which was recently made available to us 
for a period of several days. Following is a summary of the 
results of these experiments and of a series of room tem- 
perature experiments done with a small electromagnet. 

Proton or fluorine resonances have been observed in all 
substances investigated. Among these have been: benzene, 
ether, paraffin, lithium hydroxide, lithium fluoride, and 
hydrofluoric acid at room temperature; glycerine down to 
195°K (solid CO2); ammonium chloride, water and ice both 
as pure H,0 and in the form of a mixture containing 95 
percent D,O+5 percent H,O down to 77°K (liquid nitro- 
gen); liquid hydrogen at its boiling point; and solid methane 
down to liquid helium temperature. 

In general, broad resonance lines, many gauss wide, were 
found in solids, with the exception of CH,, while lines of 
the order of one gauss and less in width were observed in 
liquids. In reference to the distilled water and heavy water 
solutions, the DO was found to have no appreciable effect 
on the width of the resonance line in both the liquid and 
frozen states. The glycerine resonance line broadened from 
a width of the order of one gauss to a width of about 16 
gauss as the temperature was reduced to 195°K. The LiF 
gave the broadest resonance observed, greater than 25 
gauss wide at room temperature. The resonance line in 
liquid hydrogen was slightly broader than for protons in 
water at room temperature. In most cases, narrow lines 
were observed on an oscilloscope. Because of the lower 
amplitude of broad resonances, these were investigated on 
the output meter of a narrow-band 30-cycle amplifier in 
order to obtain a better signal to noise ratio. 

Our observations of narrow lines in liquids and broad 
lines in solids are in agreement with the theory of line 
widths outlined by Bloembergen, Pound, and Purcell.' In 
this theory the rotational states in liquids are the key to 
narrow lines. In solids the internuclear fields serve to 
broaden the observed resonance line. However, when rota- 
tions occur with frequencies greater than the Larmor 
frequency involved (about 30 Mc for these experiments), 
the internuclear fields practically average out to zero during 
one Larmor cycle, thus resulting in a narrow line. 

Since rotational states generally lead to narrow lines, it 
was decided to study the possibility of rotational states 
above the A-points of ammonium chloride and solid meth- 
ane, which have specific heat anomalies at —30°C and 
20°K, respectively. Although some slight broadening of the 
proton line occurred in cooling ammonium chloride from 
room temperature to 77°K, it was a consistently broad line, 
such as is characteristic of normal solids. This indicates 
that the transition across the A-point is due to an order- 
disorder phenomenon, in agreement with conclusions drawn 
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by Lawson? from an analysis of thermodynamic data. In 
contrast to this are the observations of proton resonances 
in solid methane. In a liquid helium bath a broad line, 
about 12 gauss wide, was found. When the sample was 
allowed to warm up through the A-point, the broad line 
changed to a narrow line, presumably in the Vicinity of 
the A-point. At liquid hydrogen and nitrogen tempera- 
tures, narrow lines, somewhat less than 1 gauss wide, were 
observed in solid methane. 

It should be noted that those molecules with H atoms 
showing A-point transitions are among the easiest to study 
by this method but present considerable difficulty in studies 
by means of x-rays. More detailed investigations of rota. 
tional states in solids are now being undertaken. 

Larger signals were observed at low temperatures, ap- 
proximately following the inverse temperature variation in 
Curie’s formula for the equilibrium magnetization. This is 
in agreement with results obtained at hydrogen tempera. 
tures by Rollin and Hatton* and indicates that, in the 
substances examined, relaxation times do not undergo 
drastic variations at low temperatures. In the case of 
resonances in liquid hydrogen, signal-to-noise ratios of from 
300 to over 1000 were obtained on an oscilloscope without 
the use of a narrow-band amplifier. These results indicate 
that low temperatures may be useful in the study of smal] 
samples or of substances with low concentrations of the 
nucleus under observation, particularly if narrow lines can 
be retained as in methane. However, in most substances, 
broad resonance lines are to be expected, and this will tend 
to reduce the gain caused by the low temperatures. If the 
structure of resonance lines is to be studied, or if the 
natural width of narrow lines is to be measured, stringent 
requirements must be placed on the uniformity of the 
magnetic field over the sample. Otherwise the line width 
and shape will be determined by the magnet and not by 
the sample. 

We are indebted to Professor E. M. Purcell of Harvard 
for his suggestions in reference to the r.f. circuit which we 
are using. We also wish to thank many members of the 
M.I.T. staff, particularly Professor C. F. Squire for his 
assistance in the low temperature experiments, and Dr. 
E. T. Clarke, through whose cooperation we were able to 
use the cyclotron magnet. 


* This research has been supported in part through the Joint Service 
Contract, No. W-36-039 PC-32037. ; 

** On leave from the Institute of Physics, Academia Sinica, China. 

1N. Bloembergen, R. V. Pound, and E. M. Purcell, Phys. Rev. 71, 
466(A) (1947). 

2 Rollin and Hatton, Nature 156, 201 (1947). 

* Andrew W. Lawson, Phys. Rev. 57, 417 (1940). 





The Unboundedness of Quantum-Mechanical 
Matrices 


AUREL WINTNER 


Department of Mathematics, The Johns Hopkins University, 
Baltimore, Maryland 
April 9, 1947 


ILBERT’S definition of a bounded infinite matrix A is 
equivalent to the requirement that the linear substi- 


tution y= Ax should define a vector y= (41, ¥2, « - +) of finite 
length whenever x = (x1, x2, ---) is a vector of finite length 
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data. In (ie., a point of Hilbert’s space). Direct inspections show 
sOnances that all those quantum-mechanical matrices Q, P which 
ad line, have been calculated explicitly (linear oscillator, hydrogen 
Ple was atom, etc.) fail to be bounded in Hilbert’s sense. 
oad line It is natural to expect that this failure is independent of 
inity of the explicit choice of the particular Hamiltonians H in 
empera- question. In fact, the unboundedness of Q, P is indicated 
de, were by heuristic considerations concerning the location of 
possible energy levels. However, the literature consulted 
1 atoms does not contain a general proof. The purpose of this note 
© study is the construction of such a proof. 
| Studies If Planck’s constant is expressed in appropriate units, 
of rota- Heisenberg’s relation can be written in the form 
res, ap. QP—PQ=HE, (1) 
ation in where E denotes the unit matrix and Q, P are supposed to 
This is be Hermitian. Without any reference to a Hamiltonian H, 
"mpera- it will be shown that these assumptions prevent the bound- 
in the edness of Q and P. 
indergo Let A* denote the transposed matrix of the complex- 
case of conjugate of an arbitrary matrix A. Thus A is Hermitian 
of from if and only if A*=A. Since Q and P are Hermitian, it 
without follows that A* is Q—dP, if A denotes Q+iP. Hence, 
— Eq. (1) can be written in the form 
rp A*A—AA*=E. (2) 
nes Cp Furthermore, if P and Q are bounded, then A is bounded. 
tances, It will be shown that these assumptions lead to a con- 
ll tend tradiction. 
? If the By the spectrum of a bounded matrix A is meant the set 
if the of those real or complex values \ for which the bounded 
rungent matrix XE—A fails to have a (unique) bounded reciprocal 
of the (AE—A)~. It is known that the spectrum of every bounded 
width matrix is contained in a sufficiently large circle about the 
not by origin of the \-plane, and that every bounded matrix has 
a spectrum! (i.e., that the spectrum contains at least one 
arvard \-value). Because of the first of these two properties, 
ich we A—}E must have a bounded reciprocal matrix if the abso- 
of the ) lute value of \ is large enough. On the other hand, it is 
or his clear that the expression to the left side of Eq. (2) remains 
id Dr. unaltered if A is replaced by A—AE, where ) is any real 
ble to value. Hence, if A —XE, where ) is real and large enough, 
is denoted by A, there results a bounded matrix A satisfy- 
Service ing Eq. (2) and possessing a bounded reciprocal A~. But 
Chia. Eq. (2) can then be written in the form 
ev. 71, 
B-ABA“™=E, (3) 


where B= A*A. In fact, the product of bounded matrices is 
bounded, and their multiplication is associative. In par- 

. ticular, B is bounded. 

ical It is clear from the above definition of a spectrum that 
if A, A and B are bounded, then B and its transform 
ABA™ have the same spectrum. Since Eq. (3) means that 
ABA~ js identical with B—E, it follows that B and B—E 
have the same spectrum. Finally, the definition of a spec- 
trum also shows that the spectrum of B—E results from 


x A is that of B by a rigid motion, representing a translation (of 
absti- unit span) of the \-plan into itself. 
finite Accordingly, the spectrum of B remains invariant under 


ength a certain translation. Consequently, either the spectrum of 
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B contains every point of a line in the \-plane or the 
spectrum of B has no point at all. However, as mentioned 
above,’ both of these possibilities are ruled out by the 
boundedness of B. This contradiction proves the unbound- 
edness of B= A*A, and therefore that of A =Q+éP. 


_1 A. Wintner, Spektraltheorie der Unendlichen Matrisen (Hirzel, Leip- 
zig, 1929) pp. 143-44. 





On the Dependence of Nuclear Radius on 
the Mass Number 


E. AMALDI AND B. N. Cacctapuort 


Istituto di Fisica della Universita, Centro di Fisica Nucleare 
del C.N.R., Rome, Italy 


April 17, 1947 


N an earlier paper' we reported some measurements on 

nuclear radius for 19 elements with Z between 4 and 

83, and sought a linear representation of our results by an 
expression of the form 


R=b+r,A. (1) 


By graphical interpolation, in which we intentionally 
enhanced the importance of results corresponding to ele- 
ments with high atomic number, we obtained for 6 and ro 
the following values: 


b=1.5 to2.0X10-"% cm, ro=1.3 to 1.0XK10-"% cm. (2) 


The quantity 5 represents the surface effect of the nu- 
cleus, that is, the contribution of nuclear force range of 
superficial particles to the value of R; in fact our measure- 
ments refer to a nucléar radius R defined as the distance 
from the center of a nucleus at which a fast neutron inter- 
acts with the nucleus. Other definitions may be given; for 
instance,? when the binding energies of contiguous and 
specular isobar pairs are compared, the nuclear radius can 
be defined as the radius of the proton-distribution sphere. 
In all cases in which a comparison was possible, the values 
obtained using this second definition were found to be 
smaller than ours by about 1 X 10-" cm. 

In our former paper’ it was mentioned that, notwith- 
standing the accuracy of the values which we found for R, 
it is not possible to determine the values of b and ro with 
great precision, owing to the appreciable deviations of the 
individual points from the linear representation of Eq. (1). 

Therefore we thought it worth while to treat the question 
again by statistical methods in order to test the correctness 
of fit of a theoretical curve to a set of experimental values, 
although such procedures are generally open to some 
criticism. 

The second column of Table I gives the values of R;+4R,; 
for 17 elements. The results obtained for C and Au have 
been omitted; the geometrical conditions for the C measure- 
ments were not quite the same as for the other elements, 
and the results obtained with Au were much less accurate 
than for other elements. 

The third column gives the ratio of the statistical error 
to the standard deviation; the values 5R; of the second 
column are in all cases the greater of the two errors. A 
detailed account of the calculation of these two errors was 
given in reference 1; we shall only stress the fact that in 


























TABLE I. 








Statistical error 





Element (Ri+é6Ri) -10-" cm Standard deviation 
4 Be 3.22 +0.074 1.5 
5B 4.30 +0.24 — 

12 Mg $.40+0.13 0.55 
13 Al 5.54+0.11 0.63 
16S 5.02 +0.14 2 
26 Fe 6.63 +0.06 1,2 
28 Ni 6.47 +0.074 1.5 
29 Cu 6.75 +0.15 0.3 
30 Zn 6.96 +0.17 3 
34 Se 7.31+0.20 0.55 
47 Ag 7.81+1.09 1.3 
48 Cd 8.23 +0.07 2.3 
50 Sn 8.50 +0.085 4 

51 Sb 8.33 +0.096 10 
80 Hg 9.50 +0.16 0.58 
82 Pb 9.00 +0.071 0.87 
83 Bi 9.08 +0.097 3.7 








about half of the cases we found that this ratio was smaller, 
and in the other half was greater than unity, so we believe 
that the main error is purely statistical. This means that 
the probability for the real value to deviate by an amount 
AR; from the observed value may be represented by means 
of a normal distribution (Gaussian law). 

Under this reasonable assumption, we may calculate the 
probability that a given distribution will deviate from the 
interpolated straight line as much, or more, than our 
observed points. 

Such a probability is obviously given by 


re x2 n—l © txt nl 
P= fox dx/ fre xn—Idy, (3) 


where 


r= 3, (AR; /6R:? (4) 


and v is the number of measured points used for comparison 
with the interpolated straight line of Eq. (1); in our case 
n=17. The deviations of the observed points from the 
interpolated line are indicated with AR; and the corre- 
sponding experimental errors are 5R; (see the second column 
of Table I). 

The function defined by (3) is formally identical to the 
function P introduced by Pearson for his ‘x? test,’ 
although in our case the definition of x? is slightly different 
from his. In order to calculate P, the most convenient linear 
representation (1) is the one that makes x? a minimum. 
According to the definition (4) of x?, this obviously means 
that the constants 5 and ro of (1) must be determined by 
the least-squares method, considering the different preci- 
sion 5R; of the single points. 

In this way we get 


b= (0.696 +0.082) x 10- cm, 
ro= (1.52+0.02) X 10-4 cm, (5) 


that is, two values quite different from those obtained with 
the graphical interpolation (2), in which we intentionally 
over-rated the points corresponding to heavy elements. 
Applying definition (4), with the experimental values of 
Table I, and through the representations (1) and (5), we 
obtain for x? the value x*270 which corresponds, for 
n=17, to a very small value of P defined by (3). From 
Pearson's tables,’ it is seen that, for »=17, P is smaller 
than 10-*, even for x?=70. 
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Thus we must argue that representation (1) is not a good 
one if we assume purely accidental errors; in fact, the prob. 
ability of obtaining a distribution such as the one we found 
for the deviations of the experimental values from a linear 
law (1) is exceedingly small. However, it is well known that 
a small value of P does not necessarily indicate that the 
experimental deviations from law (1) are systematic, be. 
cause the “goodness of fit’ test has meaning only in the 
case of purely accidental errors. 

We believe that the very small value obtained by us for 
P gives good evidence for a dependence of the nuclear 
radius on the number of particles in a nucleus, which js 
more complex than the simple law (1). 

 E. Amaldi, D. Bocciarelli, B. N. Cacciapuoti, and G. C. Trabacchi, 
Nuovo. Cimento, in press. Presented at the conference on elementary 
particles and low temperature at Cambridge, England, July, 1946, 

2H. A. Bethe, Phys. Rev. 54, 436 (1938); E. P Wigner, Phys. Rey. 


56, 519 (1939). 
3K. Pearson, Tables for Statisticians and Biometricians, Part I. 





Search for Element 87 


Horta HULUBEI 
Bucarest, Rumania 
March 31, 1947 


HE article by F. R. Hirsh, Jr.' “The Search for 
Element 87” came to my notice only recently because 
of the interruption of postal relations between America and 
Rumania during the war. I wish to comment on the nu- 
merical data employed in his research and on the research 
itself. He adopts 0.812A and 1.045A for the presumptive 
wave-lengths of the absorption line L;;; and the emission 
line La, respectively. The former departs by 0.011A from 
the value deduced, by Moseley’s law, from the most reliable 
experimental values obtained for neighboring elements. 
The latter departs by 0.017A from the value to be expected 
for Le;, and by 0.011A from the average of the to-be- 
expected values of La; and Laz. Convincing evidence for 
the element 87 could not be found, even assuming the 
element to be present in the sample, with such expectations 
as to the wave-lengths. I suggest that his categorical objec- 
tions to the inferences of other authors are based on very 
inexact estimates of the emission-lines which are to be 
expected. 

I find the experiments of Hirsh far from convincing. He 
used a tungsten filament and two molybdenum anticath- 
odes, and had to clean the anticathodes frequently to 
remove the pulverized tungsten from them. A photograph 
taken with an exposure of 36 hours, and having given no 
indication whatever of the tungsten L lines under the con- 
ditions of his observations, should not be expected to reveal 
the presence of an element present very probably in much 
smaller quantities than the tungsten remaining on the anti- 
cathode even after hourly cleanings. In addition, the Ly 
lines of tungsten would mask the presumptive line La: of 
element 87 if they were both emitted. Similarly, the lines 
of molybdenum should have made their appearance: 
whether or not they were sufficiently excited by the narrow 
residual band of continuous x-rays springing from the pri- 
mary anticathode under a 25-kv bombardment and falling 
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on the secondary anticathode, they should have appeared 
as scattered radiation. What is one to think of the scatter- 
ing with change of frequency considered by Hirsh and 
s tly the object of an inexact formula in his text? 
Myself, I should not have melted the sample even “gently”’ 
in order to insert it into the hole in the anticathode, in view 
of the presumptive properties of element 87. 

The search for a stable isotope of element 87 in a mixture 
of heavy, volatile, and absorbing alkaline substances is 
dificult, even if the element should be present in propor- 
tions of 10 or 10-*. The methods of spectrum analysis and 
the conditions of excitation must be the best attainable. 
Our observations were generally made by primary excita- 
tion, with tubes constructed specially for great stability of 
operation even when volatile substances are present. My 
spectra do not allow me to come to the same conclusion as 
that reached by Dr. Hirsh. For the reasons expounded 
above, I feel that his investigation does not help in the 
least to settle the question of the existence of a durable 
isotope of element 87. 


iF, R. Hirsh, Jr., Phys. Rev. 63, 93 (1943). 





The Magnetic Quenching of Superconductivity 


J. W. Stout 


Institute for the Study of Metals, the University of Chicago, 
Chicago, Illinois 


April 12, 1947 


N a letter with the above title Sienko and Ogg' suggest 
that for the “‘soft’’ superconductors Pb, Hg, Sn, In, 
Tl, CuS, AusBi, Zn, and Cd the threshold magnetic field, 
Hr, for the destruction of superconductivity, is better 
represented by the expression Hr = A(T.!— T) than by the 
parabolic relation Hr = B(T2—T*). T. is the temperature 
at which the metal becomes superconducting in zero field. 
Careful measurements of the threshold magnetic fields have 
been made on Sn**5, Hg**, In*5, TI*® and Pb**. Within 
the accuracy of the measurements neither the parabolic nor 
the #-power law exactly fits the data. However, the para- 
bolic law is in all cases the better approximation. The 
variation in the constants A (}-power law) and B (para- 
bolic law) between about 1°K, and the transition tempera- 
ture are: Sn (A 30 percent; B 10 percent), Hg (A 18 per- 
cent; B 3 percent), In (A 12 percent; B 8 percent), Tl (A 9 
percent; B 5 percent), Pb (A 19 percent; B 8 percent). The 
values quoted are the averages obtained by the various 
experimenters. In general the variations among the meas- 
urements from different laboratories are less than the 
deviation from either the parabolic or $-power relation. 

It is true that the constant A is more nearly the same for 
different elements than is B, but A as determined from the 
initial slope of Hr vs. T varies from 90 for Cd® to 46 for Sn, 
a variation well outside the experimental error and suffi- 
cient to make it doubtful that A is a combination of uni- 
versal constants as Sienko and Ogg suggest. 

From an expression for the threshold curve and the well- 
known thermodynamic relation 


C.—C,=(VT/4r){ (dHr/dT+Hr(@Hr/dT*)}, 


one can calculate the difference in heat capacity, C,— Cn, 


LETTERS TO 








THE EDITOR 741 


between the superconducting and normal states. V is the 
volume. For the }-power relation one obtains 


C.—C, = (VA*/4e)(3T?— 37.479), 
and for the parabolic relation 
C.—C, = (VB*/2r)(3T*—T2T). 


Measurements on Sn’ and TI® indicate that the heat ca- 
pacity of these two metals in the superconducting state 
obeys a T* law. For metals in the non-superconducting 
state the experimental heat capacities at low temperatures 
are in agreement with the theoretical prediction of an elec- 
tronic heat capacity proportional to the temperature plus 
a lattice heat capacity which in simple cases varies as 7°. 
The parabolic relation predicts a difference in heat capacity 
of superconducting and normal metal in agreement with 
these considerations. The absence of a term linear in tem- 
perature from C,—C, obtained from the §-power law is 
inconsistent with the experimental and theoretical results 
for the heat capacity of metals in the non-superconducting 
state. 


1M. J. Sienko and R. A. . Phys. Rev. 71, 319 (1947). 

2W. J. de Haas and A. D. Engelkes, Physica 4, 325 (1937). 

3A. D. Misener, Proc. Roy. Soc. A174, 262 (1940). 

4 J. G. Daunt and K. Mendelssohn, Proc. Roy. Soc. A160, 127 (1937). 
(1938)" Daunt, A. Horseman, and K. Mendelssohn, Phil. Mag. 27, 754 

*N. Kurti and F. Simon, Proc. Roy. Soc. A151, 610 (1935). 

7W. H. Keesom and P. H. Van Laer, Physica 5, 193 (1938). 

*W. H. Keesom and J. A. Kok, Physica 1, 175 (1933). 












Bounded Linear Harmonic Oscillator and 
Phase Transitions of Second Order 


B. SuRYAN 
Bangalore City, South India 
March 28, 1947 


N a recent letter in the Physical Review, Corson and 
Kaplan! have put forth the idea that the concept of a 
linear bounded harmonic oscillator in the theory of solids 
may be of help in understanding anomalies in the specific 
heat-curve of a crystalline solid. Similar ideas were de- 
veloped by the author in 1945 and were reported in a lec- 
ture at Central College, Bangalore City, South India. These 
ideas were not published at the time, but in view of the 
considerations of Corson and Kaplan, it has seemed desir- 
able to place them on record. 

The problem arose in connection \with specific heat 
anomalies at low temperatures in crystalline solids. It was 
thought that the concept of the linear bounded harmonic 
oscillator might be of service in studying these phenomena. 
The relevant theory has been worked out by Auluck and 
Kothari.? The allowed energy values are no longer (m-+-})hyv, 
where » is the frequency of the oscillator, but increase 
rapidly with decreasing amplitudes of oscillation. However, 
the present author applied the above concepts to the 
second-order phase transition or \-anomalies of Ehrenfest.* 
It is well known that such phase transitions are all ac- 
companied by marked changes in the expansion coefficient 
of the solid. Therefore, if at the \-temperature the different 
oscillators in the crystal are bounded to different extents, 
the oscillators will be distributed among two well-defined 








sets of energy states, and these two sets may be thought of 
as forming the “‘phases’’ above mentioned. This may con- 
veniently be taken as a picture of the second-order phase 
transition, i.e., the crystal oscillators going from one state 
of bounding to another. Also because the quantum of 
transition betyveen two states of the oscillator is no longer 
hv, a shift or a new line may be possible in the Raman spec- 
trum or infra-red absorption of the crystal near the )- 
temperature. As there are many substances whose \-tem- 
peratures are near room temperature, the investigation of 
the Raman spectrum may provide a test of the above ideas. 
Also, one must consider the current theory that the specific 
heat change may arise from molecular rotation within the 
crystal. Experiments on the Raman or infra-red absorption 
spectrum may lead to a decision between the two con- 
ceptions. 
E. M. Corson and I. Kaolin. Phys. Rev. 71, 130 (1947). 

atk S. Auluck and D. S. Kothari, Proc. Camb. Phil. Soc. 41, 175 


sP. Ehrenfest, Proc. K. Acad. Weten. Amsterdam 36, 153 (1933); 
Comm. No. 75b from the K. Onnes Laboratory, Leyden. 





On the Distribution of Electric 
Charge in Nuclei 


B. Davison AND W. H. WATSON 
National Research Council, Chalk River, Ontario, Canada 
April 21, 1947 


N a recent letter! by J. G. Winans entitled ‘A Classical 
Model for the Nucleus,” attention was drawn to the 
striking relation between the atomic numbers (Z) and 
Corresponding masses (A) of stable nuclei. Winans sug- 
gested a crystalline model in which the protons occupy sur- 
face positions as an obvious way of explaining the simple 
relation between Z and A. 

The existence of this relation and its possible explanation 
by surface concentration of protons was already known to 
us, and we had gone so far as to test the hypothesis that 
protons are concentrated on or near the surface in the 
following simple way. In fission, if protons and neutrons 
are uniformly mixed the (charge/mass) ratios of the frag- 
ments should be roughly equal; whereas we should expect 
a higher relative concentration of charge on the light frag- 
ments, if the protons are concentrated near the surface of 
the nucleus before fission. We have investigated this on the 
basis of Coryell’s tables? and have found that the ratio 
(charge/mass) is very closely the same for the light and 
heavy fragments. This fact, therefore, coupled with the 
hydrodynamical explanation’ of the asymmetry of fission 
in analogy with the splitting of a classical liquid drop, 
seems to invalidate the idea that the protons are concen- 
trated on the surface. In principle it is possible that some 
mechanism might exist to ensure constancy of Z/A for the 
fragments, but we know of no other nuclear phenomenon 
which would require such a mechanism. 


1J. G. Winans, Phys. Rev. 71, 379 (1947). 
2 J. A. Chem. Soc. 68, 2411 (1946). 
3 J. A. Wheeler, unpublished lecture notes, August, 1946. 
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A Note on the Long-Lived Radio-Iodine 


L. E. GLENDENIN AND R. R. Epwarps 
Massachusetts Institute of Technology,* Cambridge, Massachusetts 
April 9, 1947 


es radioactive iodine isotope of half-life about 56 
days, reported by Reid and Keston,! has been ob. 
served in an iodine sample separated from deuteron. 
bombarded Te in April, 1946, furnished by the M.LT. 
Radioactivity Center. 

Chemical identity of the activity was established bya 
series of extractions with iodine carrier. Iodide in aqueous 
solution was oxidized with nitrite, and the iodine extracted 
into CCl; re-extraction to water was performed with hj. 
sulfite as reducing agent. Through three such cycles, the 
specific activity as AgI remained constant within experi. 
mental error. 

By use of a Kr-filled Geiger counter, critical absorption 
measurements of the x-radiations were made in Ag, Cd, In, 
and Sn foils. The energy of the radiations was observed to 
lie between the K-edges of Cd and In, characteristic of 
tellurium Ka x-rays. Absorption in Al showed a small 
amount of an extremely soft component, and a half. 
thickness of 485 mg/cm? for the major fraction of the 
radiations. This corresponds to an x-ray energy of 275 
kev, within experimental error of that expected for Ke 
x-rays of Te. The soft component showed a half-thickness 
of 1.9 mg/cm? in Al, and 19 mg/cm? in C (polystyrene), 
The corresponding x-ray energies are 3.8 and 3.7 key, re. 
spectively. The La x-ray energy for Te is 3.8 kev. No 
particle radiations were observable from 20 mg samples 
(AgI or Nal) with a thin window counter, and no y-rays 
of energy greater than 30 kev were apparent in the Al 
absorption curve. The mode of decay of the 56d I, therefore, 
appears to be pure orbital electron capture, with no ac- 
companying particle or y-radiations. 

It can be shown by various considerations that the mass 
number of the long-lived iodine is more probably 125 than 
129 as proposed by Reid and Keston.' The formation of 
I'25 directly in the Te targets at the cyclotron is to be 
expected from the d, 2 reaction on Te™® (6 percent abund- 
ance) and the (d, m) reaction on Te™ (4.5 percent abund- 
ance). Considerations which render improbable the mass 
assignment of 129 are: (a) the proposed decay scheme! 
requires either that Te! have a stable isomer, or (for the 
chain to reach a stable end product) that a unique chain of 
K-capture followed by two 8 decays occur (i.e., I™—K-= 
Te!®—¢-—+'®—_8--+Xe"™ stable); (b) according to Bohr 
and Wheeler? the most stable nuclear charge (Z.) for mass 
129 is near 54 (Xe); thus I should have no appreciable 
probability for decay by K-capture, but should decay by 
B~ emission to stable Xe™; (c) no radio-iodine of haif-life 
near 56 days is reported among the U™* fission products! 
whereas I! must certainly be formed; the latter is listed 
as a very long-lived, hypothetical descendant of 70m Te*. 

The nucleus I'*5, however, would certainly decay by K- 
capture or positron emission to stable Te!® and would not 
be formed in appreciable yield as a U** fission product since 
it is shielded by stable Te5 from formation by §~ decay. 
Its predicted decay energy for K-capture, calculated from 






































the equation of Bohr and Wheeler,’ is about 0.2 Mev. This 
is concordant with the fact that no energetic radiations are 
observed. 

Possible mass assignments‘ other than 125 may be tenta- 
tively ruled out. On the basis of the Bohr-Wheeler? picture, 
[8 and isomeric states of the nuclei of mass 124, 126, 128, 
and 130 would be expected to have high decay energies 
(>1 Mev). A radioactive isomer of stable I’ would decay 
by isomeric transition, and thus would yield no Te x-rays. 
Iodine isotopes of mass number greater than 130 have 
been conclusively characterized as energetic 8~ emitters.’ 

The yield of the 56d I in deuteron-bombarded Te targets 
at the M.I.T. cyclotron is about 15yc/u amp. hr. of 14- 
Mev deuterons, calculated on an assumption of one K 
x-ray per disintegration, and an experimentally determined 
counting efficiency of about 1 percent in the usual argon- 
filled Geiger counter. The corresponding yield of 12h I™ is 
750 pc/u amp. hr. (for short bombardments). The yield 
ratio of 12h I to 56d I at saturation is 4.5. The abundance 
ratio of Te™ to Te, from which I and I are formed 
by the (d, 2m) reaction, is 5.5. Assuming the (d, m) reaction 
on Te™ to be less probable than (d, 2m) on Te™®, the pro- 
posed mass assignment of 125 for the long-lived I is con- 
sistent with the observed yield. 

We wish to thank Professor M. Deutsch of the Physics 
Department and Professors C. D. Coryell and J. W. Irvine, 
Jr., of the Chemistry Department for valuable suggestions 
and criticisms, and the Radioactivity Center for supplying 
the sample of radioactive iodine. 

* Contribution from the Eugectment of Chemistry and the Labora- 
tory for Nuclear Science and Engineering; eapoerted in part by U. S. 
Navy tment, contract NSori78 Task VI. 
1A. F. Reid and A. S. Keston, Phys. Rev. 70, 987 (1946). 

?N. Bohr and J. A. Wheeler, Phys. Rev. 56, 426 (1939). 


“Nuclei Formed in Fission,"” J. Am. Chem. Soc. 68, 2411 (1946). 
4G. T. Seaborg, Rev. Mod. Phys. 16, 1 (1944). 





55-Hour Element 61 Formed in Fission* 

Mark G. INGHRAM, Davin C. Hess, Jr., AND RICHARD 

J. HAYDEN, Argonne National Laboratory, Chicago, Illinois 
AND 


GrorGE W. PARKER, Monsanto Chemical Company, 
Clinton Laboratory, Oak Ridge, Tennessee 


April 15, 1947 


CHEMICAL separation of element 61 from fission 
products was made to verify the mass assignment and 
to study the activity of 49-hr. 61, as summarized in fission 
tables.! This activity has been assumed to have a mass of 
149 from fission yield data. 
An aliquot of the element 61 sample was placed on the 
filament of a mass spectrograph, and normal neodymium 
was added to serve as a mass standard. By operation of the 
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spectrograph the isotopes were separated according to mass 
and deposited on a photographic plate. A second “transfer” 
photographic plate was placed face to face with the first 
plate, and four developable images were obtained from its 
radioactive isotopes which by comparison with the normal 
neodymium spectrum on the first plate were shown to have 
masses 147, 149, 163, and 165. The metal lines at 147 and 
149 were very weak in comparison with the metal oxide 
lines at 163 and 165 but appeared in approximately the 
same ratio. Since the ratios of metal to metal oxide are 
different for the different elements which might have given 
these masses, this indicates that all lines are due to the 
same element. 

To verify that the activity at the mass 149 position was 
49-hr. 61' the following technique was used. The photo- 
graphic plate containing the active isotopes was placed 
successively against various parts of a larger photographic 
plate for times calculated to give equal intensity at the 
mass 149 position if the half-life were 49 hours. Five ex- 
posures were taken for successive times of 13.4 hours, 16.5 
hours, 21.6 hours, 31.0 hours, and 57.3 hours. These trans- 
fers showed the 61’ to be decaying with a half-life of 
slightly more than 49 hours, while the 147 showed no de- 
tectable decay. A decay curve taken simultaneously on 
another aliquot of the element 61 sample showed it to 
contain a 55-hour activity and a long-lived activity. Thus 
the 147 and 163 lines are due to the metal and metal oxide 
ions of 3.7-year 617 as shown previously by Hayden,* and 
the lines at 149 and 165 are due to the metal and metal 
oxide ions of 55-hr. 61”. 

The samarium-europium fraction of the original fission 
sample was also analyzed by this method. It showed ac- 
tivities at masses 151, 153, 155, 156, 167, and 169. The 
active lines at mass numbers 153 and 169 were found to 
decay with a half-life of 47 hours. Since the ratio of intensi- 
ties of the 153 to 169 is that characteristic of Sm* and 
SmO? this activity is 47-hr. Sm'™* as was shown previously 
by neutron irradiation of samarium.* The lines at 151 and 
167 were also ascribed to samarium since the ratio of the 
intensities was the same as the 153 and 169. This is the 
activity first found in fission by Lewis and Hayden. From 
intensity considerations the half-life of this activity has 
been estimated at roughly 20 years. The lines at 155 and 156 
were assigned to europium since only the metal ion lines 
appeared. The 156 decayed with a half-life of about 15.4 
days and the 155 showed no detectable decay. Thus these 
activities are due to 2 to 3 year Eu'®* and 15.4-day Eu'®*. 


* This document is based on work performed under Contract No. 
W-7401-eng-37 for the Manhattan Project at the Argonne National 
Laboratory. 

1 Rev. Mod. Phys. 18, 513 (1946). 

? Richard J. Hayden, Phys. Rev. (1947), to be published. 

* Richard J. Hayden and Mark G. Inghram, Phys. Rev. 70, 89 (1946). 
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MINUTES OF THE MEETING OF THE NEW YorRK STATE SECTION AT 
ROCHESTER UNIVERSITY, APRIL 12, 1947 


HE Spring Meeting of the New York State 
Section of the American Physical Society 
was held at Rochester University on Saturday, 
April 12, 1947. An attendance of over 130 mem- 
bers and visitors was recorded. 
The following program of invited papers was 
presented : 


Nuclear Physics and Machines. S. N. VAN VooruHIs, 
University of Rochester. 


The 130 Cyclotron at Rochester. S. W. BARNEs, University 


of Rochester. 


Photographic Plates for Use in Nuclear Physics. J. H, 
Wess, Eastman Kodak Company. 

Separation of Stable Tracer Isotopes. D. W. Srewarr, 
Eastman Kodak Company. 

Biological Effects of Nuclear Radiation. W. F. Baz, 
University of Rochester. 

Plans of the Brookhaven National Laboratory. G. B. 
Coxiitns, University of Rochester. 


The Association of Rochester Scientists, The 
Executive Committee Rochester Association. 


C. L. ANDREWs, 
Secretary 
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